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 883 records · Page 49

HERMES impact for the access of Compton form factors

We utilize the DVCS asymmetry measurements of the HERMES collaboration for access to Compton form factors in the deeply virtual regime and to generalized parton distributions. In particular, the (almost) complete measurement of DVCS observables allows us to map various asymmetries into the space of Compton form factors, where we still rely in this analysis on dominance of twist-two associated Compton form factors. We compare this one-to-one map with local Compton form factor fits and a model dependent global fit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

FuelsGen

This toolset generates heterogeneous spatial layouts, aimed at representing forest fuels applicable to wildland fire modeling. Methods are included for calibrating (fitting) this model to data extracted from plot survey information.

Gattiker, James↗

Using heterogeneous data to identify signatures of dengue outbreaks at fine spatio-temporal scales across Brazil

Dengue virus remains a significant public health challenge in Brazil, and seasonal preparation efforts are hindered by variable intra- and interseasonal dynamics. Here, we present a framework for characterizing weekly dengue activity at the Brazilian mesoregion level from 2010–2016 as time series properties that are relevant to forecasting efforts, focusing on outbreak shape, seasonal timing, and pairwise correlations in magnitude and onset. In addition, we use a combination of 18 satellite remote sensing imagery, weather, clinical, mobility, and census data streams and regression methods to identify a parsimonious set of covariates that explain each time series property. The models explained 54% of the variation in outbreak shape, 38% of seasonal onset, 34% of pairwise correlation in outbreak timing, and 11% of pairwise correlation in outbreak magnitude. Regions that have experienced longer periods of drought sensitivity, as captured by the “normalized burn ratio,” experienced less intense outbreaks, while regions with regular fluctuations in relative humidity had less regular seasonal outbreaks. Both the pairwise correlations in outbreak timing and outbreak trend between mesoresgions were best predicted by distance. Our analysis also revealed the presence of distinct geographic clusters where dengue properties tend to be spatially correlated. Forecasting models aimed at predicting the dynamics of dengue activity need to identify the most salient variables capable of contributing to accurate predictions. Our findings show that successful models may need to leverage distinct variables in different locations and be catered to a specific task, such as predicting outbreak magnitude or timing characteristics, to be useful. This advocates in favor of “adaptive models” rather than “one-size-fits-all” models. The results of this study can be applied to improving spatial hierarchical or target-focused forecasting models of dengue activity across Brazil.

59 BASIC BIOLOGICAL SCIENCES↗

Calibration and validation of the foundation for a multiphase strength model for tin

In this work, the Common Model of Multi-phase Strength and Equation of State (CMMP) model was applied to tin. Specifically, calibrations of the strength-specific elements of the CMMP foundation were developed with a combination of experiments and theory, and then the model was validated experimentally. The first element of the foundation is a 10 multi-phase analytic treatment of the melt temperature and the shear modulus for the solid phases. These models were parameterized for each phase based on ab initio calculations using the software VASP (Vienna Ab initio Simulations Package) based on density functional theory (DFT). The shear modulus model for the ambient β phase was validated with ultrasonic sound speed measurements as a function of pressure and temperature. The second element of the foundation is a viscoplastic strength model for the β phase, upon which strength for inaccessible higher-pressure phases can be scaled as necessary. The stress-strain response of tin was measured at strain rates of 10 -3 to 3 x 10 3 s -1 and temperatures ranging from 87 to 373 K. The Preston-Tonks-Wallace (PTW) strength model was fit to that data using Bayesian model calibration. For validation, six forward and two reverse Taylor impact experiments were performed at different velocities to measure large plastic deformation of tin at strain rates up to 10 5 s -1 . The PTW model accurately predicted the deformed shapes of the cylinders, with modest discrepancies attributed to the inability 20 of PTW to capture the effects of twinning and dynamic recrystallization. Some material in the simulations of higher velocity Taylor cylinders reached the melting temperature, thus testing the multiphase model because of the presence of a second phase, the liquid. In simulations using a traditional modeling approach, the abrupt reduction of strength upon melt resulted in poor predictions of the deformed shape and non-physical temperatures. With CMMP, the most deformed material points evolved gradually to a mixed solid-liquid but never fully liquid state, never fully lost strength, 25 stayed at the melt temperature as the latent heat of fusion was absorbed, and predicted the deformed shape well.

36 MATERIALS SCIENCE↗

Overview of Algorithms for Using Particle Morphology in Pre-Detonation Nuclear Forensics

A major goal in pre-detonation nuclear forensics is to infer the processing conditions and/or facility type that produced radiological material. This review paper focuses on analyses of particle size, shape, texture (“morphology”) signatures that could provide information on the provenance of interdicted materials. For example, uranium ore concentrates (UOC or yellowcake) include ammonium diuranate (ADU), ammonium uranyl carbonate (AUC), sodium diuranate (SDU), magnesium diuranate (MDU), and others, each prepared using different salts to precipitate U from solution. Once precipitated, UOCs are often dried and calcined to remove adsorbed water. The products can be allowed to react further, forming uranium oxides UO3, U3O8, or UO2 powders, whose surface morphology can be indicative of precipitation and/or calcination conditions used in their production. This review paper describes statistical issues and approaches in using quantitative analyses of measurements such as particle size and shape to infer production conditions. Statistical topics include multivariate t tests (Hotelling’s T 2 ), design of experiments, and several machine learning (ML) options including decision trees, learning vector quantization neural networks, mixture discriminant analysis, and approximate Bayesian computation (ABC). ABC is emphasized as an attractive option to include the effects of model uncertainty in the selected and fitted forward model used for inferring processing conditions.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

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↗

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↗

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↗

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↗

