Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “model 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,369 records · Page 76

On the Implementation and Further Validation of a Time Domain Boundary Element Method Broadband Impedance Boundary Condition

A time domain boundary integral equation with Burton-Miller reformulation is presented for acoustic scattering by surfaces with liners in a uniform mean flow. The Ingard-Myers impedance boundary condition is implemented using a broadband multipole impedance model and converted into time domain differential equations to augment the boundary integral equation. The coupled integral-differential equations are solved numerically by a March-On-in-Time (MOT) scheme. While the Ingard-Myers condition is known to support Kelvin-Helmholtz instability due to its use of a vortex sheet interface between the flow and the liner surface, it is found that by neglecting a second derivative term in the current time domain impedance boundary condition formulation, the instability an be effectively suppressed in computation. The proposed formulation and implementation are validated using a NASA Langley Research Center Grazing Flow Impedance Tube (GFIT) experimental dataset with satisfactory results. Moreover, a minimization procedure for finding the poles and coefficients of the broadband multiple impedance model is formulated in this paper by which, unlike the commonly used vector-fitting method, passivity of the model is ensured. Numerical tests show the proposed minimization approach is effective for modeling liners that are commonly used in aeroacoustic applications.

Time Domain Impedance Boundary Condition↗

Meteorological adjustment of yearly mean values for air pollutant concentration comparison

Using multiple linear regression analysis, models which estimate mean concentrations of Total Suspended Particulate (TSP), sulfur dioxide, and nitrogen dioxide as a function of several meteorologic variables, two rough economic indicators, and a simple trend in time are studied. Meteorologic data were obtained and do not include inversion heights. The goodness of fit of the estimated models is partially reflected by the squared coefficient of multiple correlation which indicates that, at the various sampling stations, the models accounted for about 23 to 47 percent of the total variance of the observed TSP concentrations. If the resulting model equations are used in place of simple overall means of the observed concentrations, there is about a 20 percent improvement in either: (1) predicting mean concentrations for specified meteorological conditions; or (2) adjusting successive yearly averages to allow for comparisons devoid of meteorological effects. An application to source identification is presented using regression coefficients of wind velocity predictor variables.

Sidik, S. M.↗

RXTE Observation of Cygnus X-1: III. Implications for Compton Corona and ADAF Models: Implications for Compton Corona and ADAF Models - Report 3

We have recently shown that a 'sphere + disk' geometry Compton corona model provides a good description of Rossi X-ray Timing Explorer (RXTE) observations of the hard/low state of Cygnus X-1. Separately, we have analyzed the temporal data provided by RXTE. In this paper we consider the implications of this timing analysis for our best-fit 'sphere + disk' Comptonization models. We focus our attention on the observed Fourier frequency-dependent time delays between hard and soft photons. We consider whether the observed time delays are: created in the disk but are merely reprocessed by the corona; created by differences between the hard and soft photon diffusion times in coronae with extremely large radii; or are due to 'propagation' of disturbances through the corona. We find that the time delays are most likely created directly within the corona; however, it is currently uncertain which specific model is the most likely explanation. Models that posit a large coronal radius [or equivalently, a large Advection Dominated Accretion Flow (ADAF) region] do not fully address all the details of the observed spectrum. The Compton corona models that do address the full spectrum do not contain dynamical information. We show, however, that simple phenomenological propagation models for the observed time delays for these latter models imply extremely slow characteristic propagation speeds within the coronal region.

Nowak, Michael A.↗

Numerical grid generation; Proceedings of the Symposium on Numerical Generation of Curvilinear Coordinate Systems and Their Use in the Numerical Solution of Partial Differential Equations, Nashville, TN, April 13-16, 1982

