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 955 records · Page 53

Measuring neutrino masses with a future galaxy survey

We perform a detailed forecast on how well a EUCLID-like photometric galaxy and cosmic shear survey will be able to constrain the absolute neutrino mass scale. Adopting conservative assumptions about the survey specifications and assuming complete ignorance of the galaxy bias, we estimate that the minimum mass sum of Σm{sub ν} ≅ 0.06 eV in the normal hierarchy can be detected at 1.5σ to 2.5σ significance, depending on the model complexity, using a combination of galaxy and cosmic shear power spectrum measurements in conjunction with CMB temperature and polarisation observations from PLANCK. With better knowledge of the galaxy bias, the significance of the detection could potentially reach 5.4σ. Interestingly, neither PLANCK+shear nor PLANCK+galaxy alone can achieve this level of sensitivity; it is the combined effect of galaxy and cosmic shear power spectrum measurements that breaks the persistent degeneracies between the neutrino mass, the physical matter density, and the Hubble parameter. Notwithstanding this remarkable sensitivity to Σm{sub ν}, EUCLID-like shear and galaxy data will not be sensitive to the exact mass spectrum of the neutrino sector; no significant bias ( < 1σ) in the parameter estimation is induced by fitting inaccurate models of the neutrino mass splittings to the mock data, nor does the goodness-of-fit of these models suffer any significant degradation relative to the true one (Δχ{sub eff}{sup 2} < 1)

79 ASTRONOMY AND ASTROPHYSICS↗

Contrasts in compliant fault zone properties inferred from geodetic measurements in the San Francisco Bay area

In crustal fault zones, regions of damaged rock characterized by reduced elastic shear modulus can influence patterns of near-field interseismic deformation. In order to study these compliant fault zones (CFZs) and how they might develop over the lifetimes of faults, we compare two fault segments with contrasting fault age and lithology along the San Andreas Fault in the San Francisco Bay Area. New geodetic measurements of the interseismic velocity fields at each location are used to constrain fault zone parameters through a Markov chain Monte Carlo method. At Black Mountain, in the Santa Cruz Mountains of the San Francisco Peninsula, we do not find evidence for a compliant fault zone; instead, we find that the geodetic data are more consistent with a model of a single fault in a homogeneous elastic half-space. At Lake San Andreas, a younger fault segment 35 km farther north, we find evidence for a compliant fault zone about 3.4 +1.1/-1.4 km wide, containing a shear modulus of about 40% of the shear modulus of the surrounding rock. We also find that the best fitting CFZ model at this location, unlike the best fitting homogeneous half-space model, has a locking depth that agrees well with the observed depth of microseismicity. Based on differences in fault age, cumulative displacement, and lithology between Black Mountain and Lake San Andreas, we infer that lithology plays an important and, in this case, perhaps a dominant role in the accumulation of fault zone damage structures and the development of CFZs over the lifetime of a fault.

58 GEOSCIENCES↗

Scalable computations for nonstationary Gaussian processes

Nonstationary Gaussian process models can capture complex spatially varying dependence structures in spatial datasets. However, the large number of observations in modern datasets makes fitting such models computationally intractable with conventional dense linear algebra. In addition, derivative-free or even first-order optimization methods can be very slow to converge when estimating many spatially varying parameters. In this paper, we present a computational framework which couples an algebraic block diagonal plus low-rank covariance matrix approximation with stochastic trace estimation to facilitate the efficient use of second-order solvers for maximum likelihood estimation of Gaussian process models with many parameters. We demonstrate the effectiveness of these methods by simultaneously fitting 192 parameters in the popular nonstationary model of Paciorek and Schervish using 107,600 sea surface temperature anomaly measurements.

97 MATHEMATICS AND COMPUTING↗

Solid and gas thermal conductivity models improvement and validation in various porous insulation materials

