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 937 records · Page 52

Numerical analysis of flow non-uniformity in the hot gas manifold of the Space Shuttle main engine

Three-dimensional viscous flow in a conceptual hot gas manifold (HGM) for the Space Shuttle Main Engine High Pressure Fuel Turbopump (SSME HPFTP) was numerically analyzed. A finite difference scheme was used to solve the Navier-Stokes equations. The exact geometry of the SSME HGM was modeled using boundary fitted curvilinear coordinates and the General Interpolants Method (GIM) code. Slight compressibility of the subsonic flow was modeled using a linearized equation of state with artificial compressibility. A time relaxation method was used to obtain a steady state solution. The feasibility and potential usefulness of computational methods in assisting the design of SSME components which involves the flow of fluids within complex geometrical shapes is demonstrated.

Thoenes, J.↗

Gravity waves in Titan's atmosphere

Scintillations (high frequency variations) observed in the radio signal during the occultation of Voyager 1 by Titan (Hinson and Tyler, 1983) provide information concerning neutral atmospheric density fluctuations on scales on hundreds of meters to a few kilometers. Those seen at altitudes higher than 25 km above the surface were interpreted by Hinson and Tyler as being caused by linear, freely propagating (energy-conserving) gravity waves, but this interpretation was found to be inconsistent with the scintillation data below the 25-km altitude level. Here an attempt is made to interpret the entire scintillation profile between the surface and the 90-km altitude level in terms of gravity waves generated at the surface. Numerical calculations of the density fluctuations caused by two-dimensional, nonhydrostatic, finite-amplitude gravity waves propagating vertically through Titan's atmosphere are performed to produce synthetic scintillation profiles for comparison with the observations. The numerical model accurately treats the effects of wave transience, nonlinearity, and breakdown due to convective instability in the overturned part of the wave. The high-altitude scintillation data were accurately recovered with a freely propagating wave solution, confirming the analytic model of Hinson and Tyler. It is found that the low-altitude scintillation data can be fit by a model where a component of the gravity waves becomes convectively unstable and breaks near the 15 km level. The large-scale structure of the observed scintillation profile in the entire altitude range between 5 and 85 km can be simulated by a model where the freely propagating and breaking waves are forced at the surface simultaneously. Further analysis of the Voyager 1 Titan low-altitude scintillation data, using inversion theory appropriate for strong scattering, could potentially remove some of the ambiguities remaining in this analysis and allow a better determination of the strength and source of the waves.

Friedson, A. James↗

Physical Models for Accreting Pulsars at High Luminosity

A new window for better understanding the accretion onto strongly magnetized neutron stars in X-ray binaries is opening. In these systems the accreted material follows the magnetic field lines as it approaches the neutron star, forming accretion columns above the magnetic poles. The plasma falls toward the neutron star surface at near-relativistic speeds, losing energy by emitting X-rays. The X-ray spectral continua are commonly described using phenomenological models, i.e., power laws with different types of curved cut-offs at higher energies. Here we consider high luminosity pulsars. In these systems the mass transfer rate is high enough that the accreting plasma is thought to be decelerated in a radiation-dominated radiative shock in the accretion columns. While the theory of the emission from such shocks had already been developed by 2007, a model for direct comparison with X-ray continuum spectra in xspec or isis has only recently become available. Here we analyze the broadband X-ray spectra of the accreting pulsars Centaurus X-3 and 4U1626-67 obtained withNuSTAR. We present results from traditional empirical modeling as well as successfully apply the radiation-dominated radiative shock model. We also fit the energy-dependent pulse profiles of 4U 1626-67 using a new relativistic light bending model.

Pottschmidt, Katja↗

The ROSAT spectrum of 3C 351 - A warm absorber in an X-ray-'quiet' quasar?

