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 361 records · Page 20

Assessment of Response Surface Models using Independent Confirmation Point Analysis

This paper highlights various advantages that confirmation-point residuals have over conventional model design-point residuals in assessing the adequacy of a response surface model fitted by regression techniques to a sample of experimental data. Particular advantages are highlighted for the case of design matrices that may be ill-conditioned for a given sample of data. The impact of both aleatory and epistemological uncertainty in response model adequacy assessments is considered.

DeLoach, Richard↗

The SAGA Survey. V. Modeling Satellite Systems around Milky Way–Mass Galaxies with Updated UniverseMachine

Abstract Environment plays a critical role in shaping the assembly of low-mass galaxies. Here, we use the U niverse M achine (UM) galaxy–halo connection framework and Data Release 3 of the Satellites Around Galactic Analogs (SAGA) Survey to place dwarf galaxy star formation and quenching into a cosmological context. UM is a data-driven forward model that flexibly parameterizes galaxy star formation rates (SFRs) using only halo mass and assembly history. We add a new quenching model to UM, tailored for galaxies with m ⋆ ≲ 10 9 M ⊙ , and constrain the model down to m ⋆ ≳ 10 7 M ⊙ using new SAGA observations of 101 satellite systems around Milky Way (MW)–mass hosts and a sample of isolated field galaxies in a similar mass range from the Sloan Digital Sky Survey. The new best-fit model, “UM-SAGA,” reproduces the satellite stellar mass functions, average SFRs, and quenched fractions in SAGA satellites while keeping isolated dwarfs mostly star-forming. The enhanced quenching in satellites relative to isolated field galaxies leads the model to maximally rely on halo assembly to explain the observed environmental quenching. Extrapolating the model down to m ⋆ ∼ 10 6.5 M ⊙ yields a quenched fraction of ≳30% for isolated field galaxies and ≳80% for satellites of MW-mass hosts at this stellar mass. Spectroscopic surveys can soon test this specific prediction to reveal the relative importance of internal feedback, cessation of mass and gas accretion, satellite-specific gas processes, and reionization for the evolution of faint low-mass galaxies.

Wang, Yunchong (ORCID:000000018913626X)↗

Disc dichotomy signature in the vertical distribution of [Mg/Fe] and the delayed gas infall scenario

Context. Analysis of the Apache Point Observatory Galactic Evolution Experiment project (APOGEE) data suggests the existence of a clear distinction between two sequences of disc stars in the [α/Fe] versus [Fe/H] abundance ratio space, known as the high- and low-α sequence, respectively. This dichotomy also emerges from an analysis of the vertical distribution of the [α/Fe] abundance ratio. Aims. We aim to test whether the revised two-infall chemical evolution models designed to reproduce the low- and high-α sequences in the [α/Fe] versus [Fe/H] ratios in the solar neighbourhood are also capable of predicting the disc bimodality observed in the vertical distribution of [Mg/Fe] in APOGEE DR16 data. Methods. Along with the chemical composition of the simple stellar populations born at different Galactic times predicted by our reference chemical evolution models in the solar vicinity, we provide their maximum vertical height above the Galactic plane |zmax| computed assuming the relation between the vertical action and stellar age in APOGEE thin-disc stars. Result. The vertical distribution of the [Mg/Fe] abundance ratio predicted by the reference chemical evolution models is in agreement with that observed when combining the APOGEE DR16 data (chemical abundances) with the astroNN catalogue (stellar ages, orbital parameters) for stars younger than 8 Gyr (only low-α sequence stars). Including the high-α disc component, the dichotomy in the vertical [Mg/Fe] abundance distribution is reproduced considering the observational cut in the Galactic height of |z|< 2 kpc. However, our model predicts an overly flat (almost constant) growth of the maximum vertical height |zmax| quantity as a function of [Mg/Fe] for high-α objects in contrast with the median values from APOGEE data. Possible explanations for such a tension are that: (i) the APOGEE sample with |z|< 2 kpc is more likely than ours to be contaminated by halo stars, causing the median values to be kinematically hotter, and (ii) external perturbations – such as minor mergers – that the Milky Way experienced in the past could have heated up the disc, and the heating of the orbits cannot be modeled by only scattering processes. Assuming a disc dissection based on chemistry for APOGEE-DR16 stars (|z|< 2 kpc), the observed |zmax| distributions for high-α and low-α sequences are in good agreement with our model predictions if we consider the errors in the vertical action estimates in the calculation. Moreover, a better agreement between predicted and observed stellar distributions at different Galactic vertical heights is achieved if asteroseismic ages are included as a constraint in the best-fit model calculations. Conclusions. The signature of a delayed gas infall episode, which gives rise to a hiatus in the star formation history of the Galaxy, are imprinted both in the [Mg/Fe] versus [Fe/H] relation and in vertical distribution of [Mg/Fe] abundances in the solar vicinity.