In the past few decades, significant efforts have been made to improve the theoretical understanding of thermal transport mechanisms in thermal insulation materials and push the thermal conductivity's lower limits. However, most works focused singularly on specific types of materials, and the models used for thermal conductivity predictions are diverse - a model that fits one material might not fit others. Here, we improve and unify the gas and solid thermal conductivity models for porous materials. Through experimental characterization of several different materials as well as literature data for other materials, these models are validated. Further, we have also found that the pressure-dependent gas thermal conductivity of most materials can be well fitted by using one or two pore sizes without using a complex pore size distribution. With the refined models, we decompose the effective thermal conductivity of several thermal insulation materials into gas, solid, and radiation contributions. For cellular (polystyrene and polyurethane) foams, the relative contributions from air, solid, and radiation are 58–75%, 3–11%, 16–38%, respectively. For granular porous materials (polyurethane and silica in this work), the contributions from air, solid, and radiation are 45–66%, 34–46%, and 0–8%, respectively. This work is expected to provide guidance on the design and optimization of the next generation of thermal insulation materials, for example, through the effort of reducing gas conduction and radiation in foams and suppressing gas and solid conduction in aerogels.

36 MATERIALS SCIENCE↗

Overcoming significant challenges in extracting off-stoichiometric thermodynamics using the compound energy formalism through complementary use of experimental and first principles data: A case study of Ba 1-x Sr x FeO 3-δ

The compound energy formalism (CEF) is a powerful framework to describe the thermodynamics of metal oxides as a function of off-stoichiometry, temperature, and composition. The thermodynamic properties are crucial materials design attributes in metal oxide-based oxygen-exchange chemical processes. Despite the richness of information an accurate CEF model provides, a method to determine a unique and accurate fit for oxygen exchange materials remains elusive. This contribution details a method for fitting the CEF model that overcomes the current fitting challenges through three innovations: 1) the combination of density functional theory calculations with experimental data decorrelates excess terms and delineates the enthalpic/entropic contributions to the Gibbs free energy; 2) a systematic determination of the important CEF model terms, removing thermodynamic predetermining human intervention; 3) a self-consistent solution of the starting oxygen offstoichiometry (δ 0 ) of thermogravimetric measurements. Thus, our method enables the reliable extraction of off-stoichiometric metal oxide thermodynamic properties and facilitates rapid materials compositional screening, and reliable process design of systems dependent on off-stoichiometric redox-active metal oxides. We apply this method to a Ba x Sr 1-x FeO 3-δ test case. We find by systematically examining the performance of the CEF model fit with and without each innovation that all three innovations are necessary for an accurate fit. We determined that reduction enthalpy is higher and more sensitive to off-stoichiometry when the Sr fraction is large (139.5 and 185.3 kJ/mol O 2 for SrFeO 3 at δ = 0 and δ = 0.5, respectively vs. nearly constant 83 kJ/mol O 2 for BaFeO 3 ). However, the reduction entropy is mostly insensitive to Sr fraction, but highly dependent on δ suggesting larger contributions in the non-configurational entropy parameters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Timing the r -process Enrichment of the Ultra-faint Dwarf Galaxy Reticulum II

The ultra-faint dwarf galaxy Reticulum II (Ret II) exhibits a unique chemical evolution history, with ${72}_{-12}^{+10}$% of its stars strongly enhanced in r-process elements. We present deep Hubble Space Telescope photometry of Ret II and analyze its star formation history. As in other ultra-faint dwarfs, the color–magnitude diagram is best fit by a model consisting of two bursts of star formation. If we assume that the bursts were instantaneous, then the older burst occurred around the epoch of reionization, forming ~80% of the stars in the galaxy, while the remainder of the stars formed ~3 Gyr later. When the bursts are allowed to have nonzero durations, we obtain slightly better fits. The best-fitting model in this case consists of two bursts beginning before reionization, with approximately half the stars formed in a short (100 Myr) burst and the other half in a more extended period lasting 2.6 Gyr. Considering the full set of viable star formation history models, we find that 28% of the stars formed within 500 ± 200 Myr of the onset of star formation. The combination of the star formation history and the prevalence of r-process-enhanced stars demonstrates that the r-process elements in Ret II must have been synthesized early in its initial star-forming phase. We therefore constrain the delay time between the formation of the first stars in Ret II and the r-process nucleosynthesis to be less than 500 Myr. This measurement rules out an r-process source with a delay time of several Gyr or more, such as GW170817.

79 ASTRONOMY AND ASTROPHYSICS↗

Time-integrated activity coefficient estimation for radionuclide therapy using PET and a pharmacokinetic model: A simulation study on the effect of sampling schedule and noise

