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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS1989

Application of statistical distribution theory to launch-on-time for space construction logistic support

The ability to launch-on-time and to send payloads into space has progressed dramatically since the days of the earliest missile and space programs. Causes for delay during launch, i.e., unplanned 'holds', are attributable to several sources: weather, range activities, vehicle conditions, human performance, etc. Recent developments in space program, particularly the need for highly reliable logistic support of space construction and the subsequent planned operation of space stations, large unmanned space structures, lunar and Mars bases, and the necessity of providing 'guaranteed' commercial launches have placed increased emphasis on understanding and mastering every aspect of launch vehicle operations. The Center of Space Construction has acquired historical launch vehicle data and is applying these data to the analysis of space launch vehicle logistic support of space construction. This analysis will include development of a better understanding of launch-on-time capability and simulation of required support systems for vehicle assembly and launch which are necessary to support national space program construction schedules. In this paper, the author presents actual launch data on unscheduled 'hold' distributions of various launch vehicles. The data have been supplied by industrial associate companies of the Center for Space Construction. The paper seeks to determine suitable probability models which describe these historical data and that can be used for several purposes such as: inputs to broader simulations of launch vehicle logistic space construction support processes and the determination of which launch operations sources cause the majority of the unscheduled 'holds', and hence to suggest changes which might improve launch-on-time. In particular, the paper investigates the ability of a compound distribution probability model to fit actual data, versus alternative models, and recommends the most productive avenues for future statistical work.

Morgenthaler, George W.↗

NASA NTRS1995

Elliptical accretion disks in active galactic nuclei

We present a calculation of the profiles of emission lines originating in a relativistic, eccentric disk, and show examples of the resulting model profiles. Our calculations are motivated by the fact that in about one-quarter of the double-peaked emission lines observed in radio-loud active galactic nuclei (and in the mildly active nucleus of NGC 1097), the red peak is stronger than the blue peak, which is contrary to the prediction of relativistic, circular disk models. Using the eccentric disk model we fit some of the observed profiles that cannot be fitted with a circular disk model. We propose two possible scenarios for the formation of an eccentric disk in an active galactic nucleus: (a) tidal perturbation of the disk around a supermassive black hole by a smaller binary companion, and (b) formation of an elliptical disk from the debris resulting from the tidal disruption of a star by the central black hole. In the former case we show that the eccentricity can be long-lived because of the presence of the binary companion. In the latter case, although the inner parts of the disk may circularize quickly, we estimate that the outer parts will maintain their eccentricity for times much longer than the local viscous time. We suggest that it may be possible to detect profile variability on much shorter timescales than those ranging from a decade to several centuries by comparing the evolution of the line profile with detailed model predictions. We argue that line-profile variability may also be the most promising discriminant among competing models for the origin of asymmetric, double-peaked emission lines.

Eracleous, Michael↗

NASA NTRS2013

Modeling the HD 32297 Debris Disk With Far-Infrared Herschel Data

HD 32297 is a young A-star (approx. 30 Myr) 112 pc away with a bright edge-on debris disk that has been resolved in scattered light. We observed the HD 32297 debris disk in the far-infrared and sub-millimeter with the Herschel Space Observatory PACS and SPIRE instruments, populating the spectral energy distribution (SED) from 63 to 500 micron..We aimed to determine the composition of dust grains in the HD 32297 disk through SED modeling, using geometrical constraints from the resolved imaging to break the degeneracies inherent in SED modeling. We found the best fitting SED model has two components: an outer ring centered around 110 AU, seen in the scattered light images, and an inner disk near the habitable zone of the star. The outer disk appears to be composed of grains>2 micron consisting of silicates, carbonaceous material, and water ice with an abundance ratio of 1:2:3 respectively and 90% porosity. These grains appear consistent with cometary grains, implying the underlying planetesimal population is dominated by comet-like bodies. We also discuss the 3.7 sigma detection of [C ii] emission at 158 micron with the Herschel PACS instrument, making HD 32297 one of only a handful of debris disks with circumstellar gas detected

infrared: stars↗

NASA NTRS2020

Evaluation of Fatigue Damage Accumulation Functions for Delamination Initiation and Propagation

