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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

Low-temperature transonic cooling flows in galaxy clusters

Calculations are presented which demonstrate that cooling flow models with large sonic radii may be consistent with observed cluster gas properties. It is found that plausible cluster parameters and cooling flow mass accretion rates can produce sonic radii of 10-20 kpc for sonic point temperatures of 1-3 x 10 to the 6th K. The numerical calculations match these cooling flows to hydrostatic atmosphere solutions for the cluster gas beyond the cooling flow region. The cooling flows produce no appreciable 'holes' in the surface brightness toward the cluster center, and the model can be made to match the observed X-ray surface brightness of three clusters in which cooling flows had been believed to be absent. It is suggested that clusters with low velocity dispersion may be the natural location for such 'cool' cooling flows, and fits of these models to the X-ray surface brightness profiles for three clusters are presented.

Sulkanen, Martin E.↗
NASA NTRS2005

The Effects of Atmospheric Opacity on the Seasonal Variation of Martian Surface Temperature

The daily and seasonal variation of surface temperature is a central element in the description of martian climate. Surface thermal inertia and albedo are critical boundary inputs for simulating surface temperature in Mars general circulation models (MGCMs). Thermal inertia (TI) is also of intrinsic interest as it may be related to regolith properties such as particle size and surface character and so high spatial resolution is desirable. The recent mapping of TI at very high (0.25 deg) spatial resolution was achieved by fitting a thermal model to surface temperature observations obtained over a broad range of several martian years. However, varying atmospheric opacity (dust and water ice clouds) can significantly influence the estimated TI field and this effect was not fully compensated for. Opacity leads to an increase in morning temperature and a decrease in afternoon temperature, thus increasing the apparent thermal inertia.

Wilson, R. J.↗
NASA NTRS1985

Seyfert galaxy narrow-line regions. II - Kinematic models

A class of kinematic models for the narrow-line regions of active galactic nuclei is presented. The models incorporate continuity-constrained radial flow of emitting material in a bounded, spherical geometry. They are designed to represent the narrow-line regions in as simple a manner as is consistent with current knowledge. The systematic comparison of line profiles arising from the models to the spectra in the forbidden O III lambda 5007 emission line which have been described in Paper I of this series is used as a probe of the line-emitting region kinematics: the model parameter areas of interest are identified and tested for consistency with the data. The results indicate heterogeneous kinematic conditions; the most common pattern is outflow of emitting material in a radially decelerating velocity field in the presence of distributed obscuration. However, some observed line profiles can be fitted by several model profiles from very different kinematic conditions, indicating that the information in a spatially unresolved line profile is insufficient to constrain uniquely the geometry and kinematics of the line-emitting region. Although the specific conclusions are model dependent, a predominantly radial pattern of velocities and a distribution of obscuring material are probably essential features of any successful kinematic model of the narrow-line regions.

Vrtilek, J. M.↗
NASA NTRS1994

Simulations of deep galaxy fields. 1: Monte Carlo simulations of optical and near-infrared counts