79 ASTRONOMY AND ASTROPHYSICS↗

Origin of gaplike behaviors in URu 2 Si 2 : Combined study via quasiparticle scattering spectroscopy and resistivity measurements

In this work, we address two long-standing questions regarding the hidden order in URu 2 Si 2 : Is it associated with the hybridization process, and what are the distinct roles played by the localized and itinerant electrons? Our quasiparticle scattering spectroscopy reveals a hybridization gap ubiquitous in the entire phase space spanned by P and Fe substitutions in URu 2 Si 2 , including the no-order and antiferromagnetic regions, with minimal change upon crossing the phase boundary. This indicates its opening is not associated with the ordering, and thus localized electrons must be the major player. Towards a consistent understanding of all the other gaplike behaviors observed only below transition temperatures, we analyze the electrical resistivity using a model in which gapped bosonic excitations are the dominant scattering source. With their stiffness set to follow an unusual temperature dependence (decreasing with decreasing temperature), this model fits all of our resistivity data well including the jump at the transition. Remarkably, the extracted gap increases slowly with increasing Fe content, similarly to the gap detected by inelastic neutron scattering at $Q_1$ = (1.4, 0, 0), suggesting a common origin. Such a model can also naturally explain the Hall effect temperature dependence without invoking Fermi surface gapping.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Kepler Data Validation I: Architecture, Diagnostic Tests, and Data Products for Vetting Transiting Planet Candidates

The Kepler Mission was designed to identify and characterize transiting planets in the Kepler Field of View and to determine their occurrence rates. Emphasis was placed on identification of Earth-size planets orbiting in the Habitable Zone of their host stars. Science data were acquired for a period of four years. Long-cadence data with 29.4 min sampling were obtained for approx. 200,000 individual stellar targets in at least one observing quarter in the primary Kepler Mission. Light curves for target stars are extracted in the Kepler Science Data Processing Pipeline, and are searched for transiting planet signatures. A Threshold Crossing Event is generated in the transit search for targets where the transit detection threshold is exceeded and transit consistency checks are satisfied. These targets are subjected to further scrutiny in the Data Validation (DV) component of the Pipeline. Transiting planet candidates are characterized in DV, and light curves are searched for additional planets after transit signatures are modeled and removed. A suite of diagnostic tests is performed on all candidates to aid in discrimination between genuine transiting planets and instrumental or astrophysical false positives. Data products are generated per target and planet candidate to document and display transiting planet model fit and diagnostic test results. These products are exported to the Exoplanet Archive at the NASA Exoplanet Science Institute, and are available to the community. We describe the DV architecture and diagnostic tests, and provide a brief overview of the data products. Transiting planet modeling and the search for multiple planets on individual targets are described in a companion paper. The final revision of the Kepler Pipeline code base is available to the general public through GitHub. The Kepler Pipeline has also been modified to support the Transiting Exoplanet Survey Satellite (TESS) Mission which is expected to commence in 2018.

data analysis↗

Abundance matching analysis of the emission-line galaxy sample in the extended Baryon Oscillation Spectroscopic Survey