3C 351 is one of the most X-ray-quiet radio quasars (alpha(ox) about 1.6). We have observed 3C 351 with the ROSAT position sensitive proportional counter (PSPC) and find a complex X-ray spectrum which is not well reproduced by a power law plus low-energy cut-off model. Soft excess, partial covering, and 'warm absorber' models can all produce acceptable fits, although only the warm absorber model gives typical values for the high-energy continuum slope. The alpha(ox) measured by using quasi-simultaneous ROSAT, MMT, and HST observations is in the range 1.5-1.6, significantly above the average of 1.37 for a complete sample of 33 3CR quasars. If the soft excess or partial covering models are correct, 3C 351 appears X-ray-quiet in the PSPC band because it has an extremely steep or flat intrinsic high-energy spectral slope. However, if the warm absorber model is correct, the quasar is intrinsically X-ray-quiet; the normalization of the intrinsic (unabsorbed) X-ray emission is unusually low relative to the optical luminosity. We investigate the properties of our warm absorber model in some detail. The apparently complicated behavior of the fit parameters may be understood by considering the effects of changing absorbing column and ionization parameter on intrinsic power-law-spectra of different slopes.

Fiore, Fabrizio↗

Model Robust Calibration: Method and Application to Electronically-Scanned Pressure Transducers

This article presents the application of a recently developed statistical regression method to the controlled instrument calibration problem. The statistical method of Model Robust Regression (MRR), developed by Mays, Birch, and Starnes, is shown to improve instrument calibration by reducing the reliance of the calibration on a predetermined parametric (e.g. polynomial, exponential, logarithmic) model. This is accomplished by allowing fits from the predetermined parametric model to be augmented by a certain portion of a fit to the residuals from the initial regression using a nonparametric (locally parametric) regression technique. The method is demonstrated for the absolute scale calibration of silicon-based pressure transducers.

Walker, Eric L.↗

Excited state dynamics of thulium ions in yttrium aluminum garnets

The processes that take place in the excited states of a trivalent Thulium (Tm) ion in an Yttrium Aluminum Garnet (YAG) crystal, being relevant to the use of this system for laser applications, have been the object of several studies. We have reexamined this system focusing our attention on the dynamics of Tm following its excitation in the H-3(sub 4) level. Under these conditions the system relaxes through a cross-relaxation process. H-3(sub 4) yields F-3(sub 4), H-3(sub 6) yields F-3(sub 4), whose rate depends upon both the concentration of the Tm ion and the temperature of the crystal. The excitation spectrum obtained by monitoring the 1.8 micron emission of Tm (due to the F-3(sub 4) yields H-3(sub 6) transition) indicates an increase in the contribution to this emission from the H-3(sub 4) level relative to the H-3(sub 5) level as the Tm concentration increases; this shows the increased role played by the H-3(sub 4) level in pumping the infrared emission. Correspondingly, the duration of the luminescence originating in the H-3(sub 4) level is shortened as the concentration of Tm increases. The concentration quenching of this lifetime can be fit to a model which assumes that the cross-relaxation is due to a dipole-dipole interaction; from this fit, the intrinsic Tm lifetime in the absence of cross relaxation can be derived. We have used this lifetime to calculate the rate of the cross-relaxation process. We have evaluated this rate as a function of the temperature and found it to be fastest at 77 K. We have also calculated the microscopic interaction parameters for the cross-relaxation process by using two independent experimental features: (1) the time evolution of the emission from the H-3(sub 4) level; and (2) the spectral overlap between the H-3(sub 4) yields F-3(sub 4) emission and the H-3(sub 6) yields F-3(sub 4) absorption. We have also considered the migration of excitation among the Tm ions in the F-3(sub 4) level and calculated the relevant microparameter by the use of the relevant spectral overlap. The data are consistent with the model in which the Tm ions, once excited into the H-3(sub 4) level decay by cross-relaxation to the F-3(sub 4), and then transfer rapidly their energy to other Tm ions.

Armagan, G.↗

Development of a New Approach to Instrument Model Design Used by Team X

The Jet Propulsion Laboratory's Advanced Design Team was formed in April 1995 to improve the quality and reduce the cost of JPL proposals and advanced mission studies. Currently a consolidation attempt is underway to develop a Model Library for use by JPL's Advanced Projects Design Team by collecting existing instrument models for inclusion in the library. This will allow users to readily find models of interest. In addition to this, there is also an attempt underway to develop a new approach to instrument model design used by the Advanced Design Team (Team X). This new approach consists of splitting up the different model parts such as orbital parameters, instrument parameters and instrument outputs into separate searchable parts. The user can then decide between design trades and use the different pieces to construct a model that will fit their needs. As well, this will lead to the opportunity for the large variety of usable instrument models.

Andrew, Shanna E.↗

Interactions of Mean Climate Change and Climate Variability on Food Security Extremes

