Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Systematic error”

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 1,495 records · Page 83

Quantification of modeling uncertainty in the Rayleigh damping model

Understanding and accurately characterizing energy dissipation mechanisms in civil structures during earthquakes is an important element of seismic assessment and design. The most commonly used model is attributed to Rayleigh. This paper proposes a systematic approach to quantify the uncertainty associated with Rayleigh's damping model. Bayesian calibration with embedded model error is employed to treat the coefficients of the Rayleigh model as random variables using modal damping ratios. Through a numerical example, we illustrate how this approach works and how the calibrated model can address modeling uncertainty associated with the Rayleigh damping model.

42 ENGINEERING↗

Comparative Analysis of Imaging and Measurements of Micrometer-Scale Fracture Aperture Fields Within a Heterogeneous Rock Using PET and X-ray CT

Knowledge of the spatial distribution of fracture apertures is essential for reliable characterization of flow and transport processes in fractured systems and for better understanding of physicochemical matrix–fracture interactions. Here, we propose and test two image-based methods, thereby extending the current experimental capabilities to characterize aperture size distribution in structurally heterogeneous geologic porous media noninvasively. The first approach utilizes an inversion method based on the dataset acquired from positron emission tomography (PET) and the second approach considers an extension of the classic missing attenuation technique that relies on clinical X-ray computed tomography (X-ray CT). Independent sets of imaging experiments are conducted on a fractured basalt core with heterogeneous matrix properties and aperture distributions to compare the two methodologies. A repeat of each experiment is conducted to verify the proposed workflows. The performance of these two imaging techniques is systematically evaluated through the analysis of signal-to-noise ratio, minimum fracture size detectability, and measurement errors. While both approaches provide a reliable estimation of fracture aperture distributions, PET yields a signal-to-noise ratio that is substantially higher than the corresponding X-ray CT measurements. Furthermore, uncertainties of the aperture values for PET are considerably lower ($\bar {\sigma}_{\text {d}} = 15\%$) compared to those obtained from X-ray CT ($\bar {\sigma}_{\text {d}} = 29\%$), allowing for the detection of minimum aperture sizes of 20 $\mu$m with 70% confidence level. Finally, these approaches provide key experimental tools for better understanding dynamic hydromechanical fracture properties in geologic systems.

58 GEOSCIENCES↗

CORONAL MAGNETOGRAPHY OF A SIMULATED SOLAR ACTIVE REGION FROM MICROWAVE IMAGING SPECTROPOLARIMETRY

We have simulated the Expanded Owens Valley Solar Array (EOVSA) radio images generated at multiple frequencies from a model solar active region, embedded in a realistic solar disk model, and explored the resulting data cube for different spectral analysis schemes to evaluate the potential for realizing one of EOVSA’s most important scientific goals—coronal magnetography. In this paper, we focus on modeling the gyroresonance and free–free emission from an on-disk solar active region model with realistic complexities in electron density, temperature and magnetic field distribution. We compare the magnetic field parameters extrapolated from the image data cube along each line of sight after folding through the EOVSA instrumental profile with the original (unfolded) parameters used in the model. We find that even the most easily automated, image-based analysis approach (Level-0) provides reasonable quantitative results, although they are affected by systematic effects due to finite sampling in the Fourier (UV) plane. Finally, we note the potential for errors due to misidentified harmonics of the gyrofrequency, and discuss the prospects for applying a more sophisticated spectrally based analysis scheme (Level-1) to resolve the issue in cases where improved UV coverage and spatial resolution are available.

79 ASTRONOMY AND ASTROPHYSICS↗

Using angular two-point correlations to self-calibrate the photometric redshift distributions of DECaLS DR9

