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 739 records · Page 41
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Broadband Spectral Investigations of Magnetar Bursts

We present our broadband (2–250 keV) time-averaged spectral analysis of 388 bursts from SGR J1550−5418, SGR 1900+14, and SGR 1806−20 detected with the Rossi X-ray Timing Explorer ( RXTE ) here and as a database in a companion web-catalog. We find that two blackbody functions (BB+BB), the sum of two modified blackbody functions (LB+LB), the sum of a blackbody function and a power-law function (BB+PO), and a power law with a high-energy exponential cutoff (COMPT) all provide acceptable fits at similar levels. We performed numerical simulations to constrain the best fitting model for each burst spectrum and found that 67.6% of burst spectra with well-constrained parameters are better described by the Comptonized model. We also found that 64.7% of these burst spectra are better described with the LB+LB model, which is employed in the spectral analysis of a soft gamma repeater (SGR) for the first time here, than with the BB+BB and BB+PO models. We found a significant positive lower bound trend on photon index, suggesting a decreasing upper bound on hardness, with respect to total flux and fluence. We compare this result with bursts observed from SGR and AXP (anomalous X-ray pulsar) sources and suggest that the relationship is a distinctive characteristic between the two. We confirm a significant anticorrelation between burst emission area and blackbody temperature, and find that it varies between the hot and cool blackbody temperatures differently than previously discussed. We expand on the interpretation of our results in the framework of a strongly magnetized neutron star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Photometric Redshifts and Stellar Masses for Galaxies from the DESI Legacy Imaging Surveys

We present a catalog of photometric redshifts and stellar masses for galaxies from imaging surveys serving the Dark Energy Spectroscopic Instrument project. The imaging footprints cover an area of over 14,000 deg 2 . This catalog contains about 0.17 billion morphologically classified galaxies with r < 23 mag. It covers the redshift range of z < 1 and the mass range of $8.4\lt \mathrm{log}({M}_{* })\lt 11.9$. A total of five photometric bands are used, including three optical filters (g, r, and z) and two infrared filters ( Wide-field Infrared Survey Explorer W 1 and W 2). We use k-nearest neighbors and local linear regression algorithms to estimate photometric redshifts. Based on a spectroscopic sample of about 2.2 million galaxies, we achieve a redshift bias of $\overline{{\rm{\Delta }}{z}_{\mathrm{norm}}}$ = 2.4 × 10 -4 , accuracy of ${\sigma }_{{\rm{\Delta }}{z}_{\mathrm{norm}}}$ = 0.017, and outlier rate of about 5.1%. The stellar mass and other physical parameters are derived by stellar population synthesis model fitting based on our photometric redshifts. The differences in $\mathrm{log}({M}_{* })$ derived using spectroscopic and photometric redshifts present a 1σ dispersion of 0.09 dex.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Magellan/M2FS and MMT/Hectochelle Spectroscopy of Dwarf Galaxies and Faint Star Clusters within the Galactic Halo*

We present spectroscopic data for 16,369 stellar targets within and/or toward 38 dwarf spheroidal galaxies and faint star clusters within the Milky Way halo environment. All spectra come from observations with the multiobject, fiber-fed echelle spectrographs M2FS at the Magellan/Clay telescope or Hectochelle at the MMT, reaching a typical limiting magnitude G ≲ 21. Data products include processed spectra from all observations and catalogs listing estimates—derived from template model fitting—of line-of-sight velocity (median uncertainty 1.4 km s –1 ) effective temperature (255 K), (base-10 logarithm of) surface gravity (0.59 dex in cgs units), [Fe/H] (0.4 dex) and [Mg/Fe] (0.27 dex) abundance ratios. The sample contains multiepoch measurements for 3720 sources, with up to 15 epochs per source, enabling studies of intrinsic spectroscopic variability. The sample contains 6087 likely red giant stars (based on surface gravity), and 4492 likely members (based on line-of-sight velocity and Gaia-measured proper motion) of the target systems. The number of member stars per individual target system ranges from a few, for the faintest systems, to ~850 for the most luminous. For most systems, our new samples extend over wider fields than have previously been observed; of the likely members in our samples, 820 lie beyond 2 times the projected half-light radius of their host system, and 42 lie beyond 5 R half .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Production and Acceleration of Antinuclei in Supernova Shockwaves

