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 901 records · Page 50
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Blind Modeling Validation Exercises Using the Horizontal Dry Cask Simulator

The U.S. Department of Energy (DOE) established a need to understand the thermal-hydraulic properties of dry storage systems for commercial spent nuclear fuel (SNF) in response to a shift towards the storage of high-burnup (HBU) fuel (> 45 gigawatt days per metric ton of uranium, or GWd/MTU). This shift raises concerns regarding cladding integrity, which faces increased risk at the higher temperatures within spent fuel assemblies present within HBU fuel compared to low-burnup fuel (≤ 45 GWd/MTU). A dry cask simulator (DCS) was built at Sandia National Laboratories (SNL) in Albuquerque, New Mexico to produce validation-quality data that can be used to test the accuracy of the modeling used to predict cladding temperatures. These temperatures are critical to evaluating cladding integrity throughout the storage cycle of commercial spent nuclear fuel. A model validation exercise was previously carried out for the DCS in a vertical configuration. Lessons learned during the previous validation exercise have been applied to a new, blind study using a horizontal dry cask simulator (HDCS). Three modeling institutions – the Nuclear Regulatory Commission (NRC), Pacific Northwest National Laboratory (PNNL), and Empresa Nacional del Uranio, S.A., S.M.E. (ENUSA) – were granted access to the input parameters from the DCS Handbook, SAND2017-13058R, and results from a limited data set from the horizontal BWR dry cask simulator tests reported in the HDCS update report, SAND2019-11688R. With this information, each institution was tasked to calculate peak cladding temperatures and air mass flow rates for ten HDCS test cases. Axial as well as vertical and horizontal transverse temperature profiles were also calculated. These calculations were done using modeling codes (ANSYS/Fluent, STAR-CCM+, or COBRA-SFS), each with their own unique combination of modeling assumptions and boundary conditions. For this validation study, the ten test cases of the horizontal dry cask simulator were defined by three independent variables – fuel assembly decay heat (0.5 kW, 1 kW, 2.5 W, and 5 kW), internal backfill pressure (100 kPa and 800 kPa), and backfill gas (helium and air). The plots provided in Chapter 3 of this report show the axial, vertical, and horizontal temperature profiles obtained from the dry cask simulator experiments in the horizontal configuration and the corresponding models used to describe the thermal-hydraulic behavior of this system. The tables provided in Chapter 3 illustrate the closeness of fit of the model data to the experiment data through root mean square (RMS) calculations of the error in peak cladding temperatures (PCTs), PCT axial locations, axial temperature profiles, vertical and horizontal temperature profiles at two different axial locations, and air mass flow rates for the ten test cases, normalized by the experimental results. The model results are assigned arbitrary model numbers to retain anonymity. Due to the relatively flat axial temperature profiles, small temperature gradients resulted in large deviations of all models’ PCT axial location from the experimental PCT axial location. When the PCT axial location error is excluded in the calculation of the combined RMS of the normalized errors that considers PCT, the temperature profiles, and the air mass flow rates, the model data fits the experimental data to within 5%. When the vault information is excluded, the model data fits the experimental data to within 2.5%. An error analysis was developed further for one model, using the model and experimental uncertainties in each validation parameter to calculate validation uncertainties. The uncertainties for each parameter were used to define quantifiable validation criteria. For this analysis, the model was considered validated for a given comparison metric if the normalized error in that metric divided by the validation uncertainty was less than or equal to 1. When considering the combined RMS of the normalized errors of all metrics divided by their validation uncertainties, the model was found to have satisfied the criterion for model validation.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2022

The dark matter halo masses of elliptical galaxies as a function of observationally robust quantities