Calibrating the redshift distributions of photometric galaxy samples is essential in weak lensing studies. The self-calibration method combines angular auto- and cross-correlations between galaxies in multiple photometric redshift (photo-z) bins to reconstruct the scattering rates matrix between redshift bins. In this paper, we test a recently proposed self-calibration algorithm using the DECaLS Data Release 9 and investigate to what extent the scattering rates are determined. We first mitigate the spurious angular correlations due to imaging systematics by a machine learning based method. We then improve the algorithm for χ 2 minimization and error estimation. Finally, we solve for the scattering matrices, carry out a series of consistency tests, and find reasonable agreements: (1) finer photo-z bins return a high-resolution scattering matrix, and it is broadly consistent with the low-resolution matrix from wider bins; (2) the scattering matrix from the Northern Galactic Cap is almost identical to that from Southern Galactic Cap; (3) the scattering matrices are in reasonable agreement with those constructed from the power spectrum and the weighted spectroscopic subsample. We also evaluate the impact of cosmic magnification. Although it changes little the diagonal elements of the scattering matrix, it affects the off-diagonals significantly. The scattering matrix also shows some dependence on scale cut of input correlations, which may be related to a known numerical degeneracy between certain scattering pairs. In conclusion, this work demonstrates the feasibility of the self-calibration method in real data and provides a practical alternative to calibrate the redshift distributions of photometric samples.

79 ASTRONOMY AND ASTROPHYSICS↗

Direct Measurement of the Nickel Spectrum in Cosmic Rays in the Energy Range from 8.8 GeV/$\mathcal{n}$ to 240 GeV/$\mathcal{n}$ with CALET on the International Space Station

The relative abundance of cosmic ray nickel nuclei with respect to iron is by far larger than for all other transiron elements; therefore it provides a favorable opportunity for a low background measurement of its spectrum. Since nickel, as well as iron, is one of the most stable nuclei, the nickel energy spectrum and its relative abundance with respect to iron provide important information to estimate the abundances at the cosmic ray source and to model the Galactic propagation of heavy nuclei. However, only a few direct measurements of cosmic-ray nickel at energy larger than ~ 3 GeV/n are available at present in the literature, and they are affected by strong limitations in both energy reach and statistics. In this Letter, we present a measurement of the differential energy spectrum of nickel in the energy range from 8.8 to 240 GeV/n , carried out with unprecedented precision by the Calorimetric Electron Telescope (CALET) in operation on the International Space Station since 2015. The CALET instrument can identify individual nuclear species via a measurement of their electric charge with a dynamic range extending far beyond iron (up to atomic number Z = 40 ). The particle’s energy is measured by a homogeneous calorimeter (1.2 proton interaction lengths, 27 radiation lengths) preceded by a thin imaging section (3 radiation lengths) providing tracking and energy sampling. This Letter follows our previous measurement of the iron spectrum [O. Adriani et al. (CALET Collaboration), Phys. Rev. Lett. 126, 241101 (2021).], and it extends our investigation on the energy dependence of the spectral index of heavy elements. It reports the analysis of nickel data collected from November 2015 to May 2021 and a detailed assessment of the systematic uncertainties. In the region from 20 to 240 GeV/n our present data are compatible within the errors with a single power law with spectral index -2.51 ± 0.07 .

79 ASTRONOMY AND ASTROPHYSICS↗

Single Pass Collider Memo: Gradient Perturbations of the SLC arc

As the beam passes through the arcs, the gradient it encounters at each magnet differs from the design value. This deviation may be in part random and in part systematic. In this note we make estimates of the effects to be expected from both kinds of errors.

43 PARTICLE ACCELERATORS↗

A 2.4% DETERMINATION OF THE LOCAL VALUE OF THE HUBBLE CONSTANT