The present report follows on the cohesive fatigue damage model methodology proposed in NASA-TP-2018-219838. In that report, an empirical function describing the incremental damage due to cyclic loading was used to calculate fatigue damage within a cohesive formulation. The form of the function was developed such that, when integrated at a constant stress amplitude from no damage to failure, it produces a life versus load response that is consistent with an S-N diagram. Therefore, the parameters of the damage model could be obtained by fitting the model predictions to an S-N diagram. The finite element analyses performed demonstrate that the cohesive fatigue accumulation function provides a link between the S-N diagram that describes crack initiation, and the Paris law that characterizes the rate of crack propagation. However, when the model was proposed, it was not known whether the form of the damage accumulation function associated with a desired S-N diagram is unique and, if not, if the link between S-N and the Paris law is unique and independent of the fatigue function selected. In the effort described herein, several alternative forms of the damage function that reproduce the desired features of S-N diagrams were found and evaluated. The effects of each of these functions on the predicted parameters of the Paris law and the propagation threshold are discussed. The results indicate that the predicted exponent m of the Paris law is indeed independent of the damage accumulation function. However, different functions predict different values for the pre-factor C of the Paris law. Therefore, the proper damage accumulation function must be selected by comparison with experiments. One of the new damage accumulation functions proposed herein was found to be particularly useful for analysis because of the ease with which the model parameters can be determined with a minimal amount of experimental information. The effectiveness of the proposed methodology and damage function was demonstrated by conducting analyses of a double cantilever beam test, a mixed-mode bending test, and a three-point bending test of a skin/doubler specimen. The results indicate that the same set of model parameters can provide accurate predictions of the rate of fatigue crack propagation for a variety of material interfaces, mode mixities, load levels, and stress ratios.

Cohesive elements↗

NASA NTRSDate not supplied

Finite Element Modeling of Extravehicular Mobility Units for Use With Human Body Models – Motivation, Major Challenges, Use Cases and Preliminary Work

Finite element modeling of pressurized spacesuits and implementation with human body models offers many advantages over physical experimentation, but also presents significant challenges. A model of pressurized spacesuit softgoods was developed, integrated with an existing hardgoods model, and fitted to a human body model. Two representative loading scenarios were simulated: dynamic external suit loading and internal occupant-driven loading, to serve as proof-of-concept for the modeling techniques employed. The modeled interactions behaved as intended and illustrate that finite element models of pressurized spacesuits can be used in conjunction with human body models to assess the biomechanical behavior. Further model development and experimental validation are needed.

Finite Element↗

NASA NTRS2024

Finite Element Modeling of Extravehicular Mobility Units for Use With Human Body Models – Motivation, Major Challenges, Use Cases and Preliminary Work

Finite Element↗

NASA NTRS1997

Hierarchical Modeling and Robust Synthesis for the Preliminary Design of Large Scale Complex Systems

Large-scale complex systems are characterized by multiple interacting subsystems and the analysis of multiple disciplines. The design and development of such systems inevitably requires the resolution of multiple conflicting objectives. The size of complex systems, however, prohibits the development of comprehensive system models, and thus these systems must be partitioned into their constituent parts. Because simultaneous solution of individual subsystem models is often not manageable iteration is inevitable and often excessive. In this dissertation these issues are addressed through the development of a method for hierarchical robust preliminary design exploration to facilitate concurrent system and subsystem design exploration, for the concurrent generation of robust system and subsystem specifications for the preliminary design of multi-level, multi-objective, large-scale complex systems. This method is developed through the integration and expansion of current design techniques: Hierarchical partitioning and modeling techniques for partitioning large-scale complex systems into more tractable parts, and allowing integration of subproblems for system synthesis; Statistical experimentation and approximation techniques for increasing both the efficiency and the comprehensiveness of preliminary design exploration; and Noise modeling techniques for implementing robust preliminary design when approximate models are employed. Hierarchical partitioning and modeling techniques including intermediate responses, linking variables, and compatibility constraints are incorporated within a hierarchical compromise decision support problem formulation for synthesizing subproblem solutions for a partitioned system. Experimentation and approximation techniques are employed for concurrent investigations and modeling of partitioned subproblems. A modified composite experiment is introduced for fitting better predictive models across the ranges of the factors, and an approach for constructing partitioned response surfaces is developed to reduce the computational expense of experimentation for fitting models in a large number of factors. Noise modeling techniques are compared and recommendations are offered for the implementation of robust design when approximate models are sought. These techniques, approaches, and recommendations are incorporated within the method developed for hierarchical robust preliminary design exploration. This method as well as the associated approaches are illustrated through their application to the preliminary design of a commercial turbofan turbine propulsion system. The case study is developed in collaboration with Allison Engine Company, Rolls Royce Aerospace, and is based on the Allison AE3007 existing engine designed for midsize commercial, regional business jets. For this case study, the turbofan system-level problem is partitioned into engine cycle design and configuration design and a compressor modules integrated for more detailed subsystem-level design exploration, improving system evaluation. The fan and low pressure turbine subsystems are also modeled, but in less detail. Given the defined partitioning, these subproblems are investigated independently and concurrently, and response surface models are constructed to approximate the responses of each. These response models are then incorporated within a commercial turbofan hierarchical compromise decision support problem formulation. Five design scenarios are investigated, and robust solutions are identified. The method and solutions identified are verified by comparison with the AE3007 engine. The solutions obtained are similar to the AE3007 cycle and configuration, but are better with respect to many of the requirements.