Purpose: The aim of this study was to investigate the accuracy of PET-based treatment planning for predicting the time-integrated activity coefficients (TIACs). Methods: The parameters of a physiologically based pharmacokinetic (PBPK) model were fitted to the biokinetic data of 15 patients to derive assumed true parameters and were used to construct true mathematical patient phantoms (MPPs). Biokinetics of 150 MBq {sup 68}Ga-DOTATATE-PET was simulated with different noise levels [fractional standard deviation (FSD) 10%, 1%, 0.1%, and 0.01%], and seven combinations of measurements at 30 min, 1 h, and 4 h p.i. PBPK model parameters were fitted to the simulated noisy PET data using population-based Bayesian parameters to construct predicted MPPs. Therapy simulations were performed as 30 min infusion of {sup 90}Y-DOTATATE of 3.3 GBq in both true and predicted MPPs. Prediction accuracy was then calculated as relative variability v{sub organ} between TIACs from both MPPs. Results: Large variability values of one time-point protocols [e.g., FSD = 1%, 240 min p.i., v{sub kidneys} = (9 ± 6)%, and v{sub tumor} = (27 ± 26)%] show inaccurate prediction. Accurate TIAC prediction of the kidneys was obtained for the case of two measurements (1 and 4 h p.i.), e.g., FSD = 1%, v{sub kidneys} = (7 ± 3)%, and v{sub tumor} = (22 ± 10)%, or three measurements, e.g., FSD = 1%, v{sub kidneys} = (7 ± 3)%, and v{sub tumor} = (22 ± 9)%. Conclusions: {sup 68}Ga-DOTATATE-PET measurements could possibly be used to predict the TIACs of {sup 90}Y-DOTATATE when using a PBPK model and population-based Bayesian parameters. The two time-point measurement at 1 and 4 h p.i. with a noise up to FSD = 1% allows an accurate prediction of the TIACs in kidneys.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Uncertain characterization of reservoir fluids due to brittleness of equation of state regression

Equations of state (EoS) play a central role in modeling the phase equilibrium of fluid mixtures. Their parameterization involves fitting a model to experimental data, i.e., solving a nonlinear, non-convex, multivariate optimization problem. The latter requires one to select design variables, domains of definition for each variable, and weights assigned to individual measurements. We demonstrate that subjective choices of an optimization algorithm and an initial guess also impact the regression process. Consequently, EoS predictions are fundamentally uncertain even after the EoS tuning to a limited set of experimental data points. We demonstrate this observation for two hydrocarbon reservoir fluids, in which five properties of the heaviest carbon fraction are treated as design variables. While all the optimization algorithms and initial guesses match experimental data for the gas and liquid properties, the resulting EoS parameterizations lead to dramatically different predictions of the fluid’s thermophysical behavior in the unsampled pressure and temperature regions. In conclusion, we propose the probabilistic treatment of design variables to quantify the predictive uncertainty of the resulting fluid models.

15 GEOTHERMAL ENERGY↗

Surrogate modelling the Baryonic Universe II: On forward modelling the colours of individual and populations of galaxies

ABSTRACT Among the properties shaping the light of a galaxy, the star formation history (SFH) is one of the most challenging to model due to the variety of correlated physical processes regulating star formation. In this work, we leverage the stellar population synthesis model fsps, together with SFHs predicted by the hydrodynamical simulation IllustrisTNG and the empirical model universemachine, to study the impact of star formation variability on galaxy colours. We start by introducing a model-independent metric to quantify the burstiness of a galaxy formation model, and we use this metric to demonstrate that universemachine predicts SFHs with more burstiness relative to IllustrisTNG. Using this metric and principal component analysis, we construct families of SFH models with adjustable variability, and we show that the precision of broad-band optical and near-infrared colours degrades as the level of unresolved short-term variability increases. We use the same technique to demonstrate that variability in metallicity and dust attenuation presents a practically negligible impact on colours relative to star formation variability. We additionally provide a model-independent fitting function capturing how the level of unresolved star formation variability translates into imprecision in predictions for galaxy colours; our fitting function can be used to determine the minimal SFH model that reproduces colours with some target precision. Finally, we show that modelling the colours of individual galaxies with per cent-level precision demands resorting to complex SFH models, while producing precise colours for galaxy populations can be achieved using models with just a few degrees of freedom.