We present the measurements of the small-scale clustering for the emission-line galaxy (ELG) sample from the extended Baryon Oscillation Spectroscopic Survey (eBOSS) in the Sloan Digital Sky Survey IV (SDSS-IV). We use conditional abundance matching method to interpret the clustering measurements from 0.34 to 70 h -1 Mpc. In order to account for the correlation between properties of ELGs and their environment, we add a secondary connection between star formation rate of ELGs and halo accretion rate. Three parameters are introduced to model the ELG [O II ] luminosity and to mimic the target selection of eBOSS ELGs. The parameters in our models are optimized using Markov Chain Monte Carlo (MCMC) method. We find that by conditionally matching star formation rate of galaxies and the halo accretion rate, we are able to reproduce the eBOSS ELG small-scale clustering within 1σ error level. Our best-fitting model shows that the eBOSS ELG sample only consists of ~ 12 per cent of all star-forming galaxies, and the satellite fraction of eBOSS ELG sample is 19.3 per cent. We show that the effect of assembly bias is ~ 20 per cent on the two-point correlation function and ~ 5 per cent on the void probability function at scale of r ~ 20 h -1 Mpc.

79 ASTRONOMY AND ASTROPHYSICS↗

Carbon-Binder Weight Loading Optimization for Improved Lithium-Ion Battery Rate Capability

Battery performance is strongly correlated with electrode microstructure and weight loading of the electrode components. Among them are the carbon-black and binder additives that enhance effective conductivity and provide mechanical integrity. However, these both reduce effective ionic transport in the electrolyte phase and reduce energy density. Therefore, an optimal additive loading is required to maximize performance, especially for fast charging where ionic transport is essential. Such optimization analysis is however challenging due to the nanoscale imaging limitations that prevent characterizing this additive phase and thus quantifying its impact on performance. Herein, an additive-phase generation algorithm has been developed to remedy this limitation and identify percolation threshold used to define a minimal additive loading. Improved ionic transport coefficients from reducing additive loading has been then quantified through homogenization calculation, macroscale model fitting, and experimental symmetric cell measurement, with good agreement between the methods. Rate capability test demonstrates capacity improvement at fast charge at the beginning of life, from 37% to 55%, respectively for high and low additive loading during 6C CC charging, in agreement with macroscale model, and attributed to a combination of lower cathode impedance, reduced electrode tortuosity and cathode thickness.

25 ENERGY STORAGE↗

Carbon-Binder Optimization for Lithium-Ion Battery Extreme Fast Charge

Battery performance is strongly correlated with electrode microstructure and weight loading of the electrode components. Among them are the carbon-black and binder additives that enhance effective conductivity and provide mechanical integrity. However, these both reduce effective ionic transport in the electrolyte phase and reduce energy density. Therefore, an optimal additive loading is required to maximize performance, especially for fast charging where ionic transport is essential. Such optimization analysis is however challenging due to the nanoscale imaging limitations that prevent characterizing this additive phase and thus quantifying its impact on performance. Herein, an additive-phase generation algorithm has been developed to remedy this limitation and identify percolation threshold used to define a minimal additive loading. Improved ionic transport coefficients from reducing additive loading has been then quantified through homogenization calculation, macroscale model fitting, and experimental symmetric cell measurement, with good agreement between the methods. Rate capability test demonstrates capacity improvement at fast charge at the beginning of life, from 37% to 55%, respectively for high and low additive loading during 6C CC charging, in agreement with macroscale model, and attributed to a combination of lower cathode impedance, reduced electrode tortuosity and cathode thickness.

carbon-binder additives↗

Modeling the structure of H I in the galactic disk

For comparison with a detailed H I survey, we calculated models for the H I intensity distribution, T(B)(l, b, v), for a 18 x 12 deg region in the disk of the Galaxy toward l = 205 deg. In one class of models, neutral hydrogen is in the form of clouds, while in the other class of models it is a continuous medium with holes whose sizes are typical of supernova remnants and superbubbles. The cloud models are a poor representation of the data in that the amplitude fluctuations are dissimilar and the continuity of the H I observations is not reproduced. A continuous H I medium with holes provides a much better match to the data, and the best-fit models suggest that the filling factor for the H I is 20-50 percent.