Recognizing that climate change will affect agricultural systems both through mean changes and through shifts in climate variability and associated extreme events, we present preliminary analyses of climate impacts from a network of 1137 crop modeling sites contributed to the AgMIP Coordinated Climate-Crop Modeling Project (C3MP). At each site sensitivity tests were run according to a common protocol, which enables the fitting of crop model emulators across a range of carbon dioxide, temperature, and water (CTW) changes. C3MP can elucidate several aspects of these changes and quantify crop responses across a wide diversity of farming systems. Here we test the hypothesis that climate change and variability interact in three main ways. First, mean climate changes can affect yields across an entire time period. Second, extreme events (when they do occur) may be more sensitive to climate changes than a year with normal climate. Third, mean climate changes can alter the likelihood of climate extremes, leading to more frequent seasons with anomalies outside of the expected conditions for which management was designed. In this way, shifts in climate variability can result in an increase or reduction of mean yield, as extreme climate events tend to have lower yield than years with normal climate.C3MP maize simulations across 126 farms reveal a clear indication and quantification (as response functions) of mean climate impacts on mean yield and clearly show that mean climate changes will directly affect the variability of yield. Yield reductions from increased climate variability are not as clear as crop models tend to be less sensitive to dangers on the cool and wet extremes of climate variability, likely underestimating losses from water-logging, floods, and frosts.

Variability↗

The SkyView Data Model: Handling Heterogeneous Data Using Dynamic Software

The SkyView Virtual Telescope allows users to query over 75 distinct datasets of astronomical imagery including more than 190 distinct surveys. These holdings are extraordinarily diverse ranging over more than 15 orders of magnitude in frequency and5 orders of magnitude in resolution. To accommodate supporting these data, SkyView has developed a simple data model for surveys. Nor surprisingly, given the diverse origins of the datasets included, some surveys do not easily fit within this model. The conventional approach to dealing with such outliers is to build a layer on top of the actual data which more faithfully adheres to the data model. With SkyView we have used a somewhat alternative approach. We allow the data model to define new behaviors in the SkyView code itself either by customizing the use of existing capabilities or pointing to new classes that have non-standard behavior. This approach is feasible since in modern computer languages the code that is run can be selected dynamically with new classes completely overriding traditional behavior. This paper describes this approach and considers how it may become more useful as our software structures become increasingly dynamically linked.

Thomas McGlynn↗

Empirical Modeling of Plant Gas Fluxes in Controlled Environments

As humans extend their reach beyond the earth, bioregenerative life support systems must replace the resupply and physical/chemical systems now used. The Controlled Ecological Life Support System (CELSS) will utilize plants to recycle the carbon dioxide (CO2) and excrement produced by humans and return oxygen (O2), purified water and food. CELSS design requires knowledge of gas flux levels for net photosynthesis (PS(sub n)), dark respiration (R(sub d)) and evapotranspiration (ET). Full season gas flux data regarding these processes for wheat (Triticum aestivum), soybean (Glycine max) and rice (Oryza sativa) from published sources were used to develop empirical models. Univariate models relating crop age (days after planting) and gas flux were fit by simple regression. Models are either high order (5th to 8th) or more complex polynomials whose curves describe crop development characteristics. The models provide good estimates of gas flux maxima, but are of limited utility. To broaden the applicability, data were transformed to dimensionless or correlation formats and, again, fit by regression. Polynomials, similar to those in the initial effort, were selected as the most appropriate models. These models indicate that, within a cultivar, gas flux patterns appear remarkably similar prior to maximum flux, but exhibit considerable variation beyond this point. This suggests that more broadly applicable models of plant gas flux are feasible, but univariate models defining gas flux as a function of crop age are too simplistic. Multivariate models using CO2 and crop age were fit for PS(sub n), and R(sub d) by multiple regression. In each case, the selected model is a subset of a full third order model with all possible interactions. These models are improvements over the univariate models because they incorporate more than the single factor, crop age, as the primary variable governing gas flux. They are still limited, however, by their reliance on the other environmental conditions under which the original data were collected. Three-dimensional plots representing the response surface of each model are included. Suitability of using empirical models to generate engineering design estimates is discussed. Recommendations for the use of more complex multivariate models to increase versatility are included.

Cornett, Jessie David↗