We use the Wide Field Camera 3 (WFC3) on the Hubble Space Telescope (HST) to reduce the uncertainty in the local value of the Hubble constant from 3.3% to 2.4%. The bulk of this improvement comes from new near-infrared (NIR) observations of Cepheid variables in 11 host galaxies of recent type Ia supernovae (SNe Ia), more than doubling the sample of reliable SNe Ia having a Cepheid-calibrated distance to a total of 19; these in turn leverage the magnitude-redshift relation based on ∼300 SNe Ia at z < 0.15. All 19 hosts as well as the megamaser system NGC 4258 have been observed with WFC3 in the optical and NIR, thus nullifying cross-instrument zeropoint errors in the relative distance estimates from Cepheids. Other noteworthy improvements include a 33% reduction in the systematic uncertainty in the maser distance to NGC 4258, a larger sample of Cepheids in the Large Magellanic Cloud (LMC), a more robust distance to the LMC based on late-type detached eclipsing binaries (DEBs), HST observations of Cepheids in M31, and new HST -based trigonometric parallaxes for Milky Way (MW) Cepheids. We consider four geometric distance calibrations of Cepheids: (i) megamasers in NGC 4258, (ii) 8 DEBs in the LMC, (iii) 15 MW Cepheids with parallaxes measured with HST /FGS, HST /WFC3 spatial scanning and/or Hipparcos , and (iv) 2 DEBs in M31. The Hubble constant from each is 72.25 ± 2.51, 72.04 ± 2.67, 76.18 ± 2.37, and 74.50 ± 3.27 km s{sup 1} Mpc{sup 1}, respectively. Our best estimate of H {sub 0} = 73.24 ± 1.74 km s{sup 1} Mpc{sup 1} combines the anchors NGC 4258, MW, and LMC, yielding a 2.4% determination (all quoted uncertainties include fully propagated statistical and systematic components). This value is 3.4 σ higher than 66.93 ± 0.62 km s{sup 1} Mpc{sup 1} predicted by ΛCDM with 3 neutrino flavors having a mass of 0.06 eV and the new Planck data, but the discrepancy reduces to 2.1 σ relative to the prediction of 69.3 ± 0.7 km s{sup 1} Mpc{sup 1} based on the comparably precise combination of WMAP +ACT+SPT+BAO observations, suggesting that systematic uncertainties in CMB radiation measurements may play a role in the tension. If we take the conflict between Planck high-redshift measurements and our local determination of H {sub 0} at face value, one plausible explanation could involve an additional source of dark radiation in the early universe in the range of Δ N {sub eff} ≈ 0.4–1. We anticipate further significant improvements in H {sub 0} from upcoming parallax measurements of long-period MW Cepheids.

79 ASTRONOMY AND ASTROPHYSICS↗

Joint analysis of galaxy-galaxy lensing and galaxy clustering: Methodology and forecasts for Dark Energy Survey

The joint analysis of galaxy-galaxy lensing and galaxy clustering is a promising method for inferring the growth function of large scale structure. Our analysis will be carried out on data from the Dark Energy Survey (DES), with its measurements of both the distribution of galaxies and the tangential shears of background galaxies induced by these foreground lenses. We develop a practical approach to modeling the assumptions and systematic effects affecting small scale lensing, which provides halo masses, and large scale galaxy clustering. Introducing parameters that characterize the halo occupation distribution (HOD), photometric redshift uncertainties, and shear measurement errors, we study how external priors on different subsets of these parameters affect our growth constraints. Degeneracies within the HOD model, as well as between the HOD and the growth function, are identified as the dominant source of complication, with other systematic effects sub-dominant. The impact of HOD parameters and their degeneracies necessitate the detailed joint modeling of the galaxy sample that we employ. Finally, we conclude that DES data will provide powerful constraints on the evolution of structure growth in the universe, conservatively/optimistically constraining the growth function to 7.9%/4.8% with its first-year data that covered over 1000 square degrees, and to 3.9%/2.3% with its full five-year data that will survey 5000 square degrees, including both statistical and systematic uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Nucleon momentum fraction, helicity and transversity from 2+1-flavor lattice QCD

A detailed analysis of the systematic uncertainties in the calculation of the isovector momentum fraction, ( x ) u - d , helicity moment, ( x ) Δ u - Δ d , and the transversity moment, ( x ) δu - δd , of the nucleon is presented using high-statistics data on seven ensembles of gauge configurations generated by the JLab/W&M/LANL/MIT collaborations using 2 + 1-flavors of dynamical Wilson-clover quarks. The much higher statistics have facilitated better control over all systematics compared to previous lattice calculations. The least understood systematic — excited-state contamination — is quantified by studying the variation of the results as a function of different estimates of the mass gap of the first excited state, obtained from two- and three-point correlation functions, and as a function of the pion mass M π . The final results are obtained using a simultaneous fit in the lattice spacing a , pion mass M π and the finite volume parameter M π L keeping leading order corrections. The data show no significant dependence on the lattice spacing and some evidence for finite-volume corrections. Our final results, in the $ \overline{\mathrm{MS}} $ scheme at 2 GeV, are ( x ) u - d = 0 . 155(17)(20), ( x ) Δ u - Δ d = 0 . 183(14)(20) and ( x ) δu - δd = 0 . 220(18)(20), where the first error is the overall analysis uncertainty assuming excited-state contributions have been removed, and the second is an additional systematic uncertainty due to possible residual excited-state contributions. These results are consistent with phenomenological global fit values.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Volume-preserving algorithm for secular relativistic dynamics of charged particles

