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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2019

Complete the Design of Benchtop Fiber-Optic Based Photo-Thermal Radiometry System

Thermal transport properties of nuclear fuels, i.e., thermal conductivity and diffusivity, are critical parameters in the thermal design of a nuclear reactor. It is well known that thermal transport properties change temporally and spatially due to dynamic microstructural evolution as a result of neutron irradiation. A real-time instrument to monitor the local thermal property change, which can survive the harsh in-pile environment, will give researchers unprecedented access to fuel behavior during irradiation. Laser-based photothermal techniques are an ideal candidate to accomplish this specific task, because of the remote, non-contact, and non-destructive features. In this report, we present the progress in the development of a bench-top photothermal radiometry (PTR) system to measure the thermal diffusivity of solid samples, and the preliminary design of a fiber-based PTR instrument to measure the real-time, in-pile thermal diffusivity of nuclear fuels. The PTR technique measures the local temperature response from a modulated heat source by detecting the sample radiation heat flux. Initially a well-developed thermal wave analytical model was fit to the data to obtain the thermal diffusivity. However, systematic errors from the diffraction effect and non-linearity issues resulted in thermal diffusivity measurements higher than literature values. These effects have been studied and an improved analytical model has been developed, and the physics behind the phenomena are still under investigation. A procedure for optimizing measurements from the bench-top system has been provided. The experiences from the bench top measurement were then used to guide the fiber-based system.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2014

DUST AND GAS IN THE MAGELLANIC CLOUDS FROM THE HERITAGE HERSCHEL KEY PROJECT. I. DUST PROPERTIES AND INSIGHTS INTO THE ORIGIN OF THE SUBMILLIMETER EXCESS EMISSION

The dust properties in the Large and Small Magellanic clouds (LMC/SMC) are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 μm. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a power-law emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models, we investigate the origin of the submillimeter excess, defined as the submillimeter emission above that expected from SMBB models fit to observations <200 μm. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 μm submillimeter excesses of 27% and 43% for the LMC and SMC, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our fitting results show that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submillimeter excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 ± 1.7) × 10{sup 5} and (8.3 ± 2.1) × 10{sup 4} M {sub ☉} for the LMC and SMC, respectively. We find significant correlations between the submillimeter excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

THE TEMPERATURES OF RED SUPERGIANTS

We present a re-appraisal of the temperatures of red supergiants (RSGs) using their optical and near-infrared spectral energy distributions (SEDs). We have obtained data of a sample of RSGs in the Magellanic Clouds using VLT+XSHOOTER, and we fit MARCS model atmospheres to different regions of the spectra, deriving effective temperatures for each star from (1) the TiO bands, (2) line-free continuum regions of the SEDs, and (3) the integrated fluxes. We show that the temperatures derived from fits to the TiO bands are systematically lower than the other two methods by several hundred kelvin. The TiO fits also dramatically overpredict the flux in the near-IR, and imply extinctions which are anomalously low compared to neighboring stars. In contrast, the SED temperatures provide good fits to the fluxes at all wavelengths other than the TiO bands, are in agreement with the temperatures from the flux integration method, and imply extinctions consistent with nearby stars. After considering a number of ways to reconcile this discrepancy, we conclude that three-dimensional effects (i.e., granulation) are the most likely cause, as they affect the temperature structure in the upper layers where the TiO lines form. The continuum, however, which forms at much deeper layers, is apparently more robust to such effects. We therefore conclude that RSG temperatures are much warmer than previously thought. We discuss the implications of this result for stellar evolution and supernova progenitors, and provide relations to determine the bolometric luminosities of RSGs from single-band photometry.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Modeling of the blood rheology in steady-state shear flows

We undertake here a systematic study of the rheology of blood in steady-state shear flows. As blood is a complex fluid, the first question that we try to answer is whether, even in steady-state shear flows, we can model it as a rheologically simple fluid, i.e., we can describe its behavior through a constitutive model that involves only local kinematic quantities. Having answered that question positively, we then probe as to which non-Newtonian model best fits available shear stress vs shear-rate literature data. We show that under physiological conditions blood is typically viscoplastic, i.e., it exhibits a yield stress that acts as a minimum threshold for flow. We further show that the Casson model emerges naturally as the best approximation, at least for low and moderate shear-rates. We then develop systematically a parametric dependence of the rheological parameters entering the Casson model on key physiological quantities, such as the red blood cell volume fraction (hematocrit). For the yield stress, we base our description on its critical, percolation-originated nature. Thus, we first determine onset conditions, i.e., the critical threshold value that the hematocrit has to have in order for yield stress to appear. It is shown that this is a function of the concentration of a key red blood cell binding protein, fibrinogen. Then, we establish a parametric dependence as a function of the fibrinogen and the square of the difference of the hematocrit from its critical onset value. Similarly, we provide an expression for the Casson viscosity, in terms of the hematocrit and the temperature. A successful validation of the proposed formula is performed against additional experimental literature data. The proposed expression is anticipated to be useful not only for steady-state blood flow modeling but also as providing the starting point for transient shear, or more general flow modeling.