General curvilinear coordinate systems are considered along with the error induced by coordinate systems, basic differential models for coordinate generation, elliptic grid generation, conformal grid generation, algebraic grid generation, orthogonal grid generation, patched coordinate systems, and solid mechanics applications of boundary fitted coordinate systems. Attention is given to coordinate system control and adaptive meshes, the application of body conforming curvilinear grids for finite difference solution of external flow, the use of solution adaptive grids in solving partial differential equations, adaptive gridding for finite difference solutions to heat and mass transfer problems, and the application of curvilinear coordinate generation techniques to the computation of internal flows. Other topics explored are related to the solution of nonlinear water wave problems using boundary-fitted coordinate systems, the numerical modeling of estuarine hydrodynamics on a boundary-fitted coordinate system, and conformal grid generation for multielement airfoils.

Thompson, J. F.↗

Simultaneous X-ray and Far-Ultraviolet Spectra of AGN with ASCA and HUT

We obtained ASCA spectra of the Seyfert 1 galaxy NGC 3516 in March 1995. Simultaneous far-UV observations were obtained with the Hopkins Ultraviolet Telescope on the Astro-2 shuttle mission. The ASCA spectrum shows a lightly absorbed power law of energy index 0.78. The low energy absorbing column is significantly less than previously seen. Prominent 0 VII and 0 VIII absorption edges are visible, but, consistent with the much lower total absorbing column, no Fe K absorption edge is detectable. A weak, narrow Fe K(alpha) emission line from cold material is present as well as a broad Fe K(alpha) line. These features are similar to those reported in other Seyfert 1 galaxies. A single warm absorber model provides only an imperfect description of the low energy absorption. In addition to a highly ionized absorber with ionization parameter U = 1.66 and a total column density of 1.4 x 10(exp 22)/sq cm, adding a lower ionization absorber with U = 0.32 and a total column of 6.9 x 10(exp 21)/sq cm significantly improves the fit. The contribution of resonant line scattering to our warm absorber models limits the Doppler parameter to less than 160 km/s at 90% confidence. Turbulence at the sound speed of the photoionized gas provides the best fit. None of the warm absorber models fit to the X-ray spectrum can match the observed equivalent widths of all the UV absorption lines. Accounting for the X-ray and UV absorption simultaneously requires an absorbing region with a broad range of ionization parameters and column densities.

Kriss, Gerard A.↗

Spectral and Temporal Properties of the Ultra-Luminous X-Ray Pulsar in M82 from 15 Years of Chandra Observations and Analysis of the Pulsed Emission Using NuSTAR