A TCP model for external beam treatment of intermediate-risk prostate cancer

Purpose: Biological models offer the ability to predict clinical outcomes. The authors describe a model to predict the clinical response of intermediate-risk prostate cancer to external beam radiotherapy for a variety of fractionation regimes. Methods: A fully heterogeneous population averaged tumor control probability model was fit to clinical outcome data for hyper, standard, and hypofractionated treatments. The tumor control probability model was then employed to predict the clinical outcome of extreme hypofractionation regimes, as utilized in stereotactic body radiotherapy. Results: The tumor control probability model achieves an excellent level of fit, R{sup 2} value of 0.93 and a root mean squared error of 1.31%, to the clinical outcome data for hyper, standard, and hypofractionated treatments using realistic values for biological input parameters. Residuals Less-Than-Or-Slanted-Equal-To 1.0% are produced by the tumor control probability model when compared to clinical outcome data for stereotactic body radiotherapy. Conclusions: The authors conclude that this tumor control probability model, used with the optimized radiosensitivity values obtained from the fit, is an appropriate mechanistic model for the analysis and evaluation of external beam RT plans with regard to tumor control for these clinical conditions.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

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↗

Lens models and magnification maps of the six Hubble Frontier Fields clusters

We present strong-lensing models as well as mass and magnification maps for the cores of the six Hubble Space Telescope (HST) Frontier Fields galaxy clusters. Our parametric lens models are constrained by the locations and redshifts of multiple image systems of lensed background galaxies. We use a combination of photometric redshifts and spectroscopic redshifts of the lensed background sources obtained by us (for A2744 and AS1063), collected from the literature, or kindly provided by the lensing community. Using our results, we (1) compare the derived mass distribution of each cluster to its light distribution, (2) quantify the cumulative magnification power of the HST Frontier Fields clusters, (3) describe how our models can be used to estimate the magnification and image multiplicity of lensed background sources at all redshifts and at any position within the cluster cores, and (4) discuss systematic effects and caveats resulting from our modeling methods. We specifically investigate the effect of the use of spectroscopic and photometric redshift constraints on the uncertainties of the resulting models. We find that the photometric redshift estimates of lensed galaxies are generally in excellent agreement with spectroscopic redshifts, where available. However, the flexibility associated with relaxed redshift priors may cause the complexity of large-scale structure that is needed to account for the lensing signal to be underestimated. Our findings thus underline the importance of spectroscopic arc redshifts, or tight photometric redshift constraints, for high precision lens models. All products from our best-fit lens models (magnification, convergence, shear, deflection field) and model simulations for estimating errors are made available via the Mikulski Archive for Space Telescopes.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling resistance to the broadly neutralizing antibody PGT121 in people living with HIV-1

PGT121 is a broadly neutralizing antibody in clinical development for the treatment and prevention of HIV-1 infection via passive administration. PGT121 targets the HIV-1 V3-glycan and demonstrated potent antiviral activity in a phase I clinical trial. Resistance to PGT121 monotherapy rapidly occurred in the majority of participants in this trial with the sampled rebound viruses being entirely resistant to PGT121 mediated neutralization. However, two individuals experienced long-term ART-free viral suppression following antibody infusion and retained sensitivity to PGT121 upon viral rebound. Here, we develop mathematical models of the HIV-1 dynamics during this phase I clinical trial. We utilize these models to understand the dynamics leading to PGT121 resistance and to identify the mechanisms driving the observed long-term viral control. Our modeling highlights the importance of the relative fitness difference between PGT121 sensitive and resistant subpopulations prior to treatment. Specifically, by fitting our models to data, we identify the treatment-induced competitive advantage of previously existing or newly generated resistant population as a primary driver of resistance. Finally, our modeling emphasizes the high neutralization ability of PGT121 in both participants who exhibited long-term viral control.

60 APPLIED LIFE SCIENCES↗

Curve fitting methods for solar radiation data modeling

This paper studies the use of several type of curve fitting method to smooth the global solar radiation data. After the data have been fitted by using curve fitting method, the mathematical model of global solar radiation will be developed. The error measurement was calculated by using goodness-fit statistics such as root mean square error (RMSE) and the value of R{sup 2}. The best fitting methods will be used as a starting point for the construction of mathematical modeling of solar radiation received in Universiti Teknologi PETRONAS (UTP) Malaysia. Numerical results indicated that Gaussian fitting and sine fitting (both with two terms) gives better results as compare with the other fitting methods.

14 SOLAR ENERGY↗

New method for fitting coefficients in standard model effective theory

We present an alternative method for carrying out a principal-component analysis of Wilson coefficients in standard model effective field theory (SMEFT). The method is based on singular-value decomposition (SVD). The SVD method provides information about the sensitivity of experimental observables to physics beyond the standard model that is not accessible in the Fisher-information method. In principle, the SVD method can also have computational advantages over diagonalization of the Fisher information matrix. We demonstrate the SVD method by applying it to the dimension-6 coefficients for the process of top-quark decay to a b quark and a W boson and use this example to illustrate some pitfalls in widely used fitting procedures. We also outline an iterative procedure for applying the SVD method to dimension-8 SMEFT coefficients.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