Bregman, Joel N.↗

Predicting fatigue from heart rate signatures using functional logistic regression

Physical fatigue can have adverse effects on humans in extreme environments. Therefore, being able to predict fatigue using easy to measure metrics such as heart rate (HR) signatures has potential to have an impact in real-life scenarios. We apply a functional logistic regression model that uses HR signatures to predict physical fatigue, where physical fatigue is defined in a data-driven manner. Data were collected using commercially available wearable devices on 47 participants hiking the 20.7-mile Grand Canyon rim-to-rim trail in a single day. Fitted model provides good predictions and interpretable parameters for real-life application.

60 APPLIED LIFE SCIENCES↗

Heavy flavor dynamics across system size at the LHC

One of the fundamental signatures of the Quark Gluon Plasma has been the suppression of heavy flavor (specifically D mesons), which has been measured via the nuclear modification factor, R AA and azimuthal anisotropies, v n , in large systems. However, multiple competing models can reproduce the same data for R AA to v n . In this talk we break down the competing effects that conspire together to successfully reproduce R AA and v n in experimental data using Trento+v-USPhydro+DAB-MOD. Then using our best fit model we make predictions for R AA and v n across system size for 208 PbPb, 129 XeXe, 40 ArAr, and 16 OO collisions. We find that 0–10% centrality has a non-trivial interplay between the system size and eccentricities such that system size effects are masked in v 2 whereas in 30–50% centrality the eccentricities are approximately constant across system size and, therefore, is a better centrality class to study D meson dynamics across system size.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An accurate and interpretable model for antimicrobial resistance in pathogenic Escherichia coli from livestock and companion animal species

Understanding the microbial genomic contributors to antimicrobial resistance (AMR) is essential for early detection of emerging AMR infections, a pressing global health threat in human and veterinary medicine. Here we used whole genome sequencing and antibiotic susceptibility test data from 980 disease causing Escherichia coli isolated from companion and farm animals to model AMR genotypes and phenotypes for 24 antibiotics. We determined the strength of genotype-to-phenotype relationships for 197 AMR genes with elastic net logistic regression. Model predictors were designed to evaluate different potential modes of AMR genotype translation into resistance phenotypes. Our results show a model that considers the presence of individual AMR genes and total number of AMR genes present from a set of genes known to confer resistance was able to accurately predict isolate resistance on average (mean F 1 score = 98.0%, SD = 2.3%, mean accuracy = 98.2%, SD = 2.7%). However, fitted models sometimes varied for antibiotics in the same class and for the same antibiotic across animal hosts, suggesting heterogeneity in the genetic determinants of AMR resistance. We conclude that an interpretable AMR prediction model can be used to accurately predict resistance phenotypes across multiple host species and reveal testable hypotheses about how the mechanism of resistance may vary across antibiotics within the same class and across animal hosts for the same antibiotic.

Chung, Henri C.↗

The Solar Chromosphere/Corona Interface: FUV-EUV Observations and Modeling of Unresolved Coronal Funnels - I