Monte Carlo simulations of three-dimensional galaxy distributions are performed, following the 1988 prescription of Chokshi & Wright, to study the photometric properties of evolving galaxy populations in the optical and near-infrared bands to high redshifts. In this paper, the first of a series, we present our baseline model in which galaxy numbers are conserved, and in which no explicit 'starburst' population is included. We use the model in an attempt to simultaneously fit published blue and near-infrared photometric and spectroscopic observations of deep fields. We find that our baseline models, with a formation redshift, z(sub f), of 1000, and H(sub 0) = 50, are able to reproduce the blue counts to b(sub j) = 22, independent of the value of Omega(sub 0), and also to provide a satisfactory fit to the observed blue-band redshift distributions, but for no value of Omega(sub 0) do we achieve an acceptable fit to the fainter blue counts. In the K band, we fit the number counts to the limit of the present-day surveys only for an Omega(sub 0) = 0 cosmology. We investigate the effect on the model fits of varying the cosmological parameters H(sub 0), the formation red-shift z(sub f), and the local luminosity function. Changing H(sub 0) does not improve the fits to the observations. However, reducing the epoch of a galaxy formation used in our simulations has a substantial effect. In particular, a model with z(sub f) approximately equal to 5 in a low Omega(sub 0) universe improves the fit to the faintest photometric blue data without any need to invoke a new population of galaxies, substantial merging, or a significant starburst galaxy population. For an Omega(sub 0) = 1 universe, however, reducing z(sub f) is less successful at fitting the blue-band counts and has little effect at all at K. Varying the parameters of the local luminosity function can also have a significant effect. In particular the steep low end slope of the local luminosity function of Franceschini et al. allows an acceptable fit to the b(sub j) less than or equal to 25 counts for Omega(sub 0) = 1, but is incompatible with Omega(sub 0) = 0.

Chokshi, Arati↗
NASA NTRS1992

Global O VI line emission from the Cygnus Loop supernova remnant and direct kinematic measurement of the associated shock

A far-ultraviolet spectrophotometric emission-line mapping of the Cygnus Loop supernova remnant is presented. These are results from the first flight of the rocket-borne, High Resolution Emission Line Spectrometer. The spatial distribution of the emission is that of a limb-brightened shell, and similar to soft X-ray maps. The emission-line profiles, which are broader than the instrument resolution, were consistent with uniformly expanding shell models. Best-fit values give a radial expansion velocity to the emissive region of 185(+/-19) km/s and a reddening-corrected average surface brightness of 8.8(+/-3.6) x 10 exp -6 ergs/sq cm s sr in the doublet. Comparison of the observed brightness with predictions of both radiative and nonradiative shock models provides constraints for the global blast wave ram pressure as well as a "covering factor" of the intermediate velocity shock.

Rasmussen, Andrew↗
NASA NTRS2013

Erratum: "Photoionization Modeling of Oxygen K Absorption in the Interstellar Medium, the Chandra Grating Spectra of XTE J1817-330" (2013, Apj, 768, 60)

In the published version of this paper, there are some minor inaccuracies in the absorption-line wavelengths listed in Table 4 as a result of a faulty reduction procedure of the Obs6615 spectrum. The shifts have been detected in a comparison with the wavelengths listed for this spectrum in the Chandra Transmission Grating Catalog and Archive (TGCat8). They are due to incorrect centroid positions of the zero-order image in both reductions as determined by the tgdetect utility which, when disentangled, yield the improved line positions of the amended Table 4 given below. It must also be pointed out that other quantitative findings of the original paper: 1. Table 5, p. 9: the column density (NH), ionization parameter, oxygen abundance of the warmabs model and the normalization and photon index of the power-law model; 2. Table 6, p. 9: the hydrogen column density of the warmabs fit; 3. Table 7, p. 9: the present oxygen equivalent widths of XTE J1817-330; and 4. Table 8, p. 10: the present oxygen column densities of XTE J1817-330 derived from both the curve of growth and warmabs model fit have been revised in the new light and are, within the estimated uncertainty ranges, in good accord with the new rendering.

abundance of the warmabs↗
NASA NTRS1997

Molecular Hydrogen in Diffuse Interstellar Clouds of Arbitrary Three-Dimensional Geometry

