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,297 records · Page 72

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

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

Hunt, R. H.↗

The structure of clusters of galaxies observed with Einstein

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

Jones, C.↗

Extended atmospheres of comets and outer planet-satellite systems

For the hydrogen coma of comet P/Halley, both a Lyman-alpha image and extensive Lyman-alpha scan data obtained by the Pioneer Venus Orbiter Ultraviolet spectrometer as well as H-alpha ground-based spectral observations obtained by the University of Wisconsin Space Physics Group were successfully interpreted and analyzed with our Monte Carlo particle trajectory model. The excellent fit of the model and the Halley data and the water production rate determined near perihelion (9 Feb. 1986) from 13 Dec. 1985 to 13 Jan. 1986 and from 1 Feb. to 7 Mar. 1986 are discussed. Studies for the circumplanetary distribution of atomic hydrogen in the Saturn and Neptune systems were undertaken for escape of H atoms from Titan and Triton, respectively. The discovery of a new mechanism which can dramatically change the normal cylindrically symmetric distribution of hydrogen about the planet is discussed. The implications for the Titan-Saturn and Triton-Neptune are summarized.

Smyth, William H.↗

First VLTI/GRAVITY Observations of HIP 65426 b: Evidence for a Low or Moderate Orbital Eccentricity

Giant exoplanets have been directly imaged over orders of magnitude of orbital separations, prompting theoretical and observational investigations of their formation pathways. In this paper, we present new VLTI/GRAVITY astrometric data of HIP 65426 b, a cold, giant exoplanet which is a particular challenge for most formation theories at a projected separation of 92 au from its primary. Leveraging GRAVITY's astrometric precision, we present an updated eccentricity posterior that disfavors large eccentricities. The eccentricity posterior is still prior dependent, and we extensively interpret and discuss the limits of the posterior constraints presented here. We also perform updated spectral comparisons with self-consistent forward-modeled spectra, finding a best-fit ExoREM model with solar metallicity and C/O = 0.6. An important caveat is that it is difficult to estimate robust errors on these values, which are subject to interpolation errors as well as potentially missing model physics. Taken together, the orbital and atmospheric constraints paint a preliminary picture of formation inconsistent with scattering after disk dispersal. Further work is needed to validate this interpretation. Analysis code used to perform this work is available on GitHub: https://github.com/sblunt/hip65426.

Sarah Blunt↗

The impact of alternative land and yield assumptions in herbaceous biomass supply modeling: one-size-fits-all resource assessment?: One-size-fits-all resource assessment?

The Billion-ton Reports series has addressed the technical economic potential of supplying additional biomass from farmland and forests.Underlying each of the reports and supporting scenarios is a series of assumptions that drive the modeled output. The assumptions have developed over time with the support of technical experts from industry, academia, and government. Energy crops have not yet reached commodity scale, and only exist in commercial production in a limited number of regions, so many of the key assumptions of biomass potential are not supported by commercial demonstrations across multiple regions. This report examines the biomass potential of herbaceous energy crops from US farmland when alternative assumptions concerning production are modeled. In particular, the assumption of energy grasses, both perennial and annual, is modified with respect to land eligibility, annual rotations, and yield growth. The results of these changes are presented in six scenarios – three of which are modifications to the base-case scenario and three to the high-yield scenario. Finally, the impact on land enrolled and biomass supply are discussed with respect to the most recent 2016 Billion-ton Report results.

09 BIOMASS FUELS↗

Surface fitting with biharmonic and harmonic models

The physical and geometric interpretation of a surface fitting technique, multiquadric equations (MQ), is discussed. The relation of this method and its reciprocal form to biharmonic and harmonic techniques respectively, is pointed out. In both forms the solutions may be viewed as being numerical approximations of an integral equation in which an unknown density function is the physical source for disturbing potential (reciprocal-MQ) or elastic displacement (MQ).

Hardy, R. L.↗

Use of Sobol’ Variance-Based Global Sensitivity Analysis and Multidimensional Legendre Polynomial Fitting for Reduced Order Modeling

Sandia National Laboratories (SNL) has developed a novel reduced order modeling approach. Prioritization of inputs is accomplished using Sobo' indices obtained through a more efficient variance-based global sensitivity analysis. To determine the Sobo' functions, simulated input values are aligned to collocation points to permit the use of Gauss-Lobatto integration, thereby reducing the number of simulation trials needed by more than an order of magnitude compared to standard Monte Carlo approaches. Furthermore, by leveraging the orthogonality of Legendre polynomials in conjunction with those same simulations at the collocation nodes, an efficient fitting method is developed to represent the Sobo' functions from which a reduced order model (ROM) is constructed. The developed method is both more efficient computationally, and the resulting ROM is more accurate. The efficacy of this technique is demonstrated on a nonlinear polynomial test function as well as the nonlinear Ishigami and Sobo' g functions.