The recent discovery by Bachetti et al. of a pulsar in M82 that can reach luminosities of up to 10(exp 40) erg s(exp -1), a factor of approximately 100 times the Eddington luminosity for a 1.4 solar mass compact object, poses a challenge for accretion physics. In order to better understand the nature of this source and its duty cycle, and in light of several physical models that have been subsequently published, we conduct a spectral and temporal analysis of the 0.58 keV X-ray emission from this source from 15 years of Chandra observations. We analyze 19 ACIS observations where the point-spread function (PSF) of the pulsar is not contaminated by nearby sources. We fit the Chandra spectra of the pulsar with a power-law model and a disk blackbody model, subjected to interstellar absorption in M82. We carefully assess for the effect of pile-up in our observations, where four observations have a pile-up fraction of 10, which we account for during spectral modeling with a convolution model. When fitted with a power-law model, the average photon index when the source is at high luminosity (LX greater than 10(exp 39) erg s(exp -1) is equal to gamma 1.33 +/-.0.15. For the disk blackbody model, the average temperature is T(sub in) 3.24 +/- 0.65 keV, the spectral shape being consistent with other luminous X-ray pulsars. We also investigated the inclusion of a soft excess component and spectral break, finding that the spectra are also consistent with these features common to luminous X-ray pulsars. In addition, we present spectral analysis from NuSTAR over the 3-50 keV range where we have isolated the pulsed component. We find that the pulsed emission in this band is best fit by a power-law with a high-energy cutoff, where gamma is equal to 0.6 +/- 0.3 and E(sub C) is equal to 14(exp +5) (sub -3)) keV. While the pulsar has previously been identified as a transient, we find from our longer-baseline study that it has been remarkably active over the 15-year period, where for 9/19 (47%) observations that we analyzed, the pulsar appears to be emitting at a luminosity in excess of 10( exp (39) erg s (exp −1), greater than 10 times its Eddington limit.

individual (M82) - stars: neutron X rays; binaries↗

An empirical propellant response function for combustion stability predictions

An empirical response function model was developed for ammonium perchlorate propellants to supplant T-burner testing at the preliminary design stage. The model was developed by fitting a limited T-burner data base, in terms of oxidizer size and concentration, to an analytical two parameter response function expression. Multiple peaks are predicted, but the primary effect is of a single peak for most formulations, with notable bulges for the various AP size fractions. The model was extended to velocity coupling with the assumption that dynamic response was controlled primarily by the solid phase described by the two parameter model. The magnitude of velocity coupling was then scaled using an erosive burning law. Routine use of the model for stability predictions on a number of propulsion units indicates that the model tends to overpredict propellant response. It is concluded that the model represents a generally conservative prediction tool, suited especially for the preliminary design stage when T-burner data may not be readily available. The model work included development of a rigorous summation technique for pseudopropellant properties and of a concept for modeling ordered packing of particulates.

Hessler, R. O.↗

Some closing thoughts: Practical payoffs from satellite systems

The benefits-to-cost ratio of satellite remote sensing, both as a substitute for conventional methods of monitoring and assessing resources, and as a supplement to these methods is examined using a model which analyzes the cost of aerial photography versus satellite scanner for producing and interpreting an image of the Earth's surface sized to LANDSAT dimensions. Examples of cost savings are tabulated for ground surveys, aerial photos, and LANDSAT. Possible additional benefits from LANDSAT D are assessed. The way in which satellites fit into more comprehensive models for resources management is discussed. It is shown that remote sensing is but one essential component in a complex system that aggregates technical. Socioeconomics, political, cultural, and other factors in the human decision process.

Source record↗

A broad-band X-ray telescope observation of the black hole candidate LMC X-1

We present the spectrum of the black hole candidate LMC X-1 as observed by the Broad-Band X-ray Telescope (BBXRT). The spectrum cannot be fitted by a simple model, but it requires a soft disk blackbody component and a power-law tail, confirming earlier studies. The blackbody disk component is essentially unchanged since the Ginga measurement in 1987. The implied mass of the compact object is approximately 4.7 (square root of cosine i) solar mass from the soft component fit. We report a 95% confidence detection of weak emission features at approximately 5.1 keV and approximately 7.3 keV. If the 5.1 keV feature is attributed to Fe I K-alpha at 6.39 keV, then the redshift is approximately 0.19. No quasi-periodic behavior is found in the data at this epoch.

Schlegel, Eric M.↗

GRAIL Refinements to Lunar Seismic Structure

Joint interpretation of disparate geophysical datasets helps to reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co‐located at one corner. Due to this limited geographical extent, near‐surface ray coverage from moonquakes is low, but increases with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefore well suited to joint inversion, since the complementary information reduces inherent model ambiguity. Previous joint inversions of the Apollo seismic data (seismic phase arrival times) and Clementine‐ or Lunar Prospector‐derived gravity data (mass and moment of inertia) attempted to recover the subsurface structure of the Moon by focusing on hypothetical lunar compositions that explore the density/velocity relationship. These efforts typically search for the best fitting thermodynamically calculated velocity/density model, allowing variables like core size, velocity, and/or composition to vary freely. Seismic velocity profiles previously derived from the Apollo seismic data through inversion of travel times vary both in the depth of the crust and mantle layers, and the seismic velocities and densities assigned to those layers. The lunar mass and moment of inertia likewise only constrain gross variations in the density profile beyond that of a uniform density sphere. As a result, composition and structure models previously obtained by jointly inverting these data retain the original uncertainties inherent in the input data sets. We will perform a joint inversion of Apollo seismic delay times and gravity data collected by the GRAIL lunar gravity mission, in order to recover seismic velocities and density as a function of latitude, longitude, and depth within the Moon. We will relate density to seismic velocity using a linear relationship that is allowed to be depth‐dependent. The corresponding coefficient (B) can reflect a variety of material properties that vary with depth, including temperature and composition. The inversion seeks to recover the set of density, velocity, and B‐coefficient perturbations that minimize (in a least‐squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗

The P/Halley: Spatial distribution and scale lengths for C2, CN, NH2, and H2O

From P/Halley long slit spectroscopic exposures on 12 dates, extending from Oct. 1985 to May 1986, spatial profiles were obtained for emissions by C2, CN, NH2, and OI(1D). Haser model scale lengths were fitted to these data. The extended time coverage allowed the checking for consistency between the various dates. The time varying production rate of P/Halley severely affected the profiles after perihelion, which is shown in two profile sequences on adjacent dates. Because of the time varying production rate, it was not possible to obtain reliable Haser model scale lengths after perihelion. The pre-perihelion analysis yielded Haser model scale lengths of sufficient consistency that they can be used for production rate determinations, whenever it is necessary to extrapolate from observed column densities within finite observing apertures. Results of scale lengths reduced to 1 AU are given and discussed.

Fink, Uwe↗

GRAIL Refinements to Lunar Seismic Structure

Joint interpretation of disparate geophysical datasets helps reduce drawbacks that can result from analyzing them individually. The Apollo seismic network was situated on the lunar nearside surface in a roughly equilateral triangle having sides approximately 1000 km long, with stations 12/14 nearly co-located at one corner. Due to this limited geographical extent, near-surface ray coverage from moonquakes is low, but increases with depth. In comparison, gravity surveys and their resulting gravity anomaly maps have traditionally offered optimal resolution at crustal depths. Gravimetric maps and seismic data sets are therefore well suited to joint inversion, since the complementary information reduces inherent model ambiguity. Previous joint inversions of the Apollo seismic data (seismic phase arrival times) and Clementine- or Lunar Prospector-derived gravity data (mass and moment of inertia) attempted to recover the subsurface structure of the Moon by focusing on hypothetical lunar compositions that explored the density/velocity relationship. These efforts typically searched for the best fitting thermodynamically calculated velocity/density model, and allowed variables like core size, velocity, and/or composition to vary freely. Seismic velocity profiles derived from the Apollo seismic data through travel time inversion vary both in the depth of the crust and mantle layers, and the seismic velocities and densities assigned to those layers. The lunar mass and moment of inertia likewise only constrain gross variations in the density profile beyond that of a uniform density sphere. As a result, composition and structure models previously obtained by jointly inverting these data retain the original uncertainties inherent in the input data sets. We perform a joint inversion of Apollo seismic delay times and gravity data collected by the GRAIL lunar gravity mission, in order to recover seismic velocity and density as a function of latitude, longitude, and depth within the Moon. We relate density (p) to seismic velocity (v) using a depth-dependent linear relationship. The corresponding coefficient (B) can reflect a variety of material properties, including temperature and composition. The inversion seeks to recover the set of p, v, and B perturbations that minimize (in a least-squares sense) the difference between the observed and calculated data.

Weber, Renee C.↗

Application of a new criterion for modeling systems

A new criterion is proposed for modeling systems which promises to be useful in deciding how complex a model should be. The criterion is based on the expected model response error instead of the error in fitting the data used for estimating the model parameters. The new criterion also does not require withholding data to be used exclusively for testing. There remains, however, the difficulty of testing a large number of candidate models that correspond to the combinations of terms used in the dynamic equations. A computational approach is suggested which greatly reduces the computations required in searching for the best model. In the suggested approach the gradient of the response with respect to the model coefficients is held fixed and numerous combinations of terms are assessed. After determining the most promising candidate model, the gradient is updated and the process is repeated. This procedure gives greater assurance that the best model is selected and does not rely on the analyst's judgement.

Taylor, L. W., Jr.↗

Synthesis of rotor test data for real-time simulation

A mathematical model of a hingeless tilting rotor is presented. The model was obtained by a systematic curve fit procedure applied to an extensive set of model scale wind tunnel data. The math model equations were used in a real time flight simulation model of a hingeless tilt rotor XV-15 to assess changes in flying qualities compared to those obtained using a previous rotor model. Extensive plots of the rotor derivatives are given. Discussions of attempts to apply multivariable linear regression technqiues to the data and the use of an analytical rotor representation are included.

Mcveigh, M. A.↗

Non-LTE model calculations for SN 1987A and the extragalactic distance scale

This paper presents model atmospheres for the first week of SN 1987A, based on the luminosity and density/velocity structure from hydrodynamic models of Woosley (1988). The models account for line blanketing, expansion, sphericity, and departures from LTE in hydrogen and helium and differ from previously published efforts because they represent ab initio calculations, i.e., they contain essentially no free parameters. The formation of the UV spectrum is dominated by the effects of line blanketing. In the absorption troughs, the Balmer line profiles were fit well by these models, but the observed emissions are significantly stronger than predicted, perhaps due to clumping. The generally good agreement between the present synthetic spectra and observations provides independent support for the overall accuracy of the hydrodynamic models of Woosley. The question of the accuracy of the Baade-Wesselink method is addressed in a detailed discussion of its approximations. While the application of the standard method produces a distance within an uncertainty of 20 percent in the case of SN 1987A, systematic errors up to a factor of 2 are possible, particularly if the precursor was a red supergiant.

Schmutz, W.↗

Improved passive microwave sounding of the atmosphere

The effort this year focused primarily on 118-GHz transmittance experiments. The data analyzed here was collected with the Microwave Temperature Sounder (MTS) radiometer package during the CAMEX deployment of 1993 with the aim of validating current models of atmospheric microwave absorption in the O2 bands near 54 and 118 GHz. Particular attention has been paid to data collected during four flights when the MTS scanned zenith while profiles of downwelling radiances were collected through ascents and descents. These radiances, in conjunction with radiosonde temperature data, permit the retrieval of band-averaged absorption profiles for each channel. The Millimeter-wave Propagation Model (MPM92) provides theoretical expressions for the absorption of microwaves by oxygen and water vapor and accounts for the interference of pressure-broadened spectral lines'. This model is a good fit to laboratory measurements at temperatures ranging from 279-327 K, but it has been suggested that extrapolation to the conditions of the atmospheric tropopause may result in underestimation of absorption by as much as 15 percent. Preliminary results of the analysis of MTS data appear to be in general agreement with the predictions of the MPM model to within the accuracy of the measurements, which through the coldest parts of the atmosphere ranges from less than plus or minus 5 percent in the most opaque channels to greater than plus or minus 10 percent in the most transparent channels. At those altitudes where each channel is most sensitive to changes in absorption, there is some indication that the modeled absorption may be biased low relative to the observations. Accurate instrument calibration provided challenges, particularly when observed radiances were as much as 260 K below the temperatures of the cold calibration load.

Staelin, D. H.↗

Boundary layer separation and reattachment

A review of developments in a new model for boundary-layer separation is presented. The separation model encompasses timevarying as well as steady-flow conditions. Engineering criteria in terms of calculable pressure gradients are obtained for both laminar and turbulent boundary-layer separation. These criteria, together with limited information on time-dependent separations, should be of value in predicting flows in turbomachinery. The reattachment of boundary layers is also shown to fit within the separation model.

Sandborn, V. A.↗

Temperature Structure in the Perseus Cluster Core Observed with Hitomi

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8-20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region _ 7′ in diameter. The SXS was operated with an energy resolution of _5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components

Aharonian, Felix↗