We have constructed three-dimensional models for the equilibrium abundance of molecular hydrogen in diffuse interstellar clouds of arbitrary geometry that are illuminated by ultraviolet radiation. The position-dependent photodissociation rate of H2 in such clouds was computed with a 26 ray approximation to model the attenuation of the incident ultraviolet radiation field by dust and by H2 line absorption. We have applied our modeling technique to the isolated diffuse cloud G236+39, assuming that the cloud has a constant density and that the thickness of the cloud along the line of sight is at every point proportional to the 100 micron continuum intensity measured by IRAS. We find that our model can successfully account for observed variations in the ratio of 100 micron continuum intensity to H I column density, with larger values of that ratio occurring along lines of sight in which the molecular hydrogen fraction is expected to be the largest. Using a standard chi square analysis to assess the goodness of fit of our models, we find (at the 60 a level) that a three-dimensional model is more successful in matching the observational data than a one-dimensional model in which the geometrical extent of the cloud along the line of sight is assumed to be much smaller than its extent in the plane of the sky. If D is the distance to G236 + 39, and given standard assumptions about the rate of grain-catalyzed H2 formation, we find that the cloud has an extent along the line of sight that is 0.9 +/- 0.1 times its mean extent projected onto the plane of the sky and a gas density of 53 +/- 8 (100 pc/D) H nuclei/cc and is illuminated by a radiation field of 1.1 +/- 0.2 (100 pc/D) times the mean interstellar radiation field. The derived 100 micron emissivity per nucleon is 1.13 +/- 0.06 x 10(exp -20) MJy/sr sq cm.

Spaans, Marco↗
NASA NTRS2014

Selection of Thermal Worst-Case Orbits via Modified Efficient Global Optimization

Efficient Global Optimization (EGO) was used to select orbits with worst-case hot and cold thermal environments for the Stratospheric Aerosol and Gas Experiment (SAGE) III. The SAGE III system thermal model changed substantially since the previous selection of worst-case orbits (which did not use the EGO method), so the selections were revised to ensure the worst cases are being captured. The EGO method consists of first conducting an initial set of parametric runs, generated with a space-filling Design of Experiments (DoE) method, then fitting a surrogate model to the data and searching for points of maximum Expected Improvement (EI) to conduct additional runs. The general EGO method was modified by using a multi-start optimizer to identify multiple new test points at each iteration. This modification facilitates parallel computing and decreases the burden of user interaction when the optimizer code is not integrated with the model. Thermal worst-case orbits for SAGE III were successfully identified and shown by direct comparison to be more severe than those identified in the previous selection. The EGO method is a useful tool for this application and can result in computational savings if the initial Design of Experiments (DoE) is selected appropriately.

Moeller, Timothy M.↗
NASA NTRS1987

Comparing the jerk with other global models of the geomagnetic field from 1960 to 1978

About 3300 satellite values of geomagnetic intensity and about 700 observatory values of annual mean magnetic vector components from 1960 to 1978 were fitted by three global models of the geomagnetic field B. Each model includes a spatially constant external field whose time dependence is a constant plus another constant times the Dst index, and each model accepts a time-independent station correction at each observatory. The time dependence of the internal Gauss coefficients is either cubic, quintic, or biquadratic (two independent quadratics, one before and one after January 1, 1970); and g1(0) also has an induced term proportional to the Dst index. The rms residual of the data fit is the same for the cubic and biquadratic models and insignificantly smaller for the quintic model. The quintic and biquadratic models have 1164 adjustable parameters, and the cubic has 1038. At a high level of significance the parameters of the best fitting biquadratic rule out a physical model for the magnetic impulse of 1969 in which the level surfaces of electrical conductivity in the lower mantle are approximately spherical, and the radial magnetic field at the core-mantle boundary goes from one quadratic time dependence to another in a year or less.

Backus, G. E.↗
NASA NTRS1990

Accretion tori and cones of ionizing radiation in Seyfert galaxies

The photoionization of extended narrow-line regions in Seyfert galaxies by the radiation produced in a thick accretion disk is studied. The emission-line spectrum is calculated for a range of black hole masses, varying the values of the ionization parameter and the disk size. It is found that models with a million solar masses fit observations of very large accretion disk sizes, while models with 10 million solar masses fit them better with smaller disks. The latter models are preferable since they have lower super-Eddington accretion rates.

Acosta-Pulido, Jose A.↗
NASA NTRS2001