Context. The assembly history of the stellar component of a massive elliptical galaxy is closely related to that of its dark matter halo. Measuring how the properties of galaxies correlate with their halo mass can therefore help to understand their evolution. Aims. We investigate how the dark matter halo mass of elliptical galaxies varies as a function of their properties, using weak gravitational lensing observations. To minimise the chances of biases, we focus on the following galaxy properties that can be determined robustly: the surface brightness profile and the colour. Methods. We selected 2409 central massive elliptical galaxies (log M*/M ⊙ ≳ 11.4) from the Sloan Digital Sky Survey spectroscopic sample. We first measured their surface brightness profile and colours by fitting Sérsic models to photometric data from the Kilo-Degree Survey (KiDS). We fitted their halo mass distribution as a function of redshift, rest-frame r-band luminosity, half-light radius, and rest-frame u - g colour, using KiDS weak lensing measurements and a Bayesian hierarchical approach. For the sake of robustness with respect to assumptions on the large-radii behaviour of the surface brightness, we repeated the analysis replacing the total luminosity and half-light radius with the luminosity within a 10 kpc aperture, L r, 10 , and the light-weighted surface brightness slope, Γ 10 . Results. We did not detect any correlation between the halo mass and either the half-light radius or colour at fixed redshift and luminosity. Using the robust surface brightness parameterisation, we found that the halo mass correlates weakly with L r,10 and anti-correlates with Γ 10 . At fixed redshift, L r, 10 and Γ 10 , the difference in the average halo mass between galaxies at the 84th percentile and 16th percentile of the colour distribution is 0.00 ± 0.11 dex. Conclusion. Our results indicate that the average star formation efficiency of massive elliptical galaxies has little dependence on their final size or colour. This suggests that the origin of the diversity in the size and colour distribution of these objects lies with properties other than the halo mass.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The GMRT archive atomic gas survey – II. Mass modelling and dark matter halo properties across late-type spirals

Studying the kinematics and mass modelling of galaxies from H i 21 cm data provides valuable insights into the properties of both the baryonic components and the dark matter halo in nearby galaxies. Despite many observational studies, mass modelling of galaxies remains challenging due to different limitations. For example, most of the previous studies involving mass modelling are based on rotation curves derived from 2D velocity fields from H i or H α spectroscopic observation which are often affected by beam smearing and projection effect. However, kinematic modelling done by fitting the ‘Tilted ring model’ to 3D data cube is not affected by these issues. In this study, we present and compare 3D kinematic modelling of a pilot sample of 11 galaxies from the GMRT archive atomic gas survey (GARCIA) using two different publicly available pipelines. We model the observed H i rotation curve using 3.6-μm infrared data and SDSS r-band data for stellar contribution, H i surface density profile for gas, and Navarro–Frenk–White profile for dark matter halo; and employ the Markov chain Monte Carlo optimization method for parameter estimation. Further, to validate our analysis, we revisit important scaling relations, e.g. the M gas –M star relation, M star –M halo relation, M gas –M halo relation and Baryonic Tully–Fisher relation. Here, the scaling relations from our analysis are broadly consistent with that reported in the literature. A larger sample of galaxies from GARCIA in the near future will allow studying these scaling relations in greater details.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Connecting relativistic density functional theory to microscopic calculations

The development of systematic effective field theories (EFTs) for nuclear forces and advances in solving the nuclear many-body problem have greatly improved our understanding of dense nuclear matter and the structure of finite nuclei. For global nuclear calculations, density functional theories (DFTs) have been developed to reduce the complexity and computational cost required in describing nuclear systems. However, DFT often makes approximations and assumptions about terms included in the functional, which may introduce systematic uncertainties compared to microscopic calculations using EFTs. In this work, we investigate possible avenues of improving nuclear DFT using nonlinear relativistic mean-field (RMF) theory. We explore the impact of RMF model extensions by fitting the nonlinear RMF model to predictions of nuclear matter and selected closed-shell nuclei using four successful chiral EFT Hamiltonians. We find that these model extensions are impactful and important in capturing the physics present within chiral Hamiltonians, particularly for charge radii and neutron skins of closed-shell nuclei. However, there are additional effects that are not captured within the RMF model, particularly within the isoscalar sector of RMF theory. Additional model extensions and the reliability of the nonlinear RMF model are discussed.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2020