Koch, Patrick N.↗

NASA NTRS2005

Spectral Photometric Properties of the Moon

We modeled the solar phase curves of the moon at a series of wavelengths using the full disk telescopic observations [1]. We endeavored to keep the database self-contained, that is, to use the values derived for the solar magnitude and phase curves of the disk-integrated [1]. These observations were made in a suite of 10 narrowband filters between 0.315 microns and 1.06 microns, and in the broad band Johnson UBV filters, as part of a larger program to obtain photoelectric photometry of the larger planets. Two aspects of the lunar observations are unique. First, the observations cover phase angles from 6deg through 120deg. More importantly, the observers used a special 20-mm diameter f/15 fused quartz lens constructed solely for this purpose. The lens reduced the whole lunar image in the focal plane to a size comparable to the planets observed as part of the same program. This image was fed directly into the photometer. Thus, these observations constitute the only existing set of phase curves of the entire lunar disk over a range of wavelengths. Table 1 lists the values of the Hapke model parameters which fit the data. Figure 1 is an example of the model fits to the data.

Dominque, D.↗

NASA NTRS1984

Satellite ephemerides for the Voyager Uranus encounter

Uranian satellite ephemerides are needed by the Voyager project to support both navigation and acquisition of scientific data. This paper presents the approach being taken to develop the ephemerides and details the initial phase of the development. That phase involves the analytical modeling of the satellites' motion and the adjustment of the model to fit astronomical observations. The paper describes the model and gives the result of a fit to 71 years of observations.

Jacobson, R. A.↗

NASA NTRS2019

Winter Wheat Yield Assessment from Landsat 8 and Sentinel-2 Data: Incorporating Surface Reflectance, Through Phenological Fitting, into Regression Yield Models

A combination of Landsat 8 and Sentinel-2 offers a high frequency of observations (3–5 days) at moderate spatial resolution (10–30 m), which is essential for crop yield studies. Existing methods traditionally apply vegetation indices (VIs) that incorporate surface reflectances (SRs) in two or more spectral bands into a single variable, and rarely address the incorporation of SRs into empirical regression models of crop yield. In this work, we address these issues by normalizing satellite data (both VIs and SRs) derived from NASA’s Harmonized Landsat Sentinel-2 (HLS) product, through a phenological fitting. We apply a quadratic function to fit VIs or SRs against accumulated growing degree days (AGDDs), which affects the rate of crop development. The derived phenological metrics for VIs and SRs, namely peak, area under curve (AUC), and fitting coefficients from a quadratic function, were used to build empirical regression winter wheat models at a regional scale in Ukraine for three years, 2016–2018. The best results were achieved for the model with near infrared (NIR) and red spectral bands and derived AUC, constant, linear, and quadratic coefficients of the quadratic model. The best model yielded a root mean square error (RMSE) of 0.201 t/ha (5.4%) and coefficient of determination R2 = 0.73 on cross-validation.