We compute the energy spectra of antideuterons ( d-bar ) and antihelium ( (He)-bar ) in cosmic rays (CRs) in a scenario where hadronic interactions inside supernova remnants (SNRs) can produce a diffusively shock-accelerated “source component” of secondary antinuclei. The key parameters that specify the SNR environment and the interstellar CR transport are tightly constrained with the new measurements provided by the AMS experiment on the B/C ratio and on the p-bar /p ratio. The best-fit models obtained from the two ratios are found to be inconsistent with each other, as the p-bar /p data require enhanced secondary production. Thus, we derive conservative (i.e., B/C-driven) and speculative ( p-bar /p-driven) upper limits to the SNR flux contributions for the d-bar and (He)-bar spectra in CRs, along with their standard secondary component expected from CR collisions in the interstellar gas. We find that the source component of antinuclei can be appreciable at kinetic energies above a few ∼10 GeV n{sup −1}, but it is always sub-dominant below a few GeV n{sup −1}, that is the energy window where dark matter (DM) annihilation signatures are expected to exceed the level of secondary production. We also find that the total (standard + SNR) flux of secondary d-bar and (He)-bar is tightly constrained by the data. Thus, the presence of interaction processes in SNRs does not critically affect the total background for DM searches.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Halo Pressure Profile through the Skew Cross-power Spectrum of the Sunyaev–Zel’dovich Effect and CMB Lensing in Planck

We measure the cosmic microwave background (CMB) skewness power spectrum in Planck , using frequency maps of the HFI instrument and the Sunyaev–Zel’dovich (SZ) component map. The two-to-one skewness power spectrum measures the cross-correlation between CMB lensing and the thermal SZ effect. We also directly measure the same cross-correlation using the Planck CMB lensing map and the SZ map and compare it to the cross-correlation derived from the skewness power spectrum. We model fit the SZ power spectrum and CMB lensing–SZ cross-power spectrum via the skewness power spectrum to constrain the gas pressure profile of dark matter halos. The gas pressure profile is compared to existing measurements in the literature including a direct estimate based on the stacking of SZ clusters in Planck .

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

PSR J0030+0451 Mass and Radius from NICER Data and Implications for the Properties of Neutron Star Matter

Neutron stars are not only of astrophysical interest, but are also of great interest to nuclear physicists because their attributes can be used to determine the properties of the dense matter in their cores. One of the most informative approaches for determining the equation of state (EoS) of this dense matter is to measure both a star's equatorial circumferential radius R e and its gravitational mass M. Here we report estimates of the mass and radius of the isolated 205.53 Hz millisecond pulsar PSR J0030+0451 obtained using a Bayesian inference approach to analyze its energy-dependent thermal X-ray waveform, which was observed using the Neutron Star Interior Composition Explorer (NICER). This approach is thought to be less subject to systematic errors than other approaches for estimating neutron star radii. We explored a variety of emission patterns on the stellar surface. In this work, our best-fit model has three oval, uniform-temperature emitting spots and provides an excellent description of the pulse waveform observed using NICER. The radius and mass estimates given by this model are ${R}_{e}={13.02}_{-1.06}^{+1.24}$ km and $M={1.44}_{-0.14}^{+0.15}\,{M}_{\odot }$ (68%). The independent analysis reported in the companion paper by Riley et al. explores different emitting spot models, but finds spot shapes and locations and estimates of R e and M that are consistent with those found in this work. We show that our measurements of R e and M for PSR J0030+0451 improve the astrophysical constraints on the EoS of cold, catalyzed matter above nuclear saturation density.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Measurement of the Drell-Yan forward-backward asymmetry and of the effective leptonic weak mixing angle in proton-proton collisions at $\sqrt{s}$ = 13 TeV