Secular dynamics of relativistic charged particles has theoretical significance and a wide range of applications. However, conventional algorithms are not applicable to this problem due to the coherent accumulation of numerical errors. To overcome this difficulty, we develop a volume-preserving algorithm (VPA) with long-term accuracy and conservativeness via a systematic splitting method. Applied to the simulation of runaway electrons with a time-span over 10 magnitudes, the VPA generates accurate results and enables the discovery of new physics for secular runaway dynamics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Study of self-interaction errors in density functional predictions of dipole polarizabilities and ionization energies of water clusters using Perdew-Zunger and locally scaled self-interaction corrected methods

Here, we studied the effects of self-interaction error (SIE) on the static dipole polarizabilities of water clusters modelled with three increasingly sophisticated, non-empirical density functional approximations (DFAs), viz. the local spin density approximation (LDA), the Perdew-Burke-Ernzherof (PBE) generalized-gradient approximation (GGA), and the strongly constrained and appropriately normed (SCAN) meta-GGA, using the Perdew-Zunger self-interaction-correction (PZ-SIC) energy functional in the Fermi-Lowdin orbital SIC (FLO-SIC) framework. Our results show that while all three DFAs overestimate the cluster polarizabilities, the description systematically improves from LDA to PBE to SCAN. The self-correlation free SCAN predicts polarizabilities quite accurately with a mean absolute error (MAE) of 0.58 Bohr 3 with respect to coupled cluster singles and doubles (CCSD) values. Removing SIE using PZ-SIC correctly reduces the DFA polarizabilities, but over-corrects, resulting in underestimated polarizabilities in SIC-LDA, -PBE, and -SCAN. Finally, we applied a recently proposed local-scaling SIC (LSIC) method using a quasi self-consistent scheme and using the kinetic energy density ratio as an iso-orbital indicator. The results show that the LSIC polarizabilities are in excellent agreement with mean absolute error of 0.08 Bohr 3 for LSIC-LDA and 0.06 Bohr 3 for LSIC-PBE with most recent CCSD polarizabilities. Likewise, the ionization energy estimates as an absolute of highest occupied energy eigenvalue predicted by LSIC are also in excellent agreement with CCSD(T) ionization energies with MAE of 0.4 eV for LSIC-LDA and 0.04 eV for LSIC-PBE. The LSIC-LDA predictions of ionization energies are comparable to the reported GW ionization energies while the LSIC-PBE ionization energies are more accurate than reported GW results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Log-normal Star Formation Histories in Simulated and Observed Galaxies

Gladders et al. have recently suggested that the star formation histories (SFHs) of individual galaxies are characterized by a log-normal function in time, implying a slow decline rather than rapid quenching. We test their conjecture on theoretical SFHs from the cosmological simulation Illustris and on observationally inferred SFHs. While the log-normal form necessarily ignores short-lived features such as starbursts, it fits the overall shape of the majority of SFHs very well. In particular, 85% of the cumulative SFHs are fitted to within a maximum error of 5% of the total stellar mass formed, and 99% to within 10%. The log-normal performs systematically better than the commonly used delayed-τ model, and is superseded only by functions with more than three free parameters. Poor fits are mostly found in galaxies that were rapidly quenched after becoming satellites. We explore the log-normal parameter space of normalization, peak time, and full width at half maximum, and find that the simulated and observed samples occupy similar regions, though Illustris predicts wider, later-forming SFHs on average. The ensemble of log-normal fits correctly reproduces complex metrics such as the evolution of Illustris galaxies across the star formation main sequence, but overpredicts their quenching timescales. SFHs in Illustris are a diverse population not determined by any one physical property of galaxies, but follow a tight relation, where width ∝ (peaktime){sup 3/2}. We show that such a relation can be explained qualitatively (though not quantitatively) by a close connection between the growth of dark matter halos and their galaxies.