phenological fitting↗

NASA NTRS2003

WMAP Detection of Early Reionization

The Wilkinson Microwave Anisotropy Probe (WMAP) has mapped the full sky in Stokes I, Q, and U parameters at frequencies 23, 33, 41, 61, and 94 GHz. We detect correlations between the temperature and polarization maps significant at more than 10 standard deviations. The correlations are inconsistent with instrument noise and are significantly larger than the upper limits established for potential systenatic errors. The correlations are present in all WMAP frequency bands with similar amplitude from 23 to 94 GHz, and are consistent with a superposition of a CMB signal with a weak foreground. The fitted CMB component is robust against different data combinations and fitting techniques. On small angular scales the WMAP data show the temperature-polarization correlation expected from adiabatic density perturbations. The data for $\ell > 20$ agree well with the signal predicted solely from the temperature power spectra, with no additional free parameters. The existence of correlations on super-horizon scales provides significant support for inflationary cosmologies. We detect excess power on large angular scales compared to predictions based on the temperature power spectra alone. The excess power is well described by reionization at redshift $11 < z < 30$ at 95\% confidence. A model-independent fit to reionization optical depth yields results consistent with the best-fit $\Lambda$CDM model, with best fit value $\tau = 0.17 \pm 0.04$ at 68\% confidence, including systematic and foreground uncertainties. This value is larger than expected given the detection of a Gunn-Peterson trough in the absorption spectra of distant quasars, and implies that the universe has a complex ionization history. I will discuss the WMAP data and its implications for reionization in the early universe.

Kogut, A.↗

NASA NTRS2003

WMAP Detection of Early Reionization

The Wilkinson Microwave Anisotropy Probe (WMAP) has mapped the full sky in Stokes I, Q, and U parameters at frequencies 23, 33, 41, 61, and 94 GHz. We detect correlations between the temperature and polarization maps significant at more than 10 standard deviations. The correlations are inconsistent with instrument noise and are significantly larger than the upper limits established for potential systematic errors. The correlations are present in all WMAP frequency bands with similar amplitude from 23 to 94 GHz, and are consistent with a superposition of a CMB signal with a weak foreground. The fitted CMB component is robust against different data combinations and fitting techniques. On small angular scales the WMAP data show the temperature-polarization correlation expected from adiabatic density perturbations. The data for $\ell > 20$ agree well with the signal predicted solely from the temperature power spectra, with no additional free parameters. The existence of correlations on super-horizon scales provides significant support for inflationary cosmologies. We detect excess power on large angular scales compared to predictions based on the temperature power spectra alone. The excess power is well described by reionization at redshift $11 < z < 30$ at 95\% confidence. A model-independent fit to reionization optical depth yields results consistent with the best-fit $\Lambda$CDM model, with best fit value $\tau = 0.17 \pm 0.04$ at 68\% confidence, including systematic and foreground uncertainties. This trough in the absorption spectra of distant quasars, and implies that the universe has a complex ionization history. I will discuss the WMAP data and its implications for reionization in the early universe.

Kogut, A.↗

NASA NTRS2012

Modeling and Maximum Likelihood Fitting of Gamma-Ray and Radio Light Curves of Millisecond Pulsars Detected with Fermi

Pulsed gamma rays have been detected with the Fermi Large Area Telescope (LAT) from more than 20 millisecond pulsars (MSPs), some of which were discovered in radio observations of bright, unassociated LAT sources. We have fit the radio and gamma-ray light curves of 19 LAT-detected MSPs in the context of geometric, outermagnetospheric emission models assuming the retarded vacuum dipole magnetic field using a Markov chain Monte Carlo maximum likelihood technique. We find that, in many cases, the models are able to reproduce the observed light curves well and provide constraints on the viewing geometries that are in agreement with those from radio polarization measurements. Additionally, for some MSPs we constrain the altitudes of both the gamma-ray and radio emission regions. The best-fit magnetic inclination angles are found to cover a broader range than those of non-recycled gamma-ray pulsars.

Johnson, T. J.↗

NASA NTRS1995

Lithospheric bending at subduction zones based on depth soundings and satellite gravity