The weakly nonlinear response and nonaffine interpretation of the Johnson–Segalman/Gordon–Schowalter model

We derive new analytical solutions for the non-affine Johnson-Segalman/Gordon-Schowalter (JS/GS) constitutive equation with a general relaxation kernel in medium-amplitude oscillatory shear (MAOS) deformation. The results show time-strain separable (TSS) nonlinearity, therefore providing new physically-meaningful interpretation to the heuristic TSS nonlinear parameter in MAOS. The upper-convected, lower-convected, and corotational Maxwell models are all subsets of the results presented here. The model assumes that the microscale elements causing stress in the material slip compared to the continuum deformation. We introduce a visualization of the non-affine deformation field that acts on stress-generating elements to reinforce the hysical interpretation of the JS/GS class of models. Lastly, a case study is presented where previously published results, from fitting TSS models to MAOS data, can be re-interpreted based on the concept of non-affine motion of the JS/GS framework.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Stellar Mass and Stellar Mass-to-light Ratio–Color Relations for Low Surface Brightness Galaxies

We estimate the stellar mass for a sample of low surface brightness galaxies (LSBGs) by fitting their multiband spectral energy distributions (SEDs) to the stellar population synthesis model. The derived stellar masses (log M {sub *}/M {sub ⊙}) span from 7.1 to 11.1, with a mean of log M {sub *}/M {sub ⊙} = 8.5, which is lower than that for normal galaxies. The stellar mass-to-light ratio (γ*) in each band varies little with the absolute magnitude but increases with higher M {sub *}. This trend of γ* with M {sub *} is even stronger in bluer bands. In addition, the γ* for our LSBGs slightly declines from the r band to the longer-wavelength bands. The log γ{sub ∗}{sup j} (j = g, r, i, and z) have relatively tight relations with optical colors of g − r and g − i. Compared with several representative γ*–color relations (MLCRs) from the literature, our MLCRs based on LSBG data are consistently among those literature MLCRs previously defined on diverse galaxy samples, and the existing minor differences between the MLCRs are caused by the differences in the SED model ingredients (including initial mass function, star formation history, and stellar population model), line fitting techniques, galaxy samples, and photometric zero-point, rather than the galaxy surface brightness itself, which distinguishes LSBGs from high surface brightness galaxies. Our LSBGs would be very likely to follow those representative MLCRs previously defined in diverse galaxy populations, if those main ingredients were taken into account.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1979

The magnetic field of Mercury

Data from Mariner 10 observations of Mercury indicate that there exists an intrinsic magnetic field of the planet, sufficiently strong at present to deflect the solar wind flow around the planet and to form a detached bow shock wave in the super Alfvenic solar wind. Four methods used to analyze the magnetic field data and derive quantitative values for the description of the planetary field include (1) comparison of bow shock and magnetopause relative positions at Mercury to those at Earth; (2) direct spherical harmonic analysis of the data; (3) modeling of the magnetosphere by an image dipole and infinite 2-D current sheet in addition to the planetary field; and (4) scaling of a mathematical model for the terrestrial magnetosphere. The results obtained yield dipole moments ranging from 2.4 to 5.1x10 to the 22d power, with the lower values associated cw cm with certain models using partial quadrupole and octupole terms to improve the least squares fitting of models to observations.

Ness, N. F.↗
NASA NTRS1988

The debate over the Cretaceous-Tertiary boundary