79 ASTRONOMY AND ASTROPHYSICS↗

Limitations of reanalysis data for wind power applications

Wind energy resource estimates commonly depend on simulated wind speed profiles generated by reanalysis or weather models due to the lack of long time series measurements with sufficient coverage at relevant heights (roughly 90 m above ground). However, modeled data, including reanalyses, can be noisy and display a wide range of biases and errors, variously attributed to terrain effects, poor coverage of assimilated inputs, and model resolution. Wind generation records, if available at high temporal and geographical resolution, can provide a proxy for wind measurements and allow for evaluation of reanalyses and weather model wind time series. We use a 7-year-long data set of hourly, plant-level generation records from over 100 wind plants across Texas to evaluate two commonly used reanalysis data sets (MERRA2 and ERA5). Additionally, we use 1-year of records (2019) to evaluate an operational, high-resolution regional weather modeling product (HRRR v3). We find that across the region, and across all modeling products, the modeled representation of wind generation (i.e., wind speeds at hub heights passed through a power curve) has relatively small mean errors when aggregated daily, but that accuracy and hourly correlation have a strong diurnal sensitivity. Accuracy and correlation systematically decline through the evening and markedly improve after sunrise. These diurnal patterns persist even in the highest resolution model tested (HRRR v3). We hypothesize the nighttime decline in accuracy is mostly due to poorly represented boundary layer conditions, perhaps related to model representation of stability, while other uncertainties (such as wake effects) play a secondary role.

17 WIND ENERGY↗

A perspective on single-crystal layered oxide cathodes for lithium-ion batteries

As the demand for lithium-ion batteries grows exponentially to feed the nascent electric-vehicle and grid-storage markets, the need for higher energy density and longer cycle life becomes more apparent. Increasing the nickel content in the layered-oxide cathodes has been a dominant strategy to increase energy density, but this has exacerbated the surface reactivity concerns. Furthermore, the high states of charge associated with increased capacity lead to cracking of the polycrystalline cathode particles, which exposes fresh surfaces and accelerates capacity fade. Single-crystal cathode particles with low susceptibility to cracking have recently taken on intense interest due to their remarkable cyclability. Many different approaches to single-crystal synthesis have been explored, but there has been no systematic analysis of the different techniques to date. Additionally, many reports have taken on a trial-and-error approach, ignoring well established literature on grain growth outside the battery field. This perspective reviews the literature on single-crystal cathode synthesis through the lens of classic grain growth theories, providing insight and hopefully accelerating the adoption of the promising single-crystal cathode morphology.

25 ENERGY STORAGE↗

How to avoid multiple scattering in strongly scattering SANS and USANS samples

Small Angle Neutron Scattering (SANS) and Ultra Small Angle Neutron Scattering (USANS) are the only available experimental techniques to provide seamless non-destructive measurements of the geometry of the accessible and inaccessible pore structure of rocks from sub-nanopore size to the scale of macropores. They have therefore become the measurement of choice for tight reservoir rocks such as organic rich shales. A simplifying assumption in the analysis is, however, that during the path of neutrons through the sample each neutron is only scattered once. Shales are samples with a high scattering power and Multiple Scattering (MS) may occur which requires special modelling for deconvolution of the results. The approach to avoid MS is to simply reduce the sample thickness to <0.15–0.5 mm. Here, in this work, we present a systematic method on wavelength selection and preparation of samples to optimise extraction of microstructural data and minimise parasitic errors. Experimentally measured SAS transmission (TSAS) values are used as a practical criterion for estimation of the extent of MS. Generous beamtime allocations allowed robust testing revealing that sample thicknesses can be twice as thick as predicted using the standard protocol. Analysing thicker samples is particularly beneficial for statistically relevant characterisation of heterogeneous samples making the new protocol the method of choice for such samples.