Feasibility of quasi-random band model in evaluating atmospheric radiance

The use of the quasi-random band model in evaluating upwelling atmospheric radiation is investigated. The spectral transmittance and total band adsorptance are evaluated for selected molecular bands by using the line by line model, quasi-random band model, exponential sum fit method, and empirical correlations, and these are compared with the available experimental results. The atmospheric transmittance and upwelling radiance were calculated by using the line by line and quasi random band models and were compared with the results of an existing program called LOWTRAN. The results obtained by the exponential sum fit and empirical relations were not in good agreement with experimental results and their use cannot be justified for atmospheric studies. The line by line model was found to be the best model for atmospheric applications, but it is not practical because of high computational costs. The results of the quasi random band model compare well with the line by line and experimental results. The use of the quasi random band model is recommended for evaluation of the atmospheric radiation.

Tiwari, S. N.↗

Constraining the Metallicities, Ages, Star Formation Histories, and Ionizing Continua of Extragalactic Massive Star Populations

We infer the properties of massive star populations using the far-ultraviolet stellar continua of 61 star-forming galaxies: 42 at low redshift observed with the Hubble Space Telescope and 19 at z ~ 2 from the MegaSaura sample. We fit each stellar continuum with a linear combination of up to 50 single-age and single-metallicity starburst99 models. From these fits, we derive light-weighted ages and metallicities, which agree with stellar wind and photospheric spectral features, and infer the spectral shapes and strengths of the ionizing continua. Inferred light-weighted stellar metallicities span 0.05–1.5 Z(sub ⊙) and are similar to the measured nebular metallicities. We quantify the ionizing continua using the ratio of the ionizing flux at 900 Å to the non-ionizing flux at 1500 Å and demonstrate the evolution of this ratio with stellar age and metallicity using theoretical single-burst models. These single-burst models only match the inferred ionizing continua of half of the sample, while the other half are described by a mixture of stellar ages. Mixed-age populations produce stronger and harder ionizing spectra than continuous star formation histories, but, contrary to previous studies that assume constant star formation, have similar stellar and nebular metallicities. Stellar population age and metallicity affect the far-UV continua in different and distinguishable ways; assuming a constant star formation history diminishes the diagnostic power. Finally, we provide simple prescriptions to determine the ionizing photon production efficiency (ξ(sub ion)) from the stellar population properties. The ξ(sub ion) inferred from the observed star-forming galaxies has a range of log(ξ(sub ion)) = 24.4–25.7 Hz erg(exp −1) that depends on the stellar population age, metallicity, star formation history, and contributions from binary star evolution. These stellar population properties must be observationally determined to accurately determine the number of ionizing photons generated by massive stars.

Chisholm, J.↗

The ages of globular cluster stars - Effects of rotation on pre-main-sequence, main-sequence, and turnoff evolution

Evolutionary sequences for low-metallicity stars (Z ranging from 0.001 to 0.0001) to study the effects of internal stellar rotation on the evolutionary time scales in the pre-main sequence, the main sequence (MS), and around the MS turnoff. Although a substantial amount of angular momentum remains in the interior, rotation is only a minor perturbation on the structure and ages of globular cluster stars. Even models with large initial angular momenta have MS lifetimes that are within 1 percent of those of standard models of the same mass and composition. Therefore, rotation does not affect age estimates of globular clusters from isochrone fitting. Furthermore, the models suggest that because rotation is not likely to affect horizontal-branch (HB) morphology, it does not affect significantly age estimates from the Delta-V method. Nevertheless, the internal angular momentum in the models is consistent with observations of surface rotational velocities on the HB, which require the preservation of a large reservoir of internal angular momentum.

Deliyannis, Constantine P.↗

The Origin and Distribution of Heavy Elements in the CD Groups MKW 4 and AWM 4

At this point, data for MKW4 have been received and those for AWM4 are still awaited. The MKW4 data have been fully analyzed and a complete manuscript is available and in final review by co-authors before submission for refereed publication. The following is a summary of our principal findings: We examined the distribution and properties of the hot gas which makes up the group halo. The inner halo shows some signs of structure, with circular or elliptical beta-models providing a poor fit to the surface brightness profile. This may be evidence of large-scale motion in the inner halo, but we do not find evidence of sharp fronts or edges in the emission. The temperature of the halo declines in the core, with deprojected spectral fits showing a central temperature of approximately 1.3keV compared to approximately 3keV at 250 arcsec. However, cooling flow models provide poor fits to the inner regions of the group and the estimated cooling time of the gas is long except within the central dominant galaxy, NGC4073. Abundance profiles show a sharp increase in the core of the group. We conclude that MKW4 is a fairly relaxed group, which has developed a strong central temperature gradient but not a large-scale cooling flow.