Chaves-Montero, Jonás (ORCID:0000000295534261)↗

A nondestructive method to estimate the chlorophyll content of Arabidopsis seedlings

Chlorophyll content decreases in plants under stress conditions, therefore it is used commonly as an indicator of plant health. Arabidopsis thaliana offers a convenient and fast way to test physiological phenotypes of mutations and treatments. But, chlorophyll measurements with conventional solvent extraction are not applicable to Arabidopsis leaves due to their small size, especially when grown on culture dishes. We provide a nondestructive method for chlorophyll measurement whereby the red, green and blue (RGB) values of a color leaf image is used to estimate the chlorophyll content from Arabidopsis leaves. The method accommodates different profiles of digital cameras by incorporating the ColorChecker chart to make the digital negative profiles, to adjust the white balance, and to calibrate the exposure rate differences caused by the environment so that this method is applicable in any environment. We chose an exponential function model to estimate chlorophyll content from the RGB values, and fitted the model parameters with physical measurements of chlorophyll contents. As further proof of utility, this method was used to estimate chlorophyll content of G protein mutants grown on different sugar to nitrogen ratios. Our method is a simple, fast, inexpensive, and nondestructive estimation of chlorophyll content of Arabidopsis seedlings. This method lead to the discovery that G proteins are important in sensing the C/N balance to control chlorophyll content in Arabidopsis.

59 BASIC BIOLOGICAL SCIENCES↗

West Flank Coso FORGE: Magnetotelluric Inversion

The Coso Magnetotelluric (MT) dataset of which the West Flank FORGE MT data is a subset, was collected by Schlumberger / WesternGeco and initially processed by the WesternGeco GeoSolutions Integrated EM Center of Excellence in Milan, Italy. The 2011 data was based on 99 soundings that were centered on the West Flank geothermal prospect. The new soundings along with previous data from 2003 and 2006 were incorporated into a 3D inversion. Full impedance tensor data were inverted in the 1-3000 Hz range. The modelling report notes several noise sources, specifically the DC powerline that is 20,000 feet west of the survey area, and may have affected data in the 0.02 to 10 Hz range. Model cell dimensions of 450 x 450 x 65 feet were used to avoid computational instability in the 3D model. The fit between calculated and observed MT values for the final model run had an RMS value of 1.807. The included figure from the WesternGeco report shows the sounding locations from the 2011, 2006 and 2003 surveys.

15 GEOTHERMAL ENERGY↗

Probing Galaxy assembly bias in BOSS galaxies using void probabilities

ABSTRACT We measure the void probability function (VPF) of galaxies in the Baryon Oscillation Spectroscopic Survey (BOSS). The VPF provides complementary information to standard two-point statistics in that it is sensitive to galaxy bias in the most extreme underdensities in the cosmic web. Thus, the VPF is ideal for testing whether halo occupation of galaxies depends on large-scale density, an effect known as galaxy assembly bias. We find that standard HOD model – one parameterized by halo mass only – fit only to the two-point function, accurately predicts the VPF. Additionally, for HOD models where density dependence is explicitly incorporated, the best-fitting models fit to the combination of the correlation function and the VPF have zero density dependence. Thus, galaxy assembly bias is not a strong source of systematic uncertainty when modelling the clustering of massive galaxies.

Walsh, Kilian↗

Ion mass dependence of irradiation-induced damage accumulation in KTaO 3

Damage production and amorphization resulting from the interaction of medium-energy (from 40 to 480 keV) noble-gas ions (from He to Kr) with potassium tantalate (KTaO 3 ) are determined using ion channeling measurements. A disorder accumulation model has been fit to the maximum damage concentration versus ion fluence to extract the cross sections for direct-impact and defect-stimulated amorphization, and the results indicate that defect-stimulated amorphization is the dominant mechanism. These cross sections exhibit a strong dependence on the calculated cross sections for displacing lattice atoms, indicating a dominant contribution of nuclear interactions to the defect production and amorphization processes under the irradiation conditions used in this study. Furthermore, these experimental findings, along with the model fits, suggest that the difference in recoil spectra between He and the other heavier ions may be the main driving force for the decreased damage efficiency observed for He ions, which results in a reduced rate of damage accumulation.