42 ENGINEERING↗
NASA NTRS1992

The remarkable twisted disk of NGC 4753 and the shapes of galactic halos

The complex dust lanes in the S0 galaxy NGC 4753 are shown to be consistent with a disk that is strongly twisted by differential precession. Yet another peculiar S0 can therefore be explained as the result of an accretion event. An evolving disk model is fitted to the observed distribution. This disk is inclined by 15 deg relative to the galaxy's equatorial plane and twisted such that its line of nodes changes smoothly by 3.8 pi over a factor of seven in radius. The model indicates that most of the galaxy's mass is unseen, is nearly spherically distributed, and has a nearly scale-free spatial distribution. The ellipticity of the total galactic mass distribution must be constant to within 20 percent over the radial extent of the twisted disk - a conclusion which may significantly constrain galaxy formation mechanisms.

Steiman-Cameron, Thomas Y.↗
NASA NTRS1995

Observations and physical interpretations of the solar wind flow properties as obtained from white light coronagraph aboard SPARTAN 201-01

The solar corona was observed with an externally occulted White Light Coronagraph (WLC) carried on the SPARTAN 201-1 spacecraft on 11-12 Apr. 1993. With observations from WLC and the ground based Mauna Loa White Light Coronagraph, a large number of polar plumes both in the north and south polar holes were traced from 1.16 to 5.5 Rs. Flow properties of the solar wind in coronal holes have been determined (Habbal et al., 1995) by using a two fluid model constrained by density profiles and scale height temperatures from the white light observations, and interplanetary measurements of the flow speed and proton mass flux from Ulysses' south polar passage. Provisions for acceleration by Alfven waves, as well as electron and proton heating, are included in the momentum and the energy equations respectively. The model computations fit remarkably well the empirical constraints of the two different density structures (plumes and coronal holes) for a range of input parameters. In this study we investigate the physical nature of the heating function used in the two-fluid model. Alfven waves have been suggested as the possible source of heating that accelerates the solar wind (Ofman and Davila, 1995). We utilize the density contrast observed in WLC data in the plume and ambient coronal hole region to estimate the Alfven wave frequencies responsible for heating these structures. The source heating function utilized in the two fluid model of the solar wind acceleration will be compared with the resonant Alfven wave heating function.

Guhathakurta, Madhulika↗
NASA NTRS2014

Dust and Gas in the Magellanic Clouds from the Heritage Herschel Key Project. I. Dust Properties and Insights into the Origin of the Submm (Submillimeter) Excess Emission

The dust properties in the Large and Small Magellanic Clouds are studied using the HERITAGE Herschel Key Project photometric data in five bands from 100 to 500 micromillimeters. Three simple models of dust emission were fit to the observations: a single temperature blackbody modified by a powerlaw emissivity (SMBB), a single temperature blackbody modified by a broken power-law emissivity (BEMBB), and two blackbodies with different temperatures, both modified by the same power-law emissivity (TTMBB). Using these models we investigate the origin of the submillimeter excess; defined as the submillimeter (submm) emission above that expected from SMBB models fit to observations < 200 micromillimeters. We find that the BEMBB model produces the lowest fit residuals with pixel-averaged 500 micromillimeters submillimeter excesses of 27% and 43% for the Large and Small Magellanic Clouds, respectively. Adopting gas masses from previous works, the gas-to-dust ratios calculated from our fitting results show that the TTMBB fits require significantly more dust than are available even if all the metals present in the interstellar medium (ISM) were condensed into dust. This indicates that the submillimeter excess is more likely to be due to emissivity variations than a second population of colder dust. We derive integrated dust masses of (7.3 plus or minus 1.7) x 10 (sup 5) and (8.3 plus or minus 2.1) x 10 (sup 4) solar masses for the Large and Small Magellanic Clouds, respectively. We find significant correlations between the submillimeter excess and other dust properties; further work is needed to determine the relative contributions of fitting noise and ISM physics to the correlations.