(U)SANS↗

Real-Space Pseudopotential Method for the Calculation of Third-Row Elements X-ray Photoelectron Spectroscopic Signatures

X-ray photoelectron spectroscopy (XPS) is a powerful characterization technique that unveils subtle chemical environment differences via core– electron binding energy (CEBE) analysis. We extend the development of realspace pseudopotential methods to calculating 1s, 2s, and 2p 3/2 CEBEs of third-row elements (S, P, and Si) within the framework of Kohn–Sham density-functional theory (KS-DFT). The new approach systematically prevents variational collapse and simplifies core-excited orbital selection within dense energy level distributions. However, careful error cancellation analysis is required to achieve accuracy comparable to all-electron methods and experiments. Combined with real-space KS-DFT implementation, this development enables large-scale simulations with both Dirichlet boundary conditions and periodic boundary conditions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Energy-Screened Many-Body Expansion for Protein–Ligand Interactions: Examining Convergence for Metalloenzymes Through Seven–Body Interactions

Fragment-based quantum chemistry is a powerful strategy for calculating protein−ligand interaction energies using quantum chemistry methods. Rigorous convergence often requires hundreds of atoms in the protein binding-site model, especially if that model is constructed using distance-based criteria to select amino acid residues, while three- and four-body calculations exhibit instability related to combinatorial proliferation in the number of subsystem calculations. Here, we report an energy-based screening protocol for the many-body expansion applied to protein−ligand interactions, implemented in the open-source FRAGME∩T code. Using a combination of aggressive screening based on semiempirical quantum chemistry, with an improved graph-theoretical algorithm to eliminate unimportant subsystems, we are able to perform n-body calculations up to n = 7 using density functional theory in triple-ζ basis sets. Distance cutoffs further reduce the cost without compromising accuracy. Rapid and stable convergence of the many-body expansion is obtained by n = 4, for a pair of metalloenzymes in which a divalent ion coordinates directly to the ligand. As compared to previous results that relied solely on distance cutoffs, oscillations in the n-body corrections are reduced or eliminated, although residual errors remain in one case. This work demonstrates that benchmark-quality protein−ligand interaction energies can be systematically converged using a method with excellent parallel efficiency and scalability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

SPITZER OBSERVATIONS OF SHOCK-EXCITED HYDROGEN DEUTERIDE IN IC 443C, HH 7, AND HH 54: PROBING THE GAS-PHASE DEUTERIUM ABUNDANCE IN THE DENSE INTERSTELLAR MEDIUM

We present and analyze observations-carried out using the Infrared Spectrograph (IRS) on the Spitzer Space Telescope-of the R(3) and R(4) pure rotational lines of hydrogen deuteride (HD) detected from shock-heated material associated with the supernova remnant IC 443C and with the Herbig-Haro objects HH 7 and HH 54. Assuming a continuous temperature distribution for gas observed along the sight lines, we have constrained the gas density to be in the {approx}10{sup 3}-10{sup 4} cm{sup -3} range, using both spectroscopic data for H{sub 2}, HD, and CO from IRS and from the Infrared Space Observatory (ISO), as well as photometric data from Spitzer's Infrared Array Camera. The derived HD abundance relative to H{sub 2} is quite sensitive to the assumed excitation conditions in the emitting gas. Assuming that HD accounts for all gas-phase deuterium in the emitting material, and using all the available spectroscopic data to constrain the excitation conditions, we obtained gas-phase deuterium abundances [D/H]{sub gas} of 0.95{sup +0.54}{sub -0.27} Multiplication-Sign 10{sup -5} and 0.87{sup +0.31}{sub -0.27} Multiplication-Sign 10{sup -5} (statistical errors only) for IC 443C and HH 54, respectively. The uncertainties in the HD abundance are dominated by systematic effects related to the poorly known excitation conditions, and more accurate estimates of the HD abundance in shocked molecular clouds will require measurements of the emission in additional HD rotational transitions.

79 ASTRONOMY AND ASTROPHYSICS↗