36 MATERIALS SCIENCE↗

Formal Experimentation and Analysis of Handheld RFID Readers as a Tool for Nuclear Material Accounting

Commercial, off-the-shelf Radio Frequency Identification (RFID) systems have been successfully deployed for inventory tracking in numerous industries, but their viability in the tracking of complex environments containing nuclear material is less understood. Of primary interest in this setting is a RFID tracking system for nuclear material accounting which can reliably identify as many tags as possible with high accuracy, while also reducing an operator’s exposure to radiation. In this work, we develop a formal statistical approach to identify relevant handheld RFID reader settings which optimize tagging performance. To achieve this goal, we design a full factorial split-plot experiment for a static shelf configuration scene with 50 randomly placed nuclear material containers affixed with RFID tags. We use Bayesian inference to fit a second-order response surface model which expresses the probability of a successful match for each container as a function of the experimental factors. Such effects are allowed to vary across individual containers and the containers’ population in its entirety to estimate overall effects. Uncertainties of estimates and predictions are quantified via their corresponding posterior distributions. Following extensive model checking and validation, the fitted model is used to identify experimental factors which maximize matching probabilities at both the container-level and for the full shelf configuration scene. We also analyze sensitivity of performance to relevant factors.

25 ENERGY STORAGE↗

Exploring the spatial variability of soil properties in an Alfisol Catena

Detailed digital soil maps showing the spatial heterogeneity of soil properties consistent with the landscape are required for site-specific management of plant nutrients, land use planning and process-based environmental modeling. We characterized the short-scale spatial heterogeneity of soil properties in an Alfisol catena in a tropical landscape of Sri Lanka. The impact of different land-uses (paddy, vegetable and un-cultivated) was examined to assess the impact of anthropogenic activities on the variability of soil properties at the catenary level. Conditioned Latin hypercube sampling was used to collect 58 geo-referenced topsoil samples (0–30 cm) from the study area. Soil samples were analyzed for pH, electrical conductivity (EC), organic carbon (OC), cation exchange capacity (CEC) and texture. The spatial correlation between soil properties was analyzed by computing crossvariograms and subsequent fitting of theoretical model. Spatial distribution maps were developed using ordinary kriging. The range of soil properties, pH: 4.3–7.9; EC: 0.01–0.18 dS m –1 ; OC: 0.1–1.37%; CEC: 0.44– 11.51 cmol (+) kg –1 ; clay: 1.5–25% and sand: 59.1–84.4% and their coefficient of variations indicated a large variability in the study area. Electrical conductivity and pH showed a strong spatial correlation which was reflected by the cross-variogram close to the hull of the perfect correlation. Moreover, cross-variograms calculated for EC and Clay, CEC and OC, CEC and clay and CEC and pH indicated weak positive spatial correlation between these properties. Relative nugget effect (RNE) calculated from variograms showed strongly structured spatial variability for pH, EC and sand content (RNE < 25%) while CEC, organic carbon and clay content showed moderately structured spatial variability (25% < RNE < 75%). Spatial dependencies for examined soil properties ranged from 48 to 984 m. The mixed effects model fitting followed by Tukey's post-hoc test showed significant effect of land use on the spatial variability of EC. Our study revealed a structured variability of topsoil properties in the selected tropical Alfisol catena. Except for EC, observed variability was not modified by the land uses. Investigated soil properties showed distinct spatial structures at different scales and magnitudes of strength. Our results will be useful for digital soil mapping, site specific management of soil properties, developing appropriate land use plans and quantifying anthropogenic impacts on the soil system.

Geo-statistical analysis↗

Evaluating penalized logistic regression models to predict Heat-Related Electric grid stress days

Understanding the conditions associated with stress on the electricity grid is important in the development of contingency plans for maintaining reliability during periods when the grid is stressed. In this paper, heat-related grid stress and the relationship with weather conditions is examined using data from the eastern United States. Penalized logistic regression models were developed and applied to predict stress on the electric grid using weather data. The inclusion of other weather variables, such as precipitation, in addition to temperature improved model performance. Several candidate models and datasets were examined. A penalized logistic regression model fit at the operation-zone level was found to provide predictive value and interpretability. Additionally, the importance of different weather variables observed at different time scales were examined. Maximum temperature and precipitation were identified as important across all zones while the importance of other weather variables was zone specific. The methods presented in this work are extensible to other regions and can be used to aid in planning and development of the electrical grid.