A coronal funnel model, developed by Rabin (199l), was tested against a calibrated spectroheliogram recorded in 171 - 175 Angstrom bandpass. This image was recorded on board a sounding rocket experiment flown on 1994 November 3, called the Multi-Spectral Solar Telescope Array, II (MSSTA II), MSSTA, a joint project of Stanford University, the NASA Marshall Space Flight Center and the Lawrence Livermore National Laboratory, is an observing platform composed of a set of normal-incidence, multilayer-coated optics designed to obtain narrow bandpass, high resolution images (1 - 3 arc sec) at selected FUV, EUV and soft X-ray wavelengths (44 Angstroms - 1550 Angstroms). Using full-disk images centered at 1550 Angstroms (C IV) and 173 Angstroms (FE IX/X), the funnel model, which is based on coronal back-heating, was tested against the data incorporating observed constraints on global coverage and measured flux. Found, was a class of funnel models that could account for the quiescent, globally diffuse and unresolved emission seen in the 171 - 175 Angstrom bandpass, where the funnels are assumed to be rooted in the C IV supergranular network. These models, when incorporated with the Chianti spectral code, suggest that this emission is mostly of upper transition region origin and primarily composed of FE IX plasma. The funnels are found to have constrictions, Gamma approx. 6 - 20, which is in good agreement with the observations. Further, the fitted models simultaneously satisfy global areal constraints seen in both images; namely, that a global network of funnels must cover approx. 70 - 95 % of the total solar surface area seen in the 171 - 175 Angstrom image, and = 45 % of the disk area seen in the 1550 Angstrom bandpass. These findings support the configuration of the EUV magnetic network as suggested by Reeves et af. (1974) and put forth in more detail by Gabriel (1976). Furthermore, the models are in good agreement with differential emission measure estimates made of the transition region by Raymond & Doyle (1981) for temperatures, 250,000 K = T = 650,000 K, based on full-disk observations made on board by SkyLab.

Martinez-Galarce, D. S.↗

The Solar Chromosphere/Corona Interface. I. Far-Ultraviolet to Extreme-Ultraviolet Observations and Modeling of Unresolved Coronal Funnels

A coronal funnel model, developed by D. Rabin, was tested against a calibrated spectroheliogram recorded in the 170-1 75 A bandpass. This image was recorded on board a sounding-rocket experiment flown on 1994 November 3, called the Multi-Spectral Solar Telescope Array II (MSSTA II). MSSTA, a joint project of Stanford University, the NASA Marshall Space Flight Center, and the Lawrence Livermore National Laboratory' is an observing platform composed of a set of normal-incidence, multilayer-coated optics designed to obtain narrow-bandpass, high-resolution images (1 sec.- 3 sec.) at selected far-ultraviolet (FUV), extreme-ultraviolet (EUV), and soft X-ray wavelengths (44-1550 A). Using full disk images centered at 1550 A (C IV) and 173 A (Fe IX/X), the funnel model, which is based on coronal back-heating, was tested against the data incorporating observed constraints on global coverage and measured flux. Found was a class of funnel models that could account for the quiescent, globally diffuse and unresolved emission seen in the 171-175 A bandpass, where the funnels are assumed to be rooted in the C IV supergranular network. These models, when incorporated with the CHIANTI spectral code, suggest that this emission is mostly of upper transition region origin and primarily composed of Fe IX plasma. The funnels are found to have constrictions, Gamma approx. 6-20, which is in good agreement with the observations. Further, the fitted models simultaneously satisfy global areal constraints seen in both images; namely,that a global network of funnels must cover approx. 700 - 95% of the total solar surface area seen in the 171-175 A image, and a 5% of the disk area seen in the 1550 A bandpass. These findings support the configuration of the EUV magnetic network as suggested by Reeves et al. and put forth in more detail by Gabriel. Furthermore, the models are in good agreement with differential emission measure estimates made of the transition region by J. C. Raymond & J. G. Doyle for temperatures 250,000 K less than or = T less than or = 650,000 K, based on full-disk observations made on board Skylab.

Martinez-Galarce, Dennis S.↗

The Near-Earth Encounter of Asteroid 308635 (2005 YU55): Thermal IR Observations

The near-Earth approach (0.00217 AU, or 0.845 lunar distances) of the C-type asteroid 308635 (2005 YU55) in November 2011 presented a rare opportunity for detailed observations of a low-albedo NEA in this size range. As part of a multi-telescope campaign to measure visible and infrared spectra and photometry, we obtained mid-infrared (approx. 8 to 22 micron) photometry and spectroscopy of 2005 YU55 using Michelle on the Gemini North telescope on UT November 9 and 10,2011. An extensive radar campaign together with optical light-curves established the rotation state of YU55. In addition, the radar imaging resulted in a shape model for the asteroid, detection of numerous boulders on its surface, and a preliminary estimate of its equatorial diameter at 380 +/- 20 m. In a preliminary analysis, applying the radar and lightcurve-derived parameters to a rough-surface thermophysical model fit to the Gemini/Michelle thermal emission photometry results in a thermal inertia range of approximately 500 to 1500 J/sq m/0.5s/K, with the low-thermal-inertia solution corresponding to the small end of the radar size range and vice versa. Updates to these results will be presented and modeling of the thermal contribution to the measured near-infrared spectra from Palomar/Triplespec and IRTF/SpeX will also be discussed.