Magellanic↗
NASA NTRS2020

Bayesian Statistical Models for Community Annoyance Survey Data

This paper demonstrates the use of two Bayesian statistical models to analyze single-event sonic boom exposure and human annoyance data from community response surveys. Each model is fit to data from a NASA pilot study.Unlike many community noise surveys, this study used a panel sample to collect multiple observations per participant instead of a single observation. Thus, a multilevel (also known as hierarchical or mixed-effects) model is used to account for the within-subject correlation in the panel sample data. This paper describes a multilevel logistic regression model and a multilevel ordinal regression model. The paper also proposes a method for calculating a summary dose-response curve from the multilevel models that represents the population. The two models’ summary dose-response curves are visually similar. However, their estimates differ when calculating the noise dose at a fixed percent highly annoyed.

Musical instruments↗
DOE OSTI2015

THE CIRCUMSTELLAR DISK OF THE Be STAR o AQUARII AS CONSTRAINED BY SIMULTANEOUS SPECTROSCOPY AND OPTICAL INTERFEROMETRY

Omicron Aquarii is a late-type, Be shell star with a stable and nearly symmetric Hα emission line. We combine Hα interferometric observations obtained with the Navy Precision Optical Interferometer covering 2007 through 2014 with Hα spectroscopic observations over the same period and a 2008 observation of the system's near-infrared spectral energy distribution to constrain the properties of o Aqr's circumstellar disk. All observations are consistent with a circumstellar disk seen at an inclination of 75° ± 3° with a position angle on the sky of 110° ± 8° measured East from North. From the best-fit disk density model, we find that 90% of the Hα emission arises from within 9.5 stellar radii, and the mass associated with this Hα disk is ∼1.8 × 10{sup −10} of the stellar mass, and that the associated angular momentum, assuming Keplerian rotation for the disk, is ∼1.6 × 10{sup −8} of the total stellar angular momentum. The occurrence of a central quasi-emission feature in Mg ii λ4481 is also predicted by this best-fit disk model and the computed profile compares successfully with observations from 1999. To obtain consistency between the Hα line profile modeling and the other constraints, it was necessary in the profile fitting to weight the line core (emission peaks and central depression) more heavily than the line wings, which were not well reproduced by our models. This may reflect the limitation of assuming a single power law for the disk's variation in equatorial density. The best-fit disk density model for o Aqr predicts that Hα is near its maximum strength as a function of disk density, and hence the Hα equivalent width and line profile change only weakly in response to large (factor of ∼5) changes in the disk density. This may in part explain the remarkable observed stability of o Aqr's Hα emission line profile.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Predicting molecular dipole moments by combining atomic partial charges and atomic dipoles

The molecular dipole moment ( μ ) is a central quantity in chemistry. It is essential in predicting infrared and sum-frequency generation spectra as well as induction and long-range electrostatic interactions. Furthermore, it can be extracted directly—via the ground state electron density—from high-level quantum mechanical calculations, making it an ideal target for machine learning (ML). Here, we choose to represent this quantity with a physically inspired ML model that captures two distinct physical effects: local atomic polarization is captured within the symmetry-adapted Gaussian process regression framework which assigns a (vector) dipole moment to each atom, while the movement of charge across the entire molecule is captured by assigning a partial (scalar) charge to each atom. The resulting “MuML” models are fitted together to reproduce molecular μ computed using high-level coupled-cluster theory and density functional theory (DFT) on the QM7b dataset, achieving more accurate results due to the physics-based combination of these complementary terms. The combined model shows excellent transferability when applied to a showcase dataset of larger and more complex molecules, approaching the accuracy of DFT at a small fraction of the computational cost. We also demonstrate that the uncertainty in the predictions can be estimated reliably using a calibrated committee model. The ultimate performance of the models—and the optimal weighting of their combination—depends, however, on the details of the system at hand, with the scalar model being clearly superior when describing large molecules whose dipole is almost entirely generated by charge separation. These observations point to the importance of simultaneously accounting for the local and non-local effects that contribute to μ ; furthermore, they define a challenging task to benchmark future models, particularly those aimed at the description of condensed phases.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1990

Universal adsorption at the vapor-liquid interface near the consolute point

The ellipticity of the vapor-liquid interface above mixtures of methylcyclohexane (C7H14) and perfluoromethylcyclohexane (C7F14) has been measured near the consolute point T(c) = 318.6 K. The data are consistent with a model of the interface that combines a short-ranged density-vs height profile in the vapor phase with a much longer-ranged composition-versus-height profile in the liquid. The value of the free parameter produced by fitting the model to the data is consistent with results from two other simple mixtures and a mixture of a polymer and solvent. This experiment combines precision ellipsometry of the vapor-liquid interface with in situ measurements of refractive indices of the liquid phases, and it precisely locates the consolute point.