Large-body impact on the Earth is a rare but indisputable geologic process. The impact rate is approximately known from objects discovered in Earth-crossing orbits and from the statistics of craters on the Earth's surface. Tektite and microtektite strewn fields constitute unmistakable ejecta deposits that can be due only to large-body impacts. The Cretaceous-Tertiary (K-T) boundary coincides with an unusually severe biological trauma, and this stratigraphic horizon is marked on a worldwide basis by anomalous concentrations of noble metals in chondritic proportions, mineral spherules with relict quench-crystallization textures, and mineral and rock grains showing shock deformation. These features are precisely compatible with an impact origin. Although only impact explains all the types of K-T boundary evidence, the story may not be as simple as once thought. The original hypothesis envisioned one large impact, triggering one great extinction. Newer evidence hints at various complications. Different challenges are faced by the occupants of each apex of a three-cornered argument over the K-T event. Proponents of a non-impact explanation must show that the evidence fits their preferred model better than it fits the impact scenario. Proponents of the single impact-single extinction view must explain away the complications. Proponents of a more complex impact crisis must develop a reasonable scenario which honors the new evidence.

Alvarez, W.↗
NASA NTRS1993

Video transmission on ATM networks

The broadband integrated services digital network (B-ISDN) is expected to provide high-speed and flexible multimedia applications. Multimedia includes data, graphics, image, voice, and video. Asynchronous transfer mode (ATM) is the adopted transport techniques for B-ISDN and has the potential for providing a more efficient and integrated environment for multimedia. It is believed that most broadband applications will make heavy use of visual information. The prospect of wide spread use of image and video communication has led to interest in coding algorithms for reducing bandwidth requirements and improving image quality. The major results of a study on the bridging of network transmission performance and video coding are: Using two representative video sequences, several video source models are developed. The fitness of these models are validated through the use of statistical tests and network queuing performance. A dual leaky bucket algorithm is proposed as an effective network policing function. The concept of the dual leaky bucket algorithm can be applied to a prioritized coding approach to achieve transmission efficiency. A mapping of the performance/control parameters at the network level into equivalent parameters at the video coding level is developed. Based on that, a complete set of principles for the design of video codecs for network transmission is proposed.

Chen, Yun-Chung↗
NASA NTRS2016

The Impact of Truth Surrogate Variance on Quality Assessment/Assurance in Wind Tunnel Testing

Minimum data volume requirements for wind tunnel testing are reviewed and shown to depend on error tolerance, response model complexity, random error variance in the measurement environment, and maximum acceptable levels of inference error risk. Distinctions are made between such related concepts as quality assurance and quality assessment in response surface modeling, as well as between precision and accuracy. Earlier research on the scaling of wind tunnel tests is extended to account for variance in the truth surrogates used at confirmation sites in the design space to validate proposed response models. A model adequacy metric is presented that represents the fraction of the design space within which model predictions can be expected to satisfy prescribed quality specifications. The impact of inference error on the assessment of response model residuals is reviewed. The number of sites where reasonably well-fitted response models actually predict inadequately is shown to be considerably less than the number of sites where residuals are out of tolerance. The significance of such inference error effects on common response model assessment strategies is examined.

DeLoach, Richard↗
NASA NTRS2014

Alma Measurements of the HNC and HC3N Distributions in Titan's Atmosphere

We present spectrally and spatially resolved maps of HNC and HC3N emission from Titan's atmosphere, obtained using the Atacama Large Millimeter/submillimeter Array on 2013 November 17. These maps show anisotropic spatial distributions for both molecules, with resolved emission peaks in Titan's northern and southern hemispheres. TheHC3N maps indicate enhanced concentrations of this molecule over the poles, consistent with previous studies of Titan's photochemistry and atmospheric circulation. Differences between the spectrally integrated flux distributions of HNC and HC3N show that these species are not co-spatial. The observed spectral line shapes are consistent with HNC being concentrated predominantly in the mesosphere and above (at altitudes z approx.. greater than 400 km), whereas HC3N is abundant at a broader range of altitudes (z approx. equal to 70-600 km). From spatial variations in the HC3N line profile, the locations of the HC3N emission peaks are shown to be variable as a function of altitude. The peaks in the integrated emission from HNC and the line core (upper atmosphere) component of HC3N (at z approx. greater than 300 km) are found to be asymmetric with respect to Titan's polar axis, indicating that the mesosphere may be more longitudinally variable than previously thought. The spatially integrated HNC and HC3N spectra are modeled using the NEMESIS planetary atmosphere code and the resulting best-fitting disk-averaged vertical mixing ratio profiles are found to be in reasonable agreement with previous measurements for these species. Vertical column densities of the best-fitting gradient models for HNC and HC3N are 1.9 × 10(exp 13) per sq.cm and 2.3 × 10(exp 14) per sq.cm, respectively.