Production, Outflow Velocity, and Radial Distribution of H2O and OH in the Coma of Comet C/1995 O1 [Hale-Bopp] from Wide Field Imaging of OH

Observations of OH are a useful proxy of the water production rate (Q(sub H2O)) and outflow velocity (V(sub out)) in comets. We use wide field images taken on 03/28/1997 and 04/08/1997 that capture the entire scale length of the OH coma of comet C/1995O1 (Hale-Bopp) to obtain Q(sub H2O) from the model-independent method of aperture summation. We also extract the radial brightness profile of OH 3080 angstroms out to cometocentric distances of up to 10(exp 6) km using an adaptive ring summation algorithm. Radial profiles are obtained as azimuthal averages and in quadrants covering different position angles relative to the comet-Sun line. These profiles are fit using both fixed and variable velocity two-component spherical expansion models to determine VOH with increasing distance from the nucleus. The OH coma of Hale-Bopp was more spatially extended than in previous comets, and this extension is best matched by a variable acceleration of H2O and OH that acted across the entire coma, but was strongest within 1-2 x 10(exp 4) km from the nucleus. This acceleration led to VOH at 10(exp 6) km that was 2-3 times greater than that obtained from a 1P/Halleytype comet at 1 AU, a result that is consistent with gas-kinetic models, extrapolation from previous observations of OH in comets with Q(sub H2O) > 10(exp 29)/s, and radio measurements of the outer coma Hale-Bopp OH velocity profile. When the coma is broken down by quadrant, we find an azimuthal asymmetry in the radial distribution that is characterized by an increase in the spatial extent of OH in the region between the orbit-trailing and anti-sunward directions. Model fits to this area and comparison with radio OH measurements suggest greater acceleration in this region, with VOH UP to 1.5 times greater at 10(exp 6) km radial distance than elsewhere in the coma.

Harris, Walter M.↗
NASA NTRS2021

Potential origin of the state-dependent high-energy tail in the black hole microquasar Cygnus X-1 as seen with INTEGRAL

Context. 0.1–10 MeV observations of the black hole microquasar Cygnus X-1 have shown the presence of a spectral feature in the form of a power law in addition to the standard black body (0.1–10 keV) and Comptonization (10–200 keV) components observed by INTEGRAL in several black-hole X-ray binaries. This so-called “high-energy tail” was recently shown to be strong in the hard spectral state of Cygnus X-1, and, in this system, has been interpreted as the high-energy part of the emission from a compact jet. Aims. This result was nevertheless obtained from a data set largely dominated by hard state observations. In the soft state, only upper limits on the presence and hence the potential parameters of a high-energy tail could be derived. Using an extended data set, we aim to obtain better constraints on the properties of this spectral component in both states. Methods. We make use of data obtained from about 15 years of observations with the INTEGRAL satellite. The data set is separated into the different states and we analyze stacked state-resolved spectra obtained from the X-ray monitors, the gamma-ray imager, and the gamma-ray spectrometer (SPI) onboard. Results. A high-energy component is detected in both states, confirming its earlier detection in the hard state and its suspected presence in the soft state with INTEGRAL, as seen in a much smaller SPI data set. We first characterize the high-energy tail components in the two states through a model-independent, phenomenological analysis. We then apply physical models based on hybrid Comptonization (eqpair and belm). The spectra are well modeled in all cases, with a similar goodness of the fits. While in the semi-phenomenological approach the high-energy tail has similar indices in both states, the fits with the physical models seem to indicate slightly different properties. Based on this approach, we discuss the potential origins of the high-energy components in both the soft and hard states, and favor an interpretation where the high-energy component is due to a compact jet in the hard state and hybrid Comptonization in either a magnetized or nonmagnetized corona in the soft state.

accretion↗
NASA NTRS2008

How Deep and Hot was Earth's Magma Ocean? Combined Experimental Datasets for the Metal-silicate Partitioning of 11 Siderophile Elements - Ni, Co, Mo, W, P, Mn, V, Cr, Ga, Cu and Pd