Lim, Lucy F.↗

Architectures of Exoplanetary Systems: Towards a Multi-planet Model for Reproducing the Planet Radius Valley and Intra-system Size Similarity

Multi-planet systems provide crucial insights into the architectures and correlations within planetary systems, which in turn offer clues into their formation and evolution histories. Through the thousands of exoplanet candidates that it discovered, NASA's Kepler mission revealed that systems with multiple transiting planets are common. We describe a new model for inferring the exoplanet population observed by Kepler, which combines a semi-parametric description of planetary systems with a prescription for envelope mass-loss due to photoevaporation. This new model is a"hybrid" of the maximum AMD (angular momentum deficit) model in which the orbital excitations of the planets are drawn by distributing the system's critical AMD (He et al. 2020), and a photoevaporation model for sculpting the observed radius valley (Neil & Rogers 2020). This is an update to the SysSim suite of exoplanet population models derived from forward modeling the Kepler detection pipeline. The hybrid model is constrained by the observed statistics of multi-transiting systems, including additional constraints from the distribution of planet radii and their intra-system size patterns. This model is capable of producing a planet radius valley as well as a preference for planet size similarity. However, we find that the Kepler catalog exhibits features that are even more significant (i.e., it has a deeper observed radius valley and stronger size similarity) than the best-fit model. We use this model to infer the intrinsic planet radius valley as well as the distribution of primordial planet masses and core masses.

kepler↗

Measurement of the effective leptonic weak mixing angle

Using pp collision data at $\sqrt{s}$ = 13 TeV, recorded by the LHCb experiment between 2016 and 2018 and corresponding to an integrated luminosity of 5.4 fb −1 , the forward-backward asymmetry in the pp → Z/γ* → μ + μ − process is measured. The measurement is carried out in ten intervals of the difference between the muon pseudorapidities, within a fiducial region covering dimuon masses between 66 and 116 GeV, muon pseudorapidities between 2.0 and 4.5 and muon transverse momenta above 20 GeV. These forward-backward asymmetries are compared with predictions, at next-to-leading order in the strong and electroweak couplings. The measured effective leptonic weak mixing angle is ${\sin}^2{\theta}_{\textrm{eff}}^{\ell }=0.23147\pm 0.00044\pm 0.00005\pm 0.00023,$ where the first uncertainty is statistical, the second arises from systematic uncertainties associated with the asymmetry measurement, and the third arises from uncertainties in the fit model used to extract ${\sin}^2{\theta}_{\textrm{eff}}^{\ell }$ from the asymmetry measurement. This result is based on an arithmetic average of results using the CT18, MSHT20, and NNPDF31 parameterisations of the proton internal structure, and is consistent with previous measurements and with predictions from the global electroweak fit.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Enhanced mechanical property evaluation using innovative data analytics capability

This report focuses on efforts to improve Multiphysics Object-Oriented Simulation Environ-ment (MOOSE) for mechanical property evaluation using data analytics. These efforts involveimprovements to the stochastic tools module (STM) for stochastic simulations of MOOSE mul-tiphysics model and the development of inverse optimization capabilities. The report gives anoverview of the STM and describes recent updates to its core capabilities and theory on itsreduced-order model (ROM) schemes. Examples are also provided showing the impact of theseupdates and exhibits the usefulness of ROMs. An overview of the gradient based inverse opti-mziation algorithms are given along with examples of their application to source identification.Inverse optimization will provide a new methodology in STM for fitting model parameters toexperimental data.

97 MATHEMATICS AND COMPUTING↗