ALMA Telescope↗
NASA NTRS2019

Quantifying the Influence of Redox Conditions on the Seismic Properties of Olivine

A series of eight olivine specimens were fabricated by hot-pressing at 1200 degrees Celsius and 300 megapascals (MPa). Each hot-pressed specimen was then wrapped in Pt, Ni or NiFe foil to vary oxygen fugacity (fO2), and interrogated via forced torsional oscillation. Mechanical testing was conducted at 10 oscillation periods between 1 and 1000 seconds, at a confining pressure of 200 MPa, during a slow staged-cooling from a maximum temperature of 1200 degrees Celsius to room temperature. After mechanical testing, each specimen was axially sectioned and EBSD (Electron BackScatter Diffraction) was used for the determination of the representative grain size, and grain size distribution of each sample. In addition, each longitudinal section was mapped via FTIR (Fourier Transform InfraRed spectroscopy) to determine the spatial distribution and concentrations of chemically bound and molecular water. Amongst these specimens, chemically bound ‘water’ contents were observed to vary between 0 and 1150 atom ppm (parts per million) H/Si, and molecular water concentrations varied between 0 and 245 atom ppm H/Si. Our forced-oscillation results demonstrate that the measured magnitude of anelastic relaxation within the experimental ‘window’ of oscillation periods is unrelated to the water content. Rather, a relationship was observed between the magnitude of anelastic relaxation and the prevailing redox conditions, which is influenced by the choice of metal sleeving used during the mechanical test. Further, regardless of water content or metal sleeving, each specimen exhibits coupled variations in shear modulus and dissipation within the observational window, indicative of ‘high-temperature background’ behavior, that can be described by a Burgers-type model. During initial fitting of the Burgers models, the unrelaxed shear modulus at a reference temperature of 900 degrees Celsius (G (sub UR)) and the temperature derivative of the unrelaxed shear modulus (dG (sub U)/dT), were treated as adjustable parameters. For all Fe-bearing olivine samples (but not a hydrous and oxidized Fe-free sample) we observe deficits of G (sub UR), and increased values of dG (sub U)/dT, relative to the expected elastic (anharmonic) behavior of Fo90 olivine. This behavior is indicative of anelastic relaxation occurring at shorter periods than observable within the ‘window’ of oscillation periods used in the mechanical test. Moving towards a comprehensive seismologically applicable Burgers model, which includes this newly observed effect of redox conditions on anelastic relaxation, we will present our progress in reconciling truly anharmonic and elastic behavior of Fo (sub 90) olivine with our observed forced-oscillation data.

Cline, Christopher, II↗
NASA NTRS2019

The Influence of Redox Conditions on the Seismic Properties of Polycrystalline Olivine