Schmidt, James W.↗
NASA NTRS1994

Ion transport and loss in the Earth's quiet ring current. 2: Diffusion and magnetosphere-ionosphere coupling

We have studied the transport and loss of H(+), He(+), and He(++) ions in the Earth's quiet time ring current (1 to 300 keV/e, 3 to 7 R(sub E), Kp less than 2+, absolute value of Dst less than 11, 70 to 110 degs pitchangles, all LT) comparing the standard radial diffusion model developed for the higher-energy radiation belt particles with measurements of the lower energy ring current ions in a previous paper. Large deviations of that model, which fit only 50% of the data to within a factor of 10, suggested that another transport mechanism is operating in the ring current. Here we derive a modified diffusion coefficient corrected for electric field effects on ring current energy ions that fit nearly 80% of the data to within a factor of 2. Thus we infer that electric field fluctuations from the low-latitude to midlatitude ionosphere (ionospheric dynamo) dominated the ring current transport, rather than high-latitude or solar wind fluctuations. Much of the remaining deviation may arise from convective electric field transport of the E less than 30 keV particles. Since convection effects cannot be correctly treated with this azimuthally symmetric model, we defer treatment of the lowest-energy ions to a another paper. We give chi(exp 2) contours for the best fit, showing the dependence of the fit upon the internal/external spectral power of the predicted electric and magnetic field fluctuations.

Sheldon, R. B.↗
DOE OSTI2016

Rotation of X-ray polarization in the glitches of a silicon crystal monochromator

EXAFS studies on dilute samples are usually carried out by collecting the fluorescence yield using a large-area multi-element detector. This method is susceptible to the `glitches' produced by all single-crystal monochromators. Glitches are sharp dips or spikes in the diffracted intensity at specific crystal orientations. If incorrectly compensated, they degrade the spectroscopic data. Normalization of the fluorescence signal by the incident flux alone is sometimes insufficient to compensate for the glitches. Measurements performed at the state-of-the-art wiggler beamline I20-scanning at Diamond Light Source have shown that the glitches alter the spatial distribution of the sample's quasi-elastic X-ray scattering. Because glitches result from additional Bragg reflections, multiple-beam dynamical diffraction theory is necessary to understand their effects. Here, the glitches of the Si(111) four-bounce monochromator of I20-scanning just above the Ni K edge are associated with their Bragg reflections. A fitting procedure that treats coherent and Compton scattering is developed and applied to a sample of an extremely dilute (100 micromolal) aqueous solution of Ni(NO 3 ) 2 . The depolarization of the wiggler X-ray beam out of the electron orbit is modeled. The fits achieve good agreement with the sample's quasi-elastic scattering with just a few parameters. The X-ray polarization is rotated up to ±4.3° within the glitches, as predicted by dynamical diffraction. These results will help users normalize EXAFS data at glitches.

Sutter, John P.↗
NASA NTRS2002

A Hall Thruster Performance Model Incorporating the Effects of a Multiply-Charged Plasma

A Hall thruster performance model that predicts anode specific impulse, anode efficiency, and thrust is discussed. The model is derived as a function of a voltage loss parameter, an electron loss parameter, and the charge state of the plasma. Experimental data from SPT and TAL type thrusters up to discharge powers of 21.6 kW are used to determine the best fit for model parameters. General values for the model parameters are found, applicable to high power thrusters and irrespective of thruster type. Performance of a 50 kW thruster is calculated for an anode specific impulse of 2500 seconds or a discharge current of 100 A.

Hofer, Richard R.↗
DOE OSTI2015

On the pressure effect in energetic deposition of Cu thin films by modulated pulsed power magnetron sputtering: A global plasma model and experiments

The modulated pulsed power magnetron sputtering (MPPMS) discharge processes are numerically modeled and experimentally investigated, in order to explore the effect of the pressure on MPPMS discharges as well as on the microstructure of the deposited thin films. A global plasma model has been developed based on a volume-averaged global description of the ionization region, considering the loss of electrons by cross-B diffusion. The temporal variations of internal plasma parameters at different pressures from 0.1 to 0.7 Pa are obtained by fitting the model to duplicate the experimental discharge data, and Cu thin films are deposited by MPPMS at the corresponding pressures. The surface morphology, grain size and orientation, and microstructure of the deposited thin films are investigated by scanning electron microscopy, transmission electron microscopy, and x-ray diffraction. By increasing the pressure from 0.1 to 0.7 Pa, both the ion bombardment energy and substrate temperature which are estimated by the modeled plasma parameters decrease, corresponding to the observed transition of the deposited thin films from a void free structure with a wide distribution of grain size (zone T) into an underdense structure with a fine fiber texture (zone 1) in the extended structure zone diagram (SZD). The microstructure and texture transition of Cu thin films are well-explained by the extended SZD, suggesting that the primary plasma processes are properly incorporated in the model. The results contribute to the understanding of the characteristics of MPPMS discharges, as well as its correlation with the microstructure and texture of deposited Cu thin films.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

Modelling individual variability in growth of bull trout in the Walla Walla River Basin using a hierarchical von Bertalanffy growth model

Abstract We examined growth in length of fluvial bull trout ( Salvelinus confluentus ) in the Walla Walla River Basin, Washington and Oregon. Our objectives were to quantify individual variability in growth; examine growth within and among years, life history forms, life stages and sexes; and estimate von Bertalanffy growth parameters. Individual variability was evaluated by modelling asymptotic length ( L ∞ ) and the growth coefficient ( k ) as random variables. All models were fit with Bayesian methods and were evaluated for fit by the deviance information criterion. By incorporating individual variability, population‐level estimates of L ∞ and k appeared appropriate and estimated growth trajectories for specific bull trout fit individual observed patterns in growth. Growth trajectories and positive correlation between individual estimates of L ∞ and k suggest that some individuals grow at a faster rate and reach a larger maximum size than other individuals and those differences are maintained throughout life. Selected models suggest that fluvial migrants have higher estimates of L ∞ and k than residents, but there were only slight differences in parameter estimates among migrants from two adjacent spawning populations in the Walla Walla River Basin, as well as between males and females. Growth rates increased for fluvial migrants after subadult emigration. Individual variability in growth is consistent with the life history diversity assumed essential for bull trout population persistence. Quantifying this variability is important for modelling population dynamics and viability to conserve this threatened species.

Harris, Julianne E.↗
NASA NTRS1982

IUE spectroscopy, visible-band polarimetry, and radiometry of V641 Mon

This hot, double line, ellipsoidal variable member of NGC 2264 has been shown previously to be either a semi-detached or contact close binary. Low-resolution IUE spectra are best fitted to a Kurucz model atmosphere for very small (approximately 0.08 mag) E(B-V). The familiar interstellar absorption dip near lambda 2200 is apparently absent. A suitable model atmosphere can be fitted to the IUE fluxes, but flux excesses (compared to the model) appear for all the published U through L magnitudes. The spectrum of the B through L excess appears to follow a .0001 lambda dependence. It is shown that this cannot be interpreted as arising from another star fortuitously observed in the visible band or IR. Ground based polarization measures indicate V641 Mon to be a polarization variable. Previous and new V light curves show the amplitude of light variability itself to be variable by about a factor of 2. It is suggested that all these observed characteristics are best explained by postulating "third light" and identifying part of it with Rayleigh scattered starlight very near to the stars. From this same region there arise circumstellar absorptions which give rise to nontheoretical strengths for Si II and Si III lines.

Koch, R. H.↗
DOE OSTI2026

Reformulating a Gurson-based dynamic damage model and demonstrating improved predictive power and numerical robustness

The original Tepla (TEnsile PLAsticity) ductile damage model, based on the Gurson yield surface, has long been used to model damage evolution and material failure under dynamic loading. Unfortunately, Tepla suffered from mesh sensitivity, numerical instability, and limited predictive capability. Here, we theoretically reformulate Tepla to address these issues. We especially focus on the prediction of porosity, which is the key state variable used for modeling ductile damage, as compared to more easily measured surface velocities, which are at best an indirect measure of damage. Key model changes include separating the viscosity during volumetric void growth from underlying shear strength behavior and switching to an iterative bisection solver. The new Tepla is then calibrated on incipient spall experiments on half-hard copper and tantalum, which demonstrate its ability to simultaneously fit the model to recovered porosity distributions and measured surface velocities, a stringent test. Improved numerical behavior, such as greatly reduced mesh sensitivity, is also shown in those simulations. Finally, the new Tepla model is applied to several high-explosive loaded, sweeping wave experiments, showing the ability of the model to predict behavior on tests with significantly different loading conditions and histories than the calibration data.

36 MATERIALS SCIENCE↗