Vrtilek, Jan↗

Monte Carlo Model Insights into the Lunar Sodium Exosphere

Sodium in the lunar exosphere is released from the lunar regolith by several mechanisms. These mechanisms include photon stimulated desorption (PSD), impact vaporization, electron stimulated desorption, and ion sputtering. Usually, PSD dominates; however, transient events can temporarily enhance other release mechanisms so that they are dominant. Examples of transient events include meteor showers and coronal mass ejections. The interaction between sodium and the regolith is important in determining the density and spatial distribution of sodium in the lunar exosphere. The temperature at which sodium sticks to the surface is one factor. In addition, the amount of thermal accommodation during the encounter between the sodium atom and the surface affects the exospheric distribution. Finally, the fraction of particles that are stuck when the surface is cold that are rereleased when the surface warms up also affects the exospheric density. In [1], we showed the "ambient" sodium exosphere from Monte Carlo modeling with a fixed source rate and fixed surface interaction parameters. We compared the enhancement when a CME passes the Moon to the ambient conditions. Here, we compare model results to data in order to determine the source rates and surface interaction parameters that provide the best fit of the model to the data.

Hurley, Dana M.↗

Chain pooling to minimize prediction error in subset regression

Monte Carlo studies using population models intended to represent response surface applications are reported. Simulated experiments were generated by adding pseudo random normally distributed errors to population values to generate observations. Model equations were fitted to the observations and the decision procedure was used to delete terms. Comparison of values predicted by the reduced models with the true population values enabled the identification of deletion strategies that are approximately optimal for minimizing prediction errors.

Holms, A. G.↗

Are Compton-thin AGNs Globally Compton Thin?

We select eight nearby active galactic nuclei (AGNs) which, based on previous work, appear to be Compton-thin in the line of sight. We model with MYTORUS their broadband X-ray spectra from 20 individual observations with Suzaku, accounting self-consistently for Fe Kα line emission, as well as direct and scattered continuum from matter with finite column density and solar Fe abundance. Our model configuration allows us to measure the global, out of the line of sight, equivalent hydrogen column density separately from that in the line of sight. For 5 out of 20 observations (in 3 AGNs) we find that the global column density is in fact ≳1.5 × 1024 per sq.cm, consistent with the distant scattering matter being Compton-thick. For a fourth AGN, two out of five observations are also consistent with being Compton-thick, although with large errors. Some of these AGNs have been reported to host relativistically broadened Fe Kα emission. Based on our modeling, the Fe Kα emission line is not resolved in all but two Suzaku observations, and the data can be fitted well with models that only include a narrow Fe Kα emission line.

P. Tzanavaris↗

ADP study of gamma-ray bursts

This grant supported study of cyclotron scattering lines in the spectra of gamma-ray bursts through analysis of Ginga and HEAO-1 archival data, and modeling of the results in terms of radiation transfer calculations of cyclotron scattering in a strong magnetic field. A Monte Carlo radiation transfer code with which we are able to calculate the expected properties of cyclotron scattering lines in the spectra of gamma-ray bursts was developed. The extensive software necessary in order to carry out fits of these model spectra to gamma-ray burst spectral data, including folding of the model spectra through the detector response functions was also developed. Fits to Ginga satellite data on burst GB880205 were completed and fits to Ginga satellite data on burst GB870303 are being carried out. These fits have allowed us to test our software, as well as to garner new scientific results. This work has demonstrated that cyclotron resonant scattering successfully accounts for the locations, strengths, and widths of the observed line features in GB870303 and GB880205. The success of the model provides compelling evidence that these gamma-ray bursts come from strongly magnetic neutron stars and are galactic in origin, resolving longstanding controversies about the nature and distance of the burst sources. These results were reported in two papers which are in press in the proceedings of the Taos Workshop on Gamma-Ray Bursts, and in a paper submitted for publication.

Lamb, Don Q.↗