Eight olivine specimens were fabricated by use of a solgel method and hot-pressing at 1200 degrees Centigrade and 300 megapascals (MPa) inside of welded Pt capsules. Each hot-pressed specimen was then recovered, precision ground, and wrapped in Pt, Ni or NiFe foil to vary oxygen fugacity (fO2) during the subsequent forced torsional oscillation measurements. Mechanical testing was conducted at 10 oscillation periods between 1 and 1000 seconds, at a confining pressure of 200 MPa, during a slow staged-cooling from a maximum temperature of 1200 degrees Centigrade down to room temperature. After mechanical testing, each specimen was axially sectioned and EBSD (Electron BackScatter Diffraction) was used for the determination of the representative grain size, and grain size distribution of each sample. In addition, each longitudinal section was mapped via FTIR (Fourier Transform Infrared Spectroscopy) to determine the spatial distribution and concentrations of chemically bound and molecular water. Amongst these eight specimens, chemically bound ‘water’ contents were observed to vary between 0 and 1150 atom parts per million (ppm) H/Si, and molecular water concentrations varied between 0 and 245 atom ppm H/Si. Our forced-oscillation results demonstrate that the measured magnitude of anelastic relaxation within the experimental ‘window’ of oscillation periods is unrelated to the water content. Rather, a relationship was observed between the magnitude of anelastic relaxation and the prevailing redox conditions, which is influenced by the choice of metal sleeving used during the mechanical test. Further, regardless of water content or metal sleeving, each specimen exhibits coupled variations in shear modulus and dissipation within the observational window, indicative of ‘high-temperature background’ behavior, that can be described by a Burgers-type model. During initial fitting of the Burgers models, the unrelaxed shear modulus at a reference temperature of 900°C (elastic unloading/reloading shear modulus G (sub UR)) and the temperature derivative of the unrelaxed shear modulus (dGU/dT), were treated as adjustable parameters. For all Fe-bearing olivine samples we observe deficits of G (sub UR) and increased values of dGU/dT, relative to the expected elastic (anharmonic) behavior for Fo (Forsterite content percentage) (sub 90) olivine. This behavior is indicative of anelastic relaxation occurring at shorter periods than observable within the ‘window’ of oscillation periods used in the mechanical test. Moving towards a comprehensive and seismologically applicable Burgers model, which includes this newly observed effect of redox conditions on anelastic relaxation, we will present our progress on reconciling truly anharmonic and elastic behavior of Fo (sub 90) olivine with our observed forced-oscillation data.

Cline, C. J., II↗
NASA NTRS2022

Comparison of U.S. and Russian Slow Crack Growth Data and Models

Both the US and Russian space programs use similar predictive models for design of fused silica windows on the International space station. The Russian model can be derived from the power expression for slow crack growth (SCG) or “static fatigue.” The US uses both power and exponential models. Despite the similarity of models and data fitting approach (linear regression and right censoring), different SCG parameters have been derived by US and Russian parties for the same material (Russian fused silica) tested in a similar manner. The difference appears to be related to the use of short-term strength data along with the longer-term static fatigue data, with the power law parameter n being very sensitive to the conversion of strength data into equivalent static data. This hybrid approach is feasible if strength data is measured with a constant stress rate and is appropriately converted to a static equivalent. More research into the approach is needed. However, because of the nonlinear behavior of fused silica in log(v) – log(K(I) ) space and the sensitivity of parameter estimation to fit range, the exponential model is a better choice regardless of test method. Functions are given to convert parameters from the Russian model to those in the US model. However, accurate conversion is hampered by the lack of inert strength data. When the same test technique is used, US and Russian materials exhibit very similar parameters.

Silica, strength, impact, crack growth, windows, I↗
DOE OSTI2021

SANS quantification of bound water in water-soluble polymers across multiple concentration regimes