The forward-backward asymmetry in Drell-Yan production and the effective leptonic electroweak mixing angle are measured in proton-proton collisions at $\sqrt{s}$ = 13 TeV, collected by the CMS experiment and corresponding to an integrated luminosity of 138 fb$^{-1}$. The measurement uses both dimuon and dielectron events, and is performed as a function of the dilepton mass and rapidity. The unfolded angular coefficient $A_4$ is also extracted, as a function of the dilepton mass and rapidity. Using the CT18Z set of parton distribution functions, we obtain $\sinθ^\ell_\text{eff}$ = 0.23152 $\pm$ 0.00031, where the uncertainty includes the experimental and theoretical contributions. The measured value agrees with the standard model fit result to global experimental data. This is the most precise $\sinθ^\ell_\text{eff}$ measurement at a hadron collider, with a precision comparable to the results obtained at LEP and SLD.

FOS: Physical sciences↗
DOE OSTI2022

Geographically aware estimates of remotely sensed water properties for Chesapeake Bay

This software contains the code for a data-driven model to produce high-resolution (< 1 km2) estimates of temperature, salinity, and turbidity over decadal time scales as required by next-generation Earth Systems Models. Our model fits daily MODIS Aqua reflectance data to surface observations (< 1 m depth) from 2000-2021 in Chesapeake Bay, USA.

Stachelek, Jemma↗
DOE OSTI2022

Towards 2+1 Flavor Lattice QCD Results for the Heavy Quark Diffusion Coefficient

We apply and extend a novel approach to non-perturbatively estimate the heavy-quark momentum diffusion coefficient κ, which is a key input for the theoretical description of heavy quarkonium production in heavy ion collisions, and is important for the understanding of the elliptic flow and nuclear suppression factor of heavy flavor hadrons. In the heavy-quark limit, this coefficient is encoded in the spectral functions of color-electric and color-magnetic correlators that we calculate on the lattice to high pre cision by applying gradient flow. In a recent study we have considered quenched QCD at 1.5 T c , where we performed a detailed study of the lat tice spacing and flow time dependence of the color-electric correlator, and, using theoretically well-established model fits for the spectral reconstruc tion, we estimated the heavy-quark diffusion coefficient. Equipped with the experience obtained in quenched QCD, we estimate $κ$ from 2+1 flavor QCD ensembles at small but finite lattice spacing and flow time without increasing systematic errors significantly.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2012

The effect of high-level waste glass composition on spinel liquidus temperature

Spinel crystals precipitate in high-level waste glasses containing Fe, Cr, Ni, Mn, Zn, and Ru. The liquidus temperature (T{sub L}d) of spinel as the primary crystallization phase is a function of glass composition, and the spinel solubility (c{sub o}) is a function of both glass composition and temperature (T). Previously reported models of T{sub L} as a function of composition are based on T{sub L} measured directly, which requires laborious experimental procedures. Viewing the curve of c{sub o} versus T as the liquidus line allows a significant broadening of the composition region for model fitting. This paper estimates T{sub L} as a function of composition based on c{sub o} data obtained with the X-ray diffraction technique.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Nanoscale structure in AgSbTe2 determined by diffuse elastic neutron scattering

Diffuse elastic neutron scattering measurements confirm that AgSbTe2 has a hierarchical structure, with defects on length scales from nanometers to microns. While scattering from mesoscale structure is consistent with previously-proposed structures in which Ag and Sb order on a NaCl lattice, more diffuse scattering from nanoscale structure suggests a structural rearrangement in which hexagonal layers form a combination of (ABC), (ABA), and (AAB) stacking sequences. The AgCrSe2 structure is the best-fitting model for the local atomic arrangements.

Specht, Eliot D↗
DOE OSTI2013

graphMC: A package for testing the independence of graphs

graphMC consists of an implementation of a set of statistical tests that ensure that a Markov chain in the space of graphs has converged to its stationary distribution. The tests are applicable when the graphs have labelled nodes. The tests are applied to the binary time-series traced out by an edge between two selected nodes. The software fits models of a first-order Markov process and an independent draws from a binary distribution to the time-series and selects one based on the BIC scores of the two models. The package consists of (1) an implementation of a Markov chain driving a graph rewiring process that preserves the joint degree distribution of the graph, (2) implementation of the statistical test mentioned above and (3) two examples that show how the Markov chain infrastructure and the statistical tests can be used.

Source record↗
DOE OSTI2017

Flow regimes during immiscible displacement