A global study of trench flexure was performed by simultaneously modeling 117 bathymetric profiles (original depth soundings) and satellite-derived gravity profiles. A thin, elastic plate flexure model was fit to each bathymetry/gravity profile by minimization of the L(sub 1) norm. The six model parameters were regional depth, regional gravity, trench axis location, flexural wavelength, flexural amplitude, and lithospheric density. A regional tilt parameter was not required after correcting for age-related trend using a new high-resolution age map. Estimates of the density parameter confirm that most outer rises are uncompensated. We find that flexural wavelength is not an accurate estimate of plate thickness because of the high curvatures observed at a majority of trenches. As in previous studies, we find that the gravity data favor a longer-wavelength flexure than the bathymetry data. A joint topography-gravity modeling scheme and fit criteria are used to limit acceptable parameter values to models for which topography and gravity yield consistent results. Even after the elastic thicknesses are converted to mechanical thicknesses using the yield strength envelope model, residual scatter obscures the systematic increase of mechanical thickness with age; perhaps this reflects the combination of uncertainties inherent in estimating flexural wavelength, such as extreme inelastic bending and accumulated thermoelastic stress. The bending moment needed to support the trench and outer rise topography increases by a factor of 10 as lithospheric age increases from 20 to 150 Ma; this reflects the increase in saturation bending moment that the lithosphere can maintain. Using a stiff, dry-olivine rheology, we find that the lithosphere of the GDH1 thermal model (Stein and Stein, 1992) is too hot and thin to maintain the observed bending moments. Moreover, the regional depth seaward of the oldest trenches (approximately 150 Ma) exceeds the GDH1 model depths by about 400 m.

Levitt, Daniel A.↗

NASA NTRS1961

A Study of the Free Oscillations of the Earth

Published observations on the toroidal oscillations of the earth are critically reviewed. A supplementary analysis of the record obtained by the Lamont strain seismometer is presented. Eleven toroidal modes are identified, and it is concluded that the periods are known to within 1 per cent. A perturbation scheme involving the ratio of the angular velocity of the earth to the resonant frequency is used in calculating the effects due to the rotation of the earth on the resonant frequency. The free oscillations are viewed as a superposition of traveling waves. In a nonrotating system two traveling waves combine to produce a stationary standing wave. In a rotating system, the rotation distinguishes between waves that travel in the direction of rotation and those that travel in the opposite direction. Rotation removes a degeneracy and results in a splitting of a spectral peak of order l into 2 times l plus 1 peaks. The fractional displacement in frequency for the lowest-order toroidal oscillations is 1/206 and of the same order as the Q of the peak, so that splitting will probably not be observed in the toroidal oscillations. Viewed locally, rotation causes a particle to precess about a direction parallel to the axis of rotation. This precession will cause a variation of amplitude with time if the motion is recorded by an instrument with an anisotropic response function. Care is therefore needed in studying the time decay of a given spectral peak. Rotation also couples the normal coordinates so that a motion that is initially purely horizontal will develop a vertical component. It is expected that vertical seismometers should record particle motion with the toroidal frequencies. The perturbations of the toroidal oscillations due to core-mantle interaction are treated in detail. An exact expression is obtained for the rate of energy dissipated by a finitely conducting plate oscillating across a magnetic field. The energy dissipated at the core-mantle boundary due to viscous and hydromagnetic coupling is shown to be insignificant as compared with the energy dissipated within the mantle. The toroidal magnetic field leaking into the lower mantle combines with the dipole field, resulting in a stress on the mantle, tending to stiffen the lower boundary. The stress is of sufficient magnitude to produce a displacement toward higher frequency in the lower-order toroidal oscillations. Observations on the (sub 0) T (sub 2) oscillations lead to an estimate of the toroidal magnetic field in the lower mantle. A calculation of elastic energy in the low-order oscillations suggests a value of 10 (sup 18) ergs per cycles per hour for the energy density at low frequencies in the Chilean earthquake. Each mode of oscillation has a characteristic radial distribution of elastic energy associated with it. This distribution determines which parts of the earth contribute most heavily in determining a particular resonant frequency. The distribution of energy for the lower 17 modes for a homogeneous and a Gutenberg model earth is calculated. The resonant frequencies for models of the earth based on the Gutenberg and Lehmann distribution of elastic properties are presented. It is shown that the Gutenberg model earth fits the observations more closely than the Lehmann model and that a slight alteration of the Gutenberg model gives a significantly better fit to the observations. The alteration involves a lower shear-wave velocity in the lower mantle while the Gutenberg velocity distribution is maintained in the upper mantle. Various studies of the earth's oscillations coupled with surface-wave investigations substantiate Gutenberg's hypothesis of a layer of low velocity in the upper mantle. The physical conditions required for the formation of a region of low velocity are examined in detail. The results confirm Birch's earlier statement that a temperature gradient in excess of 6 degrees to 7 degrees per kilometer is needed to produce a decrease in velocity. The low-velocity layer does not require that the temperature approach or exceed the melting temperature. If tile upper mantle is homogeneous, the region of lower velocity should commence at the base of the crust and extend to 150 kilometers under the oceans and about 100 kilometers under continental regions. The distribution of thermal conductivity and radioactivity consistent with the low-velocity layer is also considered.