Since approximately 1990 high pressure and temperature (PT) experiments on metal-silicate systems have showed that partition coefficients (D) for siderophile (iron-loving) elements are much different than those measured at low PT conditions. The high PT data have been used to argue for a magma ocean during growth of the early Earth. Initial conclusions were based on experiments and calculations for a small number of elements such as Ni and Co. However, for many elements only a limited number of experimental data were available then, and they only hinted at values of metal-silicate D's at high PT conditions. In the ensuing decades there have been hundreds of new experiments carried out and published on a wide range of siderophile elements. At the same time several different models have been advanced to explain the siderophile elements in the earth's mantle: a) intermediate depth magma ocean; 25-30 GPa, b) deep magma ocean; up to 50 GPa, and c) early reduced and later oxidized magma ocean. Some studies have drawn conclusions based on a small subset of siderophile elements, or a set of elements that provides little leverage on the big picture (like slightly siderophile elements), and no single study has attempted to quantitatively explain more than 5 elements at a time. The purpose of this abstract is to update the predictive expressions outlined by Righter et al. (1997) with new experimental data from the last decade, test the predictive ability of these expressions against independent datasets (there are more data now to do this properly), and to apply the resulting expressions to the siderophile element patterns in Earth's upper mantle. The predictive expressions have the form: lnD = alnfO2 + b/T + cP/T + d(1Xs) + e(1Xc) + SigmafiXi + g These expressions are guided by the thermodynamics of simple metal-oxide equilibria that control each element, include terms that mimic the activity coefficients of each element in the metal and silicate, and quantify the effect of variable oxygen fugacity. Preliminary results confirm that D(Ni) and D(Co) converge at pressures near 25-30 GPa and approximately 2200 K, and show that D(Pd) and D(Cu) become too low at the PT conditions of the deepest models. Furthermore, models which force fit V and Cr mantle concentrations by metal-silicate equilibrium overlook the fact that at early Earth mantle fO2, these elements will be more compatible in Mg-perovskite and (Fe,Mg)O than in metal. Thus an intermediate depth magma ocean, at 25-30 GPa, 2200 K, and at IW-2, can explain more mantle siderophile element concentrations than other models.

Righter, Kevin↗
NASA NTRS1992

Independent Constituents in Observatory Time Series in Characterizing the Secular Variation

By an ansatz of a finite set of N superimposed oscillations, observatory time series are fitted. This start model is optimized for the independence of the used 2N-dimensional functional system as a criterion. Here, the contributions of each of the dimensions to the volume of the corresponding 2N-dimensional parallelpiped of the functional space is to be evaluated. Independence and an equal standard for the constituents are achieved by the optimizing procedure varying the start values of frequencies in the fitting ansatz when quotients of successively constructed Gram-determinants are beyond a definite threshold and of nearly the same numerical value. Results for globally distributed observatories characterize the secular variation.

Webers, W.↗
NASA NTRS1995

Bringing the cross-correlation method up to date

The cross-correlation (XC) method of Tonry & Davis (1979, AJ, 84, 1511) is generalized to arbitrary parametrized line profiles. In the new algorithm the correlation function itself, rather than the observed galaxy spectrum, is fitted by the model line profile: this removes much of the complication in the error analysis caused by template mismatch. Like the Fourier correlation quotient (FCQ) method of Bender (1990, A&A, 229, 441), the inferred line profiles are, up to a normalization constant, independent of template mismatch as long as there are no blended lines. The standard reduced chi(exp 2) is a good measure of the fit of the inferred velocity distribution, largely decoupled from the fit of the spectral template. The updated XC method performs as well as other recently developed methods, with the added virtue of conceptual simplicity.

Statler, Thomas↗
NASA NTRS2004

Current State of Modeling the Photochemistry of Titan's Mutually Dependent Atmosphere and Ionosphere