Fractional ow of immiscible phases occurs at the pore scale where grain surfaces and phases interfaces obstruct phase mobility. However, the larger scale behavior is described by a saturation-dependent phenomenological relationship called relative permeability. As a consequence, pore-scale parameters, such as phase topology and/ or geometry, and details of the flow regime cannot be directly related to Darcy-scale flow parameters. It is well understood that relative permeability is not a unique relationship of wetting-phase saturation and rather depends on the experimental conditions at which it is measured. Herein we use fast X-ray microcomputed tomography to image pore-scale phase arrangements during fractional flow and then forward simulate the flow regimes using the lattice-Boltzmann method to better understand the underlying pore-scale flow regimes and their influence on Darcy-scale parameters. We find that relative permeability is highly dependent on capillary number and that the Corey model fits the observed trends. At the pore scale, while phase topologies are continuously changing on the scale of individual pores, the Euler characteristic of the nonwetting phase (NWP) averaged over a sufficiently large field of view can describe the bulk topological characteristics; the Euler characteristic decreases with increasing capillary number resulting in an increase in relative permeability. Lastly, we quantify the fraction of NWP that flows through disconnected ganglion dynamics and demonstrate that this can be a significant fraction of the NWP flux for intermediate wetting-phase saturation. Furthermore, rate dependencies occur in our homogenous sample (without capillary end effect) and the underlying cause is attributed to ganglion flow that can significantly influence phase topology during the fractional flow of immiscible phases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

On-the-fly Doppler Broadening with Probability Table Interpolation for Unresolved Resonance Region in RMC

With the rapid development of computational power and parallel algorithms, Monte Carlo method has been widely investigated and used in neutron transport simulations due to its advantages in geometry modelling and usage of continuous energy point-wise cross sections. Nowadays the road has been paved to couple Monte Carlo codes with thermal-hydraulics codes to make multi-physics calculations. One of the most important aspects of thermal-hydraulics feedback is the detailed temperature distribution, which demands temperature dependent cross sections in neutron transport simulations. The conventional pre-generated cross section libraries at each temperature point are no longer feasible for detailed temperature distributions because of tremendous memory footprint. Hence, different on-the-fly Doppler broadening methods have been proposed for resolved resonance energy region and thermal energy region, such as the interpolation, regress model fitting, multipole representation (MPR), target motion sampling (TMS) and different integration algorithms based on Doppler broadening equation. Compared to resolved resonance energy region and thermal energy region, the representation of cross sections in unresolved resonance region is slightly different. As a result of the limitation of experiment measurements, the overlapped resonances cannot be distinguished from each other and the statistic average resonance parameters are given in the evaluated nuclear data file. In early applications, infinite-dilute pointwise cross sections were obtained from average parameters and were directly used in neutron simulations. However, it was proven that these infinite-dilute pointwise cross sections cannot capture the energy self-shielding effect, thus leading to an overestimate of reaction rates and usually an underestimate of k{sub eff} owing to the increased absorption rates and increased scattering probability to energy below unresolved resonance region where the capture-to-fission ratio is typically higher. This phenomenon is of great importance especially for problems where there is a considerable flux in unresolved resonance region. Then the probability table method was put forward by Levitt to take energy self-shielding effect into account, which provides probability distribution function of cross sections in unresolved resonance region instead of specific cross section values. The probability distribution function is presented in the form of tables and the mean values are equal to the infinite dilution cross sections. To obtain the temperature dependent cross sections in unresolved energy region, there are two different methods. One is to calculate cross sections directly from average resonance parameters and their statistical distributions , which is effective and precise but time-consuming. The other is to interpolate between pre-generated probability table files, which is much more efficient. In this paper, on-the-fly Doppler broadening of cross sections in unresolved resonance regions based on probability tables interpolation was implemented in RMC codes as a supplement to on-the-fly Doppler broadening for resolved resonance energy region and thermal region, which is based on the integration of Doppler broadening equation. The results of criticality calculation and depletion calculation show that on-the-fly Doppler broadening for the whole energy region is accurate and efficient. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Constraining GREA, an alternative theory accounting for the present cosmic acceleration