MacDonald, Gordon J. F.↗

NASA NTRS2016

Foreground Bias from Parametric Models of Far-IR Dust Emission

We use simple toy models of far-IR dust emission to estimate the accuracy to which the polarization of the cosmic microwave background can be recovered using multi-frequency fits, if the parametric form chosen for the fitted dust model differs from the actual dust emission. Commonly used approximations to the far-IR dust spectrum yield CMB residuals comparable to or larger than the sensitivities expected for the next generation of CMB missions, despite fitting the combined CMB plus foreground emission to precision 0.1 percent or better. The Rayleigh-Jeans approximation to the dust spectrum biases the fitted dust spectral index by (Delta)(Beta)(sub d) = 0.2 and the inflationary B-mode amplitude by (Delta)(r) = 0.03. Fitting the dust to a modified blackbody at a single temperature biases the best-fit CMB by (Delta)(r) greater than 0.003 if the true dust spectrum contains multiple temperature components. A 13-parameter model fitting two temperature components reduces this bias by an order of magnitude if the true dust spectrum is in fact a simple superposition of emission at different temperatures, but fails at the level (Delta)(r) = 0.006 for dust whose spectral index varies with frequency. Restricting the observing frequencies to a narrow region near the foreground minimum reduces these biases for some dust spectra but can increase the bias for others. Data at THz frequencies surrounding the peak of the dust emission can mitigate these biases while providing a direct determination of the dust temperature profile.

cosmic microwave background↗

NASA NTRS1983

Laboratory studies at high resolution of the infrared absorption spectra of a number of gases found in planetary atmospheres

The infrared absorption spectra of a number of gases found in planetary atmospheres were studied at high resolution. Absorption line measurements which can be of value for the interpretation of planetary spectra in terms of molecular abundances and conditions in the planetary atmospheres were provided. The high resolution spectra have yielded measurements of individual vibration rotation line parameters including positions, strengths, pressure broadened widths and, where assignments were unknown, the temperature sensitivity of the strengths. Such information allows the determinations of the absorption of a given molecular gas under planetary conditions of temperature and pressure and at the same time it provides the data necessary if the spectra are to be understood in terms of basic molecular theory. Thus this work has included spectral analysis in the form of line assignments as well as fitting of the data to Hamiltonian models. Such fitting is very useful in that it helps to confirm and extend the assignments.

Hunt, R. H.↗

NASA NTRS1984

The structure of clusters of galaxies observed with Einstein

Imaging observations from the Einstein Observatory are used to determine properties of galactic clusters and the intracluster gas. The hydrostatic-isothermal gas model is applied to obtain the cluster core radius, the limiting slope of the X-ray surface brightness distribution, and the presence and magnitude of a central luminosity excess above that predicted from the model when fit to the outer regions of the cluster. It is concluded that clusters can be divided into two basic families, one with small core radii and X-ray emission centered on a central dominant galaxy, and one with larger core radii and generally no central, stationary bright galaxy. The former type probably begins to form in the early stages of a cluster's dynamical evolution. The intracluster gas has a larger scale height than that of the galaxies based on the fit of the hydrostatic-isothermal model to the X-ray surface brightness profiles.

Jones, C.↗