Bramer, L. M.↗

WE-H-BRA-07: Mechanistic Modelling of the Relative Biological Effectiveness of Heavy Charged Particles

Purpose Uncertainty in the Relative Biological Effectiveness (RBE) of heavy charged particles compared to photons remains one of the major uncertainties in particle therapy. As RBEs depend strongly on clinical variables such as tissue type, dose, and radiation quality, more accurate individualised models are needed to fully optimise treatments. MethodsWe have developed a model of DNA damage and repair following X-ray irradiation in a number of settings, incorporating mechanistic descriptions of DNA repair pathways, geometric effects on DNA repair, cell cycle effects and cell death. Our model has previously been shown to accurately predict a range of biological endpoints including chromosome aberrations, mutations, and cell death. This model was combined with nanodosimetric models of individual ion tracks to calculate the additional probability of lethal damage forming within a single track. These lethal damage probabilities can be used to predict survival and RBE for cells irradiated with ions of different Linear Energy Transfer (LET). ResultsBy combining the X-ray response model with nanodosimetry information, predictions of RBE can be made without cell-line specific fitting. The model’s RBE predictions were found to agree well with empirical proton RBE models (Mean absolute difference between models of 1.9% and 1.8% for cells with α/β ratios of 9 and 1.4, respectively, for LETs between 0 and 15 keV/µm). The model also accurately recovers the impact of high-LET carbon ion exposures, showing both the reduced efficacy of ions at extremely high LET, as well as the impact of defects in non-homologous end joining on RBE values in Chinese Hamster Ovary cells.ConclusionOur model is predicts RBE without the inclusion of empirical LET fitting parameters for a range of experimental conditions. This approach has the potential to deliver improved personalisation of particle therapy, with future developments allowing for the calculation of individualised RBEs. SJM is supported by a Marie Curie International Outgoing Fellowship from the European Commission’s FP7 program (EC FP7 MC-IOF-623630)

60 APPLIED LIFE SCIENCES↗

Luminescence methodology to determine grain-boundary, grain-interior, and surface recombination in thin-film solar cells

We determine the grain-boundary (GB) recombination velocity, S GB , and grain-interior (GI) lifetime, t GI , parameters in superstrate CdS/CdTe thin-film solar cell technology by combining cathodoluminescence (CL) spectrum imaging and time-resolved photoluminescence (TRPL) measurements. We consider critical device formation stages, including after CdTe deposition, CdCl 2 treatment, and Cu diffusion. CL image analysis methods extract GB and GI intensities and grain size for hundreds of grains per sample. Concurrently, a three-dimensional CL model is developed to simulate the GI intensity as a function of t GI , S GB , grain size, and the surface recombination velocity, S surf . TRPL measurements provide an estimate of S surf for the CL model. A fit of GI intensity vs. grain size data with the CL model gives a self-consistent and representative set of S GB and t GI values for the samples: S GB (t GI ) = 2.6 x 10 6 cm/s (68-250 ps), S GB (t GI ) = 4.1 x 10 5 cm/s (1.5-3.3 ns), and S GB (t GI ) = 5.5 x 10 5 cm/s (1.0-3.8 ns) for as-deposited, CdCl 2 -treated, and CdCl 2- A nd Cu-treated samples, respectively. Thus, we find that the CdCl2 treatment both helps to passivate GBs and significantly increase the GI lifetime. Subsequent Cu diffusion increases GB recombination slightly and has nuanced effects on the GI lifetime. Finally, as a partial check on the S GB and t GI values, they are input to a Sentaurus device model, and the simulated performance is compared to the measured performance. The methodology developed here can be applied broadly to CdTe and CdSeTe thin-film technology and to other thin-film solar cell materials including Cu(In 1-x Ga x )Se 2 , Cu 2 ZnSnS 4 , and perovskites.

14 SOLAR ENERGY↗