97 MATHEMATICS AND COMPUTING↗

Limits in late time conversion of cold dark matter into dark radiation

Structure formation creates high temperature and density regions in the Universe that allow the conversion of matter into more stable states, with a corresponding emission of relativistic matter and radiation. An example of such a mechanism is the supernova event, that releases relativistic neutrinos corresponding to 99% of the binding energy of remnant neutron star. We take this phenomena as a starting point for an assumption that similar processes could occur in the dark sector, where structure formation would generate a late time conversion of cold dark matter into a relativistic form of dark matter. We performed a phenomenological study about the limits of this conversion, where we assumed a transition profile that is a generalized version of the neutrino production in supernovae events. With this assumption, we obtained an interesting modification for the constraint over the cold dark matter density. We show that when comparing with the standard ΛCDM cosmology, there is no preference for conversion, although the best fit is within 1σ from the standard model best fit. The methodology and the results obtained qualify this conversion hypothesis, from the large scale structure point of view, as a viable and interesting model to be tested in the future with small scale data, and mitigate discrepancies between observations at this scale and the pure cold dark matter model.

79 ASTRONOMY AND ASTROPHYSICS↗

Chrono-Validation of Near-Real-Time Landslide Susceptibility Models via Plugin Statistical Simulations

The idea behind any validation scheme in landslide susceptibility studies is to test whether a model calibrated on a certain data can predict an unknown dataset of the same nature (landslide presences/absences and covariates). Almost the entirety of landslide susceptibility studies are validated by subsetting a single dataset into a training and test sets. This dataset usually corresponds either to event-specific or to historical inventories. Very rarely, a multi-temporal inventory is available and, in the few cases where this condition is met, the validation practices involve training a model on a specific landslide inventory, deriving a single predictive equation and validating it on a subsequent landslide inventory. This commonly leads landslide predictive studies, even those with a strong statistical rigor, to neglect the uncertainty estimation in their modeling scheme. In statistics, validation can also be performed via statistical simulations. This means that after fitting a given model, one can generate any number of predictive functions and test their predictive skills on any type and number of unknown datasets. In this work, we take a similar direction and we apply it to model and validate three separate co-seismic inventories, including an uncertainty estimation phase. We mapped these inventories within the same area in Indonesia, for three earthquakes occurred in 2012, 2017 and 2018. Specifically, we build three event-specific Bayesian Generalize Additive Models of the binomial family. From each model we then simulate 1000 predictive realizations over the remaining two inventories, by using a plug-in scheme where all the morphometric covariates are kept fixed and only the ground motion is replaced according to the prediction target. By doing so, we introduce a new analytical tool for near-real-time landslide predictive purposes, which is able to produce a probabilistic model which stands in between the definitions of susceptibility and hazard. In fact, our model is able to accurately estimate “where” and “when” - although not “how frequently” - landslide have occurred by featuring the multitemporal information of the trigger. In our findings, the simulations are quite similar to the fitted models; and the nine combinations we analyse produce excellent performance. This result confirms the assumption that “the past is the key to the future”, as we show that the relative contribution of each variable and their interactions in each probabilistic model remains practically the same across temporal replicates. This information is not trivial because it supports the routines implemented in global near-real-time applications.

Temporal validation↗

Numerical evaluation of photosensitive tracers as a strategy for separating surface and subsurface transient storage in streams

For this work, we numerically evaluated photosensitive tracers as a potential strategy for separating the effects of surface and hyporheic storage zones (SSZs and HSZs, respectively) on stream corridor transport. Correctly separating HSZ and SSZ effects is critical to estimating the hydro-biogeochemical function of a stream because HSZs and SSZs expose solutes to significantly different biogeochemical conditions, like sunlight exposure, microbial processes, and oxygen availability. Our numerical experiments used a multiscale river-corridor transport model implemented in the ATS code, which accommodates multiple storage zones with distinct travel time distributions and biogeochemical reactions. For parameter inferences, we used Bayesian inverse modeling. We found that breakthrough curves for photo-decaying tracers from day and night injection can delineate surface and hyporheic transient storage contributions, but only when interpreted jointly through a two-storage zone model. Numerical experiments that used only daytime injection or interpreted breakthrough curves with a single storage zone model yielded good fit to breakthrough curves, but parameter estimates were biased and controlling processes misattributed, examples of good model fits for the wrong reasons. Using those biased parameter estimates in reactive transport simulations resulted in significantly different projections of denitrification, which underscores the potential for stream function to be mischaracterized if tracer tests are interpreted through an inappropriately simplified model for transient storage. More generally, this study highlights the role of modeling in evaluating the experimental design and identifying the potential of system mischaracterization and its implications.