In the context of recent observations, microphysical models, and laboratory data, a photochemical model of Titan's atmosphere, including updated chemistry focusing on rate coefficients and cross sections measured under appropriate conditions, has been developed to increase understanding of these processes and improve upon previous Titan photochemical models. The model employs a two-stream discrete ordinates method to characterize the transfer of solar radiation, and the effects of electron-impact, cosmic-ray deposition, and aerosol opacities from fractal and Mie particles are analyzed. Sensitivity studies demonstrate that an eddy diffusion profile with a homopause level of 850 km and a methane stratospheric mole fraction of 2.2% provides the best fit of stratospheric and upper atmosphere observations and an improved fit over previous Titan photochemical models. Lack of fits for C3H8, HC3N, and possibly C2H3CN can be resolved with adjustments in aerosol opacity. The model presents a benzene profile consistent with its detection in Titan's stratosphere [Coustenis et al., 2003], which may play an important role in the formation of Titan hazes. An electron peak concentration of 4200 cm(exp -3) is calculated, which exceeds observations by 20%, considerably lower than previous ionosphere models. With adjustments in aerosol opacities and surface fluxes the model illustrates that reasonable fits to existing observations are possible with a single eddy diffusion profile, contrary to the conclusions of previous Titan models. These results will aid in the receipt and interpretation of data from Cassini-Huygens, which will arrive at Titan in 2004 and deploy a probe into Titan's atmosphere in January 2005.

ion chemistry↗
NASA NTRS1975

On the origin of ultraviolet constrasts on Venus

Models for the origin of the contrasts in the ultraviolet images of Venus are examined in an attempt to determine the physical differences between light and dark regions fundamental to a clear understanding of the apparent cloud motions. To evaluate the meaning of the wavelength dependence of the contrasts, an improved determination of the spherical albedo curve for Venus in the 0.225-1.06 micron range is made by fitting appropriate theoretical models to the observations of monochromatic magnitudes as a function of phase angle. It is shown that, because of differences between the spectral dependences of spherical albedo and contrasts, at least one major absorber other than the one causing the contrasts is almost certainly required. Several classes of models are compatible with present observational evidence. The contrasts and the absorption associated with them may in fact be occurring below, within or above the main visible cloud layer, and thus an unambiguous interpretation of the apparent cloud motions is not possible.

Travis, L. D.↗
NASA NTRSDate not supplied

Machine Learning Emulators and Empirical Models Combining Climate and Global Crop Models for Seasonal Agricultural Production

We present results from several connected efforts to apply machine learning methods to estimates of seasonal agricultural production anomalies around the world. First, we apply the XGBoost Random Forest method to fit emulators that mimic global crop models participating in the Agricultural Model Intercomparison and Improvement Project (AgMIP) Global Gridded Crop Model Intercomparison (GGCMI). These are the same models used in the agricultural sector simulations of the Inter-Sectoral Impacts Model Intercomparison Project (ISIMIP). These emulators use 8 climate variables split across 5 sub-seasonal representations of the growing season for each ½ degree grid cell around the world for maize, wheat, rice and soybeans. Emulators are useful for estimating conditions that have not already been simulated by GGCMI (e.g., in a seasonal prediction model) and also to diagnose model differences and capabilities. For example, emulators of the pDSSAT maize model tend to be more reliant on mean temperatures than the LPJmL model, and few models have strong responses to cold extremes. Second, we use a similar XGBoost approach to fit empirical models for national production data for the top 20 producing countries according to the United Nations Food and Agricultural Organization (FAO). Models utilize both climate observations and the GGCM models as predictors, resulting in skillful models for many (but not all) top producing-countries. The patterns of climate and crop model features selected indicate regions and systems that are better or worse simulated by the GGCMs. For example, information in cold extreme predictors is often combined with GGCM output predictors to provide sensitivity that models may underrepresent.

machine learning↗