Scale Testing Utilizing Additive Manufacturing and CCO Pallet [Poster]
This project created models for fit and function checks as well as to evaluate worker ergonomics. The CAD models were split into printable sections and assembled.
42 ENGINEERING↗
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.
Full Modeling and parameter compression methods in configuration space for DESI 2024 and beyond
In the contemporary era of high-precision spectroscopic surveys, led by projects like DESI, there is an increasing demand for optimizing the extraction of cosmological information from clustering data. This work conducts a thorough comparison of various methodologies for modeling the full shape of the two-point statistics in configuration space. We investigate the performance of both direct fits (Full Modeling) and the parameter compression approaches (ShapeFit and Standard). We utilize the ABACUS-SUMMIT simulations, tailored to exceed DESI's precision requirements. Particularly, we fit the two-point statistics of three distinct tracers (LRG, ELG, and QSO), by employing a Gaussian Streaming Model in tandem with Convolution Lagrangian Perturbation Theory and Effective Field Theory. We explore methodological setup variations, including the range of scales, the set of galaxy bias parameters, the inclusion of the hexadecapole, as well as model extensions encompassing varying ns and allowing for w 0 w a CDM dark energy model. Throughout these varied explorations, while precision levels fluctuate and certain configurations exhibit tighter parameter constraints, our pipeline consistently recovers the parameter values of the mocks within 1σ in all cases for a 1-year DESI volume. Additionally, we compare the performance of configuration space analysis with its Fourier space counterpart using three models: PyBird, FOLPS and velocileptors, presented in companion papers. We find good agreement with the results from all these models.
79 ASTRONOMY AND ASTROPHYSICS↗
A critical assessment of solutions to the galaxy diversity problem
Galactic rotation curves exhibit a diverse range of inner slopes. Observational data indicates that explaining this diversity may require a mechanism that correlates a galaxy's surface brightness with the central-most region of its dark matter halo. Here, in this work, we compare several concrete models that capture the relevant physics required to explain the galaxy diversity problem. We focus specifically on a Self-Interacting Dark Matter (SIDM) model with an isothermal core and two Cold Dark Matter (CDM) models with/without baryonic feedback. In contrast to the CDM case, the SIDM model can lead to the formation of an isothermal core in the halo, and is also mostly insensitive to baryonic feedback processes, which act on longer time-scales. Using rotation curves from 90 galaxies in the Spitzer Photometry & Accurate Rotation Curves (SPARC) catalog, we perform a comprehensive model comparison that addresses issues of statistical methodology from prior works. The best-fit halo models that we recover are consistent with standard CDM concentration-mass and abundance matching relations. We find that both the SIDM and feedback-affected CDM models are better than a CDM model with no feedback in explaining the rotation curves of low and high surface brightness galaxies in the sample. However, when compared to each other, there is no strong statistical preference for either the SIDM or the feedback-affected CDM halo model as the source of galaxy diversity in the SPARC catalog.
79 ASTRONOMY AND ASTROPHYSICS↗
Robustness of cosmic birefringence measurement against Galactic foreground emission and instrumental systematics
The polarization of the cosmic microwave background (CMB) can be used to search for parity-violating processes like that predicted by a Chern-Simons coupling to a light pseudoscalar field. Such an interaction rotates E modes into E modes in the observed CMB signal through an effect known as cosmic birefringence. Even though isotropic birefringence can be confused with the rotation produced by a miscalibration of the detectors' polarization angles, the degeneracy between both effects is broken when Galactic foreground emission is used as a calibrator. Here, in this work, we use realistic simulations of the High-Frequency Instrument of the Planck mission to test the impact that Galactic foreground emission and instrumental systematics have on the recent birefringence measurements obtained through this technique. Our results demonstrate the robustness of the methodology against the miscalibration of polarization angles and other systematic effects, like intensity-to-polarization leakage, beam leakage, or cross-polarization effects. However, our estimator is sensitive to the EB correlation of polarized foreground emission. Here we propose to correct the bias induced by dust EB by modeling the foreground signal with templates produced in Bayesian component-separation analyses that fit parametric models to CMB data. Acknowledging the limitations of currently available dust templates like that of the Commander sky model, high-precision CMB data and a characterization of dust beyond the modified blackbody paradigm are needed to obtain a definitive measurement of cosmic birefringence in the future.
79 ASTRONOMY AND ASTROPHYSICS↗
Reduced-Order Modeling of Hydrogen Releases from Vent Stacks and with Wind
This report describes computational modeling in the HyRAM+ software for study of hydrogen behavior in two common scenarios. First, models of unignited plumes exiting a vent stack were considered. It was shown that entrainment and vent backpressure were major factors in plume physics. Second, HyRAM+ was extended to include wind effects on plumes by modifying the plume momentum balance and entrainment modeling. Use of the model showed that plume shape and length changed with wind speed and direction; in all cases, wind causes a shortening of the plume along the streamline. While the no-wind case in HyRAM+ has been validated and the newly developed wind model was fitted to very limited experimental data, more controlled experimental configurations would help validate the models and ensure accurate simulation of hydrogen plume behavior for vent stack releases or in wind.
08 HYDROGEN↗
A small-box approach to the local crystal structure of Y 3 NbO 7
The local crystal structure of Y 3 NbO 7 was probed using neutron total scattering. Five different structural models were fit to the experimental pair distribution function of this material up to 15 Å. Fits of defect fluorite (Fm$\overline{3}$m), pyrochlore (Fd$\overline{3}$m), and orthorhombic C weberite models (C222 1 ) were attempted based on prior structural investigations of Y 3 NbO 7 and related fluorite-based materials. Orthorhombic P weberite (P222 1 ) and monoclinic weberite models (P112 1 ) were derived via symmetry reduction of the orthorhombic C weberite structure. All models except monoclinic weberite failed to reproduce the experimental pair distribution function. Misfits were clearly visible for peaks arising from atom pairs within the first coordination shell (defect fluorite and pyrochlore) and beyond (orthorhombic C and P weberite). By contrast, monoclinic weberite (γ ≈ 91.8°) provided an excellent fit, with no major misfits observed. This unit cell features distorted YO 8 dodecahedra, YO 7 monocapped octahedra, and NbO 6 octahedra as building blocks. Comparison between the topologies of the cationic and anionic substructures in monoclinic weberite and its parent orthorhombic C structure revealed displacements of the metal atoms along the [100] direction and a less compact spatial distribution of oxygen atoms along the a axis. In conclusion, the monoclinic weberite model presented herein may serve as a starting point to revisit structural modeling of fluorite-based materials with challenging local structures.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Tape-in-tape-out (TITO): a new approach for in-situ contact resistance measurements in high temperature superconducting CORC® cables
One of the ongoing development challenges with ReBCO high-temperature superconducting (HTS) cables is normal zone initiation and local heating, which is associated with over-critical current flowing through and around local performance reductions in individual tapes. Although inter-tape contact resistances are well-reported for individual tapes and HTS cables, these measurements are a challenge in CORC® cables as current percolates through much of the cable. In this work, developments in tape-in-tape-out (TITO) automated individual tape powering experiments are presented, and a simplified modeling approach for current percolation in CORC® cables is developed. An optimization is formulated to fit the model parameters to a large set of TITO experiments on a single cable, allowing the layer-dependent inter-tape contact resistance to be extracted. Measurements are presented and discussed for a straight CORC® cable and a cable bent to a 152 mm and 76 mm radius. The approach provides both quantitative insight and facilitates better understanding of current sharing phenomena in HTS cables.
Teyber, Reed↗
A Novel Multiphysics Multiscale Multiporosity Shale Gas Transport Model for Geomechanics/Flow Coupling in Steady and Transient States
Summary A novel multiphysics multiscale multiporosity shale gas transport (M3ST) model was developed to investigate shale gas transport in both transient and steady states. The microscale model component contains a kerogen domain and an inorganic matrix domain, and each domain has its own geomechanical and gas transport properties. Permeabilities of various shale cores were measured in the laboratory using a pulse decay permeameter (PDP) with different pore pressure and confining stress combinations. The PDP-measured apparent permeability as a function of pore pressure under two effective stresses was fitted using the microscale M3ST model component based on nonlinear least squares fitting (NLSF), and the fitted model parameters were able to provide accurate model predictions for another effective stress. The parameters and petrophysical properties determined in the steady state were then used in the transient-state, continuum-scale M3ST model component, which performed history matching of the evolutions of the upstream and downstream gas pressures. In addition, a double-exponential empirical model was developed as a powerful alternative to the M3ST model to fit laboratory-measured apparent permeability under various effective stresses and pore pressures. The developed M3ST model and the research findings in this study provided critical insights into the role of the multiphysics mechanisms, including geomechanics, fluid dynamics and transport, and the Klinkenberg effect on shale gas transport across different spatial scales in both steady and transient states.
Engineering↗
Reconstructing the Universe: Testing the Mutual Consistency of the Pantheon and SDSS/eBOSS BAO Data Sets with Gaussian Processes
We test the mutual consistency between the baryon acoustic oscillation measurements from the eBOSS SDSS final release and the Pantheon supernova compilation in a model-independent fashion using Gaussian process regression. We also test their joint consistency with the ΛCDM model in a model-independent fashion. We also use Gaussian process regression to reconstruct the expansion history that is preferred by these two data sets. While this methodology finds no significant preference for model flexibility beyond ΛCDM, we are able to generate a number of reconstructed expansion histories that fit the data better than the best-fit ΛCDM model. These example expansion histories may point the way toward modifications to ΛCDM. We also constrain the parameters Ω{sub k} and H {sub 0} r {sub d} both with ΛCDM and with Gaussian process regression. We find that H {sub 0} r {sub d} = 10,030 ± 130 km s{sup −1} and Ω{sub k} = 0.05 ± 0.10 for ΛCDM and that H {sub 0} r {sub d} = 10,040 ± 140 km s{sup −1} and Ω{sub k} = 0.02 ± 0.20 for the Gaussian process case.
71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
Gemini Planet Imager Spectroscopy of the Dusty Substellar Companion HD 206893 B
We present new near-infrared Gemini Planet Imager (GPI) spectroscopy of HD 206893 B, a substellar companion orbiting within the debris disk of its F5V star. The J, H, K1, and K2 spectra from GPI demonstrate the extraordinarily red colors of the object, confirming it as the reddest substellar object observed to date. The significant flux increase throughout the infrared presents a challenging atmosphere to model with existing grids. Best-fit values vary from 1200 to 1800 K for effective temperature and from 3.0 to 5.0 for log(g), depending on which individual wavelength band is fit and which model suite is applied. The extreme redness of the companion can be partially reconciled by invoking a high-altitude layer of submicron dust particles, similar to dereddening approaches applied to the peculiar red field L dwarf population. However, reconciling the HD 206893 B spectra with even those of the reddest low-gravity L dwarf spectra still requires the contribution of additional atmospheric dust, potentially due to the debris disk environment in which the companion resides. Orbit fitting from 4 yr of astrometric monitoring is consistent with a ∼30 yr period, an orbital inclination of 147°, and a semimajor axis of 10 au, well within the estimated disk inner radius of ∼50 au. As one of a very few substellar companions imaged interior to a circumstellar disk, the properties of this system offer important dynamical constraints on companion–disk interaction and provide a benchmark for substellar and planetary atmospheric study.
47 OTHER INSTRUMENTATION↗
Kepler and TESS Observations of PG 1159-035
Abstract PG 1159-035 is the prototype of the PG 1159 hot (pre-)white dwarf pulsators. This important object was observed during the Kepler satellite K2 mission for 69 days in 59 s cadence mode and by the TESS satellite for 25 days in 20 s cadence mode. We present a detailed asteroseismic analysis of those data. We identify a total of 107 frequencies representing 32 ℓ = 1 modes, 27 frequencies representing 12 ℓ = 2 modes, and eight combination frequencies. The combination frequencies and the modes with very high k values represent new detections. The multiplet structure reveals an average splitting of 4.0 ± 0.4 μ Hz for ℓ = 1 and 6.8 ± 0.2 μ Hz for ℓ = 2, indicating a rotation period of 1.4 ± 0.1 days in the region of period formation. In the Fourier transform of the light curve, we find a significant peak at 8.904 ± 0.003 μ Hz suggesting a surface rotation period of 1.299 ± 0.002 days. We also present evidence that the observed periods change on timescales shorter than those predicted by current evolutionary models. Our asteroseismic analysis finds an average period spacing for ℓ = 1 of 21.28 ± 0.02 s. The ℓ = 2 modes have a mean spacing of 12.97 ± 0.4 s. We performed a detailed asteroseismic fit by comparing the observed periods with those of evolutionary models. The best-fit model has T eff = 129, 600 ± 11 100 K, M * = 0.565 ± 0.024 M ⊙ , and log g = 7.41 − 0.54 + 0.38 , within the uncertainties of the spectroscopic determinations. We argue for future improvements in the current models, e.g., on the overshooting in the He-burning stage, as the best-fit model does not predict excitation for all of the pulsations detected in PG 1159-035.
79 ASTRONOMY AND ASTROPHYSICS↗
CF + excitation in the interstellar medium
The detection of CF + in interstellar clouds potentially allows astronomers to infer the elemental fluorine abundance and the ionization fraction in ultraviolet-illuminated molecular gas. Because local thermodynamic equilibrium (LTE) conditions are hardly fulfilled in the interstellar medium (ISM), the accurate determination of the CF + abundance requires one to model its non-LTE excitation via both radiative and collisional processes. Here, we report quantum calculations of rate coefficients for the rotational excitation of CF + in collisions with para- and ortho-H 2 (for temperatures up to 150 K). As an application, we present non-LTE excitation models that reveal population inversion in physical conditions typical of ISM photodissociation regions (PDRs). We successfully applied these models to fit the CF + emission lines previously observed toward the Orion Bar and Horsehead PDRs. The radiative transfer models achieved with these new rate coefficients allow the use of CF + as a powerful probe to study molecular clouds exposed to strong stellar radiation fields.
79 ASTRONOMY AND ASTROPHYSICS↗
The kinetics of SARS-CoV-2 infection based on a human challenge study
Studying the early events that occur after viral infection in humans is difficult unless one intentionally infects volunteers in a human challenge study. Here, we use data about severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in such a study in combination with mathematical modeling to gain insights into the relationship between the amount of virus in the upper respiratory tract and the immune response it generates. We propose a set of dynamic models of increasing complexity to dissect the roles of target cell limitation, innate immunity, and adaptive immunity in determining the observed viral kinetics. We introduce an approach for modeling the effect of humoral immunity that describes a decline in infectious virus after immune activation. We fit our models to viral load and infectious titer data from all the untreated infected participants in the study simultaneously. We found that a power-law with a power h < 1 describes the relationship between infectious virus and viral load. Viral replication at the early stage of infection is rapid, with a doubling time of ~2 h for viral RNA and ~3 h for infectious virus. We estimate that adaptive immunity is initiated ~7 to 10 d postinfection and appears to contribute to a multiphasic viral decline experienced by some participants; the viral rebound experienced by other participants is consistent with a decline in the interferon response. Altogether, we quantified the kinetics of SARS-CoV-2 infection, shedding light on the early dynamics of the virus and the potential role of innate and adaptive immunity in promoting viral decline during infection.
59 BASIC BIOLOGICAL SCIENCES↗
Battery Performance, Thermal, and Life Modeling for Southern California Edison (Cooperative Research and Development Final Report)
The objective of this project is to develop battery performance models, thermal models and life models for SCE to use in simulation and analysis of grid services involving energy storage. NREL will fit the models to three separate chemistries using test data provided by SCE. The dataset for each chemistry will include electrical, thermal and aging response to different temperatures and cycling conditions, measured under a variety of cell-level and module-level experiments described in SCE test protocol documents. Once tuned to a test dataset, the battery system model software will provide predictions of battery energy and power loss, and thus cycle and calendar life, for any energy storage grid service power profile of interest to SCE. The automated battery life modeling and simulation tool can also be applied to any other grid battery systems as long as necessary input data is available to accelerate battery model development and battery lifetime analysis.
25 ENERGY STORAGE↗
Combining Experimental and Theoretical Techniques to Gain an Atomic Level Understanding of the Defect Binding Mechanism in Hard Carbon Anodes for Sodium Ion Batteries
Sodium ion batteries (NIBs) are an attractive alternative to lithium–ion batteries in applications that require large–scale energy storage due to sodium's high natural abundance and low cost. Hard carbon (HC) is the most promising anode material for NIBs; however, there is a knowledge gap in the understanding of the sodium binding mechanism that prevents a rational design of HC. This study tunes sucrose–derived HC via synthesis temperature then evaluates the structural, physical, and electrochemical properties. Neutron total scattering is used to generate structural models by fitting pair distribution functions (PDF) with a combination of molecular dynamics and reverse Monte Carlo methods. From this model, the number and type of structural features are identified, quantified, and correlated to the galvanostatic charge/discharge. A method of PDF “fingerprinting” binding sites using Na probe atoms is developed and analyzing these PDFs reveals an atomistic view of ion binding sites responsible for “defect” storage mechanisms. Combining these techniques results in an atomic–level study that provides a big picture of the Na–binding mechanism in NIBs, which allows for more precise tuning of the structure–property relationships in the future. Finally, the methodologies developed will also enable new strategies for the analysis of amorphous functional materials.
25 ENERGY STORAGE↗
Forward modeling approach to nuclear reaction cross sections: Applications in neutron inelastic scattering
The development of nuclear reaction models for the production of evaluated nuclear data has traditionally been performed by comparing measured cross sections with predictions from reaction model codes whose physical input parameters are adjusted to obtain the best agreement between measured and modeled results. To more directly probe reaction model inputs, this work introduces a forward modeling approach to experimental reaction cross-section determination, where the most important physical input parameters to reaction model calculations are obtained via 𝜒 2 minimization between measured and calculated observables. This was demonstrated using data collected by the Gamma Energy Neutron Energy Spectrometer for Inelastic Scattering (GENESIS) at the 88-inch cyclotron at Lawrence Berkeley National Laboratory, a detection array consisting of organic liquid scintillators and high-purity germanium (HPGe) detectors. Using a broad-spectrum neutron beam and a 99.98%-enriched 56 Fe target, GENESIS was used to perform a simultaneous measurement of 56 Fe 𝛾-ray production cross sections and secondary neutron energy and angle distributions. The results of the forward modeling approach to the determination of energy-differential 𝛾-ray production cross sections for the yrast 4 + → 2 + and 6 + → 4 + transitions, as well as eight other off-yrast transitions, were compared against those obtained using conventional techniques, and the results are in good agreement. In addition to discrete 𝛾-ray yield total scattered neutron energy-angular distributions as a function of incident neutron energy were also obtained using forward modeling and found to agree with evaluated data, with the exception of elastic scattering at small angles. The fitted reaction model parameters obtained through forward modeling were also used to calculate the cross section for the unobserved (𝑛, 2𝑛) reaction; excellent agreement with the current evaluation was obtained, providing a validation of the predictive capabilities of the forward model approach. This work bridges the gap between nuclear data experiment and evaluation by providing a new means for extracting inelastic neutron-scattering cross sections and neutron-induced 𝛾-ray production data while directly probing reaction model physics.
73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
Exploring the Whole Set of Accurate Sparse Interpretable Models
In data science applications, there are often many models that fit the data well. This phenomenon was called the Rashomon Effect by Leo Breiman. The set of good models is called the Rashomon Set, and the goal of this project is to locate, store, and study the Rashomon sets for classes of interpretable models, including decision trees and generalized additive models.
97 MATHEMATICS AND COMPUTING↗
Red Noise–based False Alarm Thresholds for Astrophysical Periodograms via Whittle’s Approximation to the Likelihood
Astronomers who search for periodic signals using Lomb–Scargle periodograms rely on false alarm level (FAL) estimates to identify statistically significant peaks. Although FALs are often calculated from white noise models, many astronomical time series suffer from red noise. Prewhitening is a statistical technique in which a continuum model is subtracted from the log power spectrum estimate, after which the observer can proceed with a white-noise treatment. Here we present a prewhitening-based method of calculating frequency-dependent FALs. We fit power laws and autoregressive models of order 1 to each Lomb–Scargle periodogram by minimizing the Whittle approximation to the negative log-likelihood (NLL), then calculate FALs based on the best-fit model power spectrum. Our technique is a novel extension of the Whittle NLL to datasets with uneven time sampling. We demonstrate FAL calculations using observations of α Cen B, GJ 581, HD 192310, synthetic data from the radial velocity (RV) fitting challenge, and Kepler observations of a differential rotator. The Kepler data analysis shows that only true rotation signals are detected by red noise FALs, while white noise FALs suggest all spurious peaks in the low-frequency range are significant. A high-frequency sinusoid injected into α Cen B logR$'$ HK observations exceeds the 1% red noise FAL despite having only 8.9% of the power of the dominant rotation signal. In a periodogram of HD 192310 RVs, peaks associated with differential rotation and planets are detected against the 5% red noise FAL without iterative model fitting or subtraction. The software for calculating red noise–based FALs is available on GitHub.
Astrostatistics (1882)↗