The origin of the Universe's late-time accelerated expansion remains unknown. The General Relativistic Entropic Acceleration (GREA) theory offers a compelling alternative to $Λ$CDM, attributing cosmic acceleration to entropy growth associated with cosmic and black hole horizons, without invoking a cosmological constant. We test GREA against the latest DESI DR2 Baryon Acoustic Oscillations (BAO), multiple Type Ia supernova compilations (Union3, Pantheon$\texttt{+}$, DES-SN5YR), and cosmic microwave background (CMB) distance measurements. While GREA is not nested within $Λ$CDM, it achieves a comparable goodness-of-fit, highlighting its potential as a theoretically motivated framework that circumvents some of the fine-tuning issues of the standard $Λ$CDM cosmology. We find that the best-fit model features a transient phantom crossing at $z \lesssim 2$, with $w_a\equiv \mathrm{d} w(a=1)/\mathrm{d}a \simeq-0.3$, in good agreement with observations. However, its present-day value $w_0\equiv w(z=0)$ is tightly constrained to be $w_0\simeq-1$. Upcoming low-redshift (i.e. $z < 1$) cosmological probes, from both background and perturbations, will offer promising avenues for further exploring the viability of the GREA theory.

Calderon, R. [Prague, Inst. Phys.] (ORCID:00000002↗
DOE OSTI2016

A method for simultaneous linear optics and coupling correction for storage rings with turn-by-turn beam position monitor data

We propose a method to simultaneously correct linear optics errors and linear coupling for storage rings using turn-by-turn (TbT) beam position monitor (BPM) data. The independent component analysis (ICA) method is used to isolate the betatron normal modes from the measured TbT BPM data. The betatron amplitudes and phase advances of the projections of the normal modes on the horizontal and vertical planes are then extracted, which, combined with dispersion measurement, are used to fit the lattice model. The fitting results are used for lattice correction. The method has been successfully demonstrated on the NSLS-II storage ring.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

A Redox-Electrodialysis Model with Zero Fitting Parameters: Insights into Process Limitations, Design, and Material Interventions

Redox-Electrodialysis (r-ED) is an electrochemical desalination cell architecture that has recently received considerable interest, due to its low energy demand relative to electrochemical desalination technologies that rely on electrode-based ion removal. To further improve the energy efficiency of r-ED, we developed a lumped mathematical model with no adjustable parameters to investigate the various sources of overpotential within the cell. Existing models of electrodialysis and r-ED cells either do not accurately incorporate all phenomena contributing to the overpotential or utilize empirical fitting parameters. The model developed here indicates that ohmic overpotentials, especially in the diluate chamber, are the most significant contributors to energy losses. Based on this insight, we hypothesized that adding an ion exchange resin wafer in the diluate compartment would increase the ionic conductivity and decrease the energy demand. Experimental results showed an 18% reduction in specific energy use while achieving the same degree of salt removal (20 mM to 12 mM). Furthermore, the resin wafer enabled complete desalination to potable drinking levels at a current density previously unachievable within practical operating voltage limits (4.93 mA cm -2 ). We also expanded the model to explore differences in r-ED energy use between configurations using multiple cells and a single cell with increased area.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2014

MODELING THE NONLINEAR CLUSTERING IN MODIFIED GRAVITY MODELS. I. A FITTING FORMULA FOR THE MATTER POWER SPECTRUM OF f(R) GRAVITY

Based on a suite of N-body simulations of the Hu-Sawicki model of f(R) gravity with different sets of model and cosmological parameters, we develop a new fitting formula with a numeric code, MGHalofit, to calculate the nonlinear matter power spectrum P(k) for the Hu-Sawicki model. We compare the MGHalofit predictions at various redshifts (z ≤ 1) to the f(R) simulations and find that the relative error of the MGHalofit fitting formula of P(k) is no larger than 6% at k ≤ 1 h Mpc{sup –1} and 12% at k in (1, 10] h Mpc{sup –1}, respectively. Based on a sensitivity study of an ongoing and a future spectroscopic survey, we estimate the detectability of a signal of modified gravity described by the Hu-Sawicki model using the power spectrum up to quasi-nonlinear scales.

79 ASTRONOMY AND ASTROPHYSICS↗