Contrast-variation small-angle neutron scattering (CV-SANS) is a widely used technique for quantifying hydration water in soft matter systems, but it is predominantly applied in the dilute regime or for systems with a well-defined structure factor. Here, CV-SANS was used to quantify the number of hydration water molecules associating with three water-soluble polymers with different critical solution temperatures and types of water–solute interactions in dilute, semidilute, and concentrated solution through the exploration of novel methods of data fitting and analysis. Multiple SANS fitting workflows with varying levels of model assumptions were evaluated and compared to give insight into SANS model selection. These fitting pathways ranged from general, model-free algorithms to more standard form and structure factor fitting. In addition, Monte Carlo bootstrapping was evaluated as a method to estimate parameter uncertainty through simulation of technical replicates. The most robust fitting workflow for dilute solutions was found to be form factor fitting without CV-SANS ( i.e. polymer in 100% D 2 O). For semidilute and concentrated solutions, while the model-free approach can be mathematically defined for CV-SANS data, the addition of a structure factor imposes physical constraints on the optimization problem, suggesting that the optimal fitting pathway should include appropriate form and structure factor models. The measured hydration numbers were consistent with the number of tightly bound water molecules associated with each monomer unit, and the concentration dependence of the hydration number was largely governed by the chemistry-specific interactions between water and polymer. Polymers with weaker water–polymer interactions ( i.e. those with fewer hydration water molecules) were found to have more bound water at higher concentrations than those with stronger water–polymer interactions due to the increase in the number of forced water–polymer contacts in the concentrated system. This SANS-based method to count hydration water molecules can be applied to polymers in any concentration regime, which will lead to improved understanding of water–polymer interactions and their impact on materials design.

36 MATERIALS SCIENCE↗
NASA NTRS1992

The anomalous extinction curve in the direction of Rho Ophiuchi from 950 to 1180 A

We report the interstellar extinction from 950 to 1180 A in the direction of Rho Oph. These data were obtained with a moderate-resolution sounding rocket observation of Ro Oph, and two lightly reddened B stars of identical spectral type. Unlike Voyager observations in this wavelength region, our spectral resolution is sufficient to directly measure the contribution of molecular hydrogen along the line of sight and unambiguously identify the dust contribution. The derived curves compare favorably with the extrapolation of functional fits to the Rho Oph derived from longer wavelengths. In addition, we have fitted a model of dust extinction to the derived curve and find that the best fits indicate that the line of sight exhibits a larger than average grain size and a preferential concentration of graphite relative to silicates. The sharp rise seen at the shortest wavelengths (below 1000 A) cannot be explained by standard silicate-graphite dust models and may be indicative of a nonstandard grain shape or a missing component in the extinction models.

Green, James C.↗
NASA NTRS1995

The FU Orionis outburst as a thermal accretion event: Observational constraints for protostellar disk models

The results of the time-dependent disk models developed in Bell & Lin are compared with observed properties of FU Orionis variables. Specific models are fit to the light curves of Fu Ori, V1515 Cyg, and V1057 Cyg. The slow risetime of V1515 Cyg can be matched by a self-regulated outburst model. The rapid risetimes of FU Ori and V1057 Cyg can be fitted with the application of modest perturbations to the disk surface density. Model disks display spectral features characteristic of observed objects. The color evolution of V1057 Cyg is naturally explained if mass flux drops in the inner disk (r less than 1/4 AU) while remaining steady in the outer disk. The decrease in optical line width (rotational velocity) observed during the decay of V1057 Cyg may be accounted for by an outward-propagating ionization front. We predict that before final decay to the quiescent phase, short-wavelength line widths (lambda less than 1.5 microns) will again increase. It is suggested that FU Orionis outbursts primarily occur to systems during the embedded phase with ages less than several times 10(exp 5) yr.

Bell, K. R.↗
NASA NTRS1989

Soft X-ray emission lines from low-mass X-ray binaries. I - Collision-dominated models

The soft X-ray line emission from five low-mass X-ray binaries observed by the Einstein Observatory (HEAO 2) Objective Grating Spectrometer is examined. Observed spectra are fitted to model spectra consisting of emission from an optically thin plasma. Fits are performed using a variety of different assumptions for elemental abundances, heating and ionization mechanisms, and the distribution of gas temperature and ionization parameter. The results are discussed in the context of current models for the accretion flow in low-mass X-ray binaries.

Kallman, T. R.↗