54 ENVIRONMENTAL SCIENCES↗

Surface daytime net radiation estimation using artificial neural networks

Net all-wave surface radiation (R n ) is one of the most important fundamental parameters in various applications. However, conventional R n measurements are difficult to collect because of the high cost and ongoing maintenance of recording instruments. Therefore, various empirical R n estimation models have been developed. This study presents the results of two artificial neural network (ANN) models (general regression neural networks (GRNN) and Neuroet) to estimate R n globally from multi-source data, including remotely sensed products, surface measurements, and meteorological reanalysis products. R n estimates provided by the two ANNs were tested against in-situ radiation measurements obtained from 251 global sites between 1991–2010 both in global mode (all data were used to fit the models) and in conditional mode (the data were divided into four subsets and the models were fitted separately). Based on the results obtained from extensive experiments, it has been proved that the two ANNs were superior to linear-based empirical models in both global and conditional modes and that the GRNN performed better and was more stable than Neuroet. The GRNN estimates had a determination coefficient (R 2 ) of 0.92, a root mean square error (RMSE) of 34.27 W·m –2 , and a bias of –0.61 W·m –2 in global mode based on the validation dataset. In conclusion, ANN methods are a potentially powerful tool for global R n estimation.

14 SOLAR ENERGY↗

Comparison of Echo 7 field line length measurements to magnetospheric model predictions

The Echo 7 sounding rocket experiment injected electron beams on central tail field lines near L = 6.5. Numerous injections returned to the payload as 'conjugate echoes' after mirroring in the southern hemisphere. Field line lengths calculated from measured conjugate echo bounce times and energies are compared to predictions made by integrating electron trajectories through various magnetospheric models. Although Kp at launch was 3-, quiet time magnetic models best fit the echo measurements. Geosynchronous satellite magnetometer measurements near the Echo 7 field lines during the flight were best modeled by the Olson-Pfitzer Dynamic Model and the Tsyganenko-Usmanov model for Kp = 3. The discrepancy between the models that best fit the Echo 7 data and those that fit the satellite data was most likely due to uncertainties in the small-scale configuration of the magnetospheric models. The field line length measured by the conjugate echoes showed some temporal variation in the magnetic field, also indicated by the satellite magnetometers. This demonstrates the utility an Echo-style experiment could have in substorm studies.

Nemzek, R. J.↗

ASSESSMENT OF MODELS OF GALACTIC THERMAL DUST EMISSION USING COBE /FIRAS AND COBE /DIRBE OBSERVATIONS

Accurate modeling of the spectrum of thermal dust emission at millimeter wavelengths is important for improving the accuracy of foreground subtraction for cosmic microwave background (CMB) measurements, for improving the accuracy with which the contributions of different foreground emission components can be determined, and for improving our understanding of dust composition and dust physics. We fit four models of dust emission to high Galactic latitude COBE /FIRAS and COBE /DIRBE observations from 3 mm to 100 μ m and compare the quality of the fits. We consider the two-level systems (TLS) model because it provides a physically motivated explanation for the observed long wavelength flattening of the dust spectrum and the anti-correlation between emissivity index and dust temperature. We consider the model of Finkbeiner et al. because it has been widely used for CMB studies, and the generalized version of this model that was recently applied to Planck data by Meisner and Finkbeiner. For comparison we have also fit a phenomenological model consisting of the sum of two graybody components. We find that the two-graybody model gives the best fit and the FDS model gives a significantly poorer fit than the other models. The Meisner and Finkbeiner model and the TLS model remain viable for use in Galactic foreground subtraction, but the FIRAS data do not have a sufficient signal-to-noise ratio to provide a strong test of the predicted spectrum at millimeter wavelengths.

79 ASTRONOMY AND ASTROPHYSICS↗

Local and Global Function Model of the Liver

Purpose: To develop a local and global function model in the liver based on regional and organ function measurements to support individualized adaptive radiation therapy (RT). Methods and Materials: A local and global model for liver function was developed to include both functional volume and the effect of functional variation of subunits. Adopting the assumption of parallel architecture in the liver, the global function was composed of a sum of local function probabilities of subunits, varying between 0 and 1. The model was fit to 59 datasets of liver regional and organ function measures from 23 patients obtained before, during, and 1 month after RT. The local function probabilities of subunits were modeled by a sigmoid function in relating to MRI-derived portal venous perfusion values. The global function was fitted to a logarithm of an indocyanine green retention rate at 15 minutes (an overall liver function measure). Cross-validation was performed by leave-m-out tests. The model was further evaluated by fitting to the data divided according to whether the patients had hepatocellular carcinoma (HCC) or not. Results: The liver function model showed that (1) a perfusion value of 68.6 mL/(100 g · min) yielded a local function probability of 0.5; (2) the probability reached 0.9 at a perfusion value of 98 mL/(100 g · min); and (3) at a probability of 0.03 [corresponding perfusion of 38 mL/(100 g · min)] or lower, the contribution to global function was lost. Cross-validations showed that the model parameters were stable. The model fitted to the data from the patients with HCC indicated that the same amount of portal venous perfusion was translated into less local function probability than in the patients with non-HCC tumors. Conclusions: The developed liver function model could provide a means to better assess individual and regional dose-responses of hepatic functions, and provide guidance for individualized treatment planning of RT.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Dynamic modeling of fluid transmission lines of the DSN 70-meter antennas by using a lumped parameter model

Fluid transmission lines and fittings were found to significantly affect the dynamic and steady state operation of the engineering system in which they are incorporated. Therefore, a better understanding of the operation of a system can be obtained by including the performance of the transmission lines and fittings within the system model. The most accurate model of a transmission line or fitting is obtained by using a distributed parameter model. However, a distributed parameter model tends to be very complex. This complexity can be avoided without significant loss of model accuracy by using a lumped parameter model when the length of the fluid path through the transmission line or fitting is short. This article develops a lumped parameter model for short fluid transmission lines and fittings, describes the conditions under which the model is valid, and presents the model parameters associated with the servo hydraulic system of the Deep Space Network (DSN) 70-meter antennas.

Bartos, R. D.↗

An Analysis of Dust Halo and Extinction Toward X Persei

Interstellar dust grain models are not sufficiently constrained by UV extinction curves to be able to distinguish between them. By testing grain models in the X-ray regime and applying elemental abundance constraints, we show to what extent the models can reproduce the observables in these regimes, and if they are capable of doing so while respecting the abundance limits. We tested the MRN and WD grain models. The fits to the X-ray data do not allow us to distinguish between MRN and WD; both models provide reasonable fits, but cannot do so while respecting the elemental abundance constraints. The situation in the UV regime is similar. Both MRN and WD underestimate the hydrogen column density NH. The model of ZDA provides promising results, as it finds NH much closer to the UV-measured value; further testing of this model is called for.

Valencic, Lynne A.↗

Assessment of Models of Galactic Thermal Dust Emission Using COBE/FIRAS and COBE/DIRBE Observations

Accurate modeling of the spectrum of thermal dust emission at millimeter wavelengths is important for improving the accuracy of foreground subtraction for cosmic microwave background (CMB) measurements, for improving the accuracy with which the contributions of different foreground emission components can be determined, and for improving our understanding of dust composition and dust physics. We fit four models of dust emission to high Galactic latitude COBE/FIRAS and COBE/DIRBE observations from 3 mm to 100m and compare the quality of the fits. We consider the two-level systems (TLS) model because it provides a physically motivated explanation for the observed long wavelength flattening of the dust spectrum and the anti-correlation between emissivity index and dust temperature. We consider the model of Finkbeiner et al. because it has been widely used for CMB studies, and the generalized version of this model that was recently applied to Planck data by Meisner and Finkbeiner. For comparison we have also fit a phenomenological model consisting of the sum of two-graybody components. We find that the two-graybody model gives the best fit and the FDS model gives a significantly poorer fit than the othermodels. The Meisner and Finkbeiner model and the TLS model remain viable for use in Galactic foreground subtraction, but the FIRAS data do not have a sufficient signal-to-noise ratio to provide a strong test of the predicted spectrum at millimeter wavelengths.

Odegard, N.↗

Statistical Emulation of Climate Model Projections Based on Precomputed GCM Runs*

The authors describe a new approach for emulating the output of a fully coupled climate model under arbitrary forcing scenarios that is based on a small set of precomputed runs from the model. Temperature and precipitation are expressed as simple functions of the past trajectory of atmospheric CO 2 concentrations, and a statistical model is fit using a limited set of training runs. The approach is demonstrated to be a useful and computationally efficient alternative to pattern scaling and captures the nonlinear evolution of spatial patterns of climate anomalies inherent in transient climates. The approach does as well as pattern scaling in all circumstances and substantially better in many; it is not computationally demanding; and, once the statistical model is fit, it produces emulated climate output effectively instantaneously. In conclusion, it may therefore find wide application in climate impacts assessments and other policy analyses requiring rapid climate projections.

54 ENVIRONMENTAL SCIENCES↗