SEARCH · Search NASA
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.
The report focuses on the Chi-Nu experiment and provides recent results with covariances for 235 U and 239 Pu and correlations between results and PFNS ratios.
A method for estimating covariance matrices which capture the uncertainties in calculated reactor spectra has been developed. This method is based on perturbing the parameters of a physics-based analytic model fitted to a calculated spectrum. The covariance of the perturbed analytic spectra imposes energy-dependent correlations due to the physics of the neutron processes in the reactor, i.e., a fission component, a 1/E down-scatting component, and a thermal Maxwellian component. An analytic model is developed which is shown to produce good fits to several reactor environments. The covariance matrices produced via this method are then used as the prior spectrum in STAYSL least squares spectrum adjustment where it is combined with integral metrics, such as activation measurements, to produce a high-fidelity neutron spectrum characterization. It was concluded that the methodology showed agreeable results for the ACRR free-field spectrum adjustment in STAYSL resulting in a 𝜒 2 value of 2.21 (per degree of freedom), but further work is needed to describe scattering and interface regions.
Fission product yields (FPY) are important for a variety of applications (reactor neutronics, spent fuel, dosimetry, radiochemistry, etc.) and are currently included in many of the evaluated libraries around the world. The FPYs in the current US evaluation, ENDF/B-VIII.0, are mainly based on the 1994 evaluation of England and Rider and have only had slight updates—such as the inclusion of a 2 MeV point for 239 Pu—since their development. Additionally, only mean values and uncertainties are included in the evaluation, not full correlations. Los Alamos National Laboratory, in collaboration with several other institutes, has been working on an updated evaluation for the FPYs of 239 Pu(n,f), 235 U(n,f), 238 U(n,f), and 252 Cf(sf) using the deterministic, Hauser-Feshbach, fission fragment decay code, BeoH. BeoH calculates the FPYs consistently with many other prompt and delayed fission observables, explicitly taking into account multi-chance fission and ensuring consistency between observables. In addition to providing updated means and uncertainties for the FPYs on a pointwise energy grid from thermal to 20 MeV, we calculate correlations between all FPYs at each incident energy and across incident energies. Here, we discuss the development of these covariance matrices, differences in the correlations between FPYs based on the parameters that are included in the model optimization, and correlations across incident energies for neutron-induced fission.
Computation of correlated ionic transport properties from molecular dynamics in the Green-Kubo formalism is expensive, as one cannot rely on the affordable mean square displacement approach. We use spectral decomposition of the short-time ionic displacement covariance to learn a set of diffusion eigenmodes that encode the correlation structure and form a basis for analyzing the ionic trajectories. This allows systematic reduction of the uncertainty and accelerate computations of ionic conductivity in systems with a steady-state correlation structure. We provide mathematical and numerical proofs of the method’s robustness and demonstrate it on realistic electrolyte materials.
We study conformal conserved currents in arbitrary irreducible representations of the Lorentz group using the embedding space formalism. With the help of the operator product expansion, we first show that conservation conditions can be fully investigated by considering only two- and three-point correlation functions. We then find an explicitly conformally-covariant differential operator in embedding space that implements conservation based on the standard position space operator product expansion differential operator ∂μ, although the latter does not uplift to embedding space covariantly. The differential operator in embedding space that imposes conservation is the same differential operator $D_{ijA}$ used in the operator product expansion in embedding space. We provide several examples including conserved currents in irreducible representations that are not symmetric and traceless. With an eye on four-point conformal bootstrap equations for four conserved vector currents $\langle JJJJ\rangle$ and four energy-momentum tensors $\langle TTTT\rangle$, we mostly focus on conservation conditions for $\langle JJ\mathcal{O}\rangle$ and $\langle TT\mathcal{O}\rangle$. Finally, we reproduce and extend the consequences of conformal Ward identities at coincident points by determining three-point coefficients in terms of charges.
Analytic formulas to calculate statistical uncertainties of v n { 2 k } harmonics extracted from the Q cumulants are presented. The Q cumulants are multivariate polynomial functions of the weighted means of 2m-particle azimuthal correlations, $\langle\langle 2 m \rangle\rangle$. Variances and covariances of the $\langle\langle 2 m \rangle\rangle$ are included in the analytic formulas of the uncertainties that can be calculated simultaneously with the calculations of the v n { 2 k } harmonics. The calculations are performed using a simple toy model, which roughly simulates elliptic flow azimuthal anisotropy with magnitudes around 0.05. The results are compared with the results obtained by the many data subsets, and by the bootstrapping method. The first one is a common way of estimation of the statistical uncertainties of the v n { 2 k } harmonics in a real experiment. In order to increase precision in the measurement of the v n { 2 k } harmonics, a large number of 15 000 subsets and the same number of the resampling in the bootstrap method is used. Unlike the other ways of the analytic calculation of the v n { 2 k } statistical uncertainties, our proposal that includes the use of squared weights in the calculation of both the variances and covariances, gives the best agreement with the results obtained from the subsets and bootstrap method. Additionally, a recurrence equation between Q cumulants of any order is also presented.
Lead exposures from legacy sources threaten children’s health. Soil in Omaha, Nebraska, was contaminated by emissions from a lead smelter and refinery. The U.S. Environmental Protection Agency excavated and replaced contaminated soil at the Omaha Lead Superfund Site between 1999 and 2016. The goal of this study was to assess the association of soil lead level (SLL) and soil remediation status with blood lead levels (BLLs) in children living near or on the site. We linked information on SLL at residential properties with children’s BLLs and assigned remediation status to children’s BLL measurements based on whether their measurements occurred during residence at remediated or unremediated properties. We examined the association of SLL and remediation status with elevated BLL (EBLL). We distinguished the roles of temporal trend and the intervention with time-by-intervention-status interaction contrasts. All analyses estimated odds ratios (ORs) with a generalized estimating equations approach to ensure robustness under the complex correlations among BLL measurements. All analyses controlled for relevant covariates including children’s characteristics. EBLL (> 5 μg/dL ) was associated with both residential SLL [e.g., OR=2.00; 95% confidence interval (CI): 1.83, 2.19; > 400–800 vs. ≤ 200 ppm] and neighborhood SLL [e.g., OR=1.85 (95% CI: 1.62, 2.11; > 400–800 vs. ≤ 200 ppm)] before remediation but only with neighborhood SLL after remediation. The odds of EBLL were higher before remediation [OR 1.52 (95% CI: 1.34, 1.72)]. Similarly, EBLL was positively associated with preremediation status in our interaction analysis [interaction OR=1.18 (95%CI: 1.02, 1.37)]. Residential and neighborhood SLLs were important predictors of EBLLs in children residing near or on this Superfund site. Neighborhood SLL remained a strong predictor following remediation. Our data analyses showed the benefit of soil remediation. Results from the interaction analyses should be interpreted cautiously due to imperfect correspondence of remediation times between remediation and comparison groups. https://doi.org/10.1289/EHP8657
Satellite quantification of aerosol effects on clouds relies on aerosol optical depth (AOD) as a proxy for aerosol concentration or cloud condensation nuclei (CCN). However, the lack of error characterization of satellite-based results hampers their use for the evaluation and improvement of global climate models. We show that the use of AOD for assessing aerosol–cloud interactions (ACIs) is inadequate over vast oceanic areas in the subtropics. Instead, we postulate that a more physical approach that consists of matching vertically resolved aerosol data from the Cloud-Aerosol Lidar and Infrared Pathfinder Satellite Observations (CALIPSO) satellite at the cloud-layer height with Moderate Resolution Imaging Spectroradiometer (MODIS) Aqua cloud retrievals reduces uncertainties in satellite-based ACI estimates. Combined aerosol extinction coefficients (σ) below cloud top (σBC) from the Cloud-Aerosol Lidar with Orthogonal Polarization (CALIOP) and cloud droplet number concentrations (N d ) from MODIS Aqua yield high correlations across a broad range of σBC values, with σBC quartile correlations ≥0.78. In contrast, CALIOP-based AOD yields correlations with MODIS Nd of 0.54–0.62 for the two lower AOD quartiles. Moreover, σBC explains 41% of the spatial variance in MODIS N d , whereas AOD only explains 17%, primarily caused by the lack of spatial covariability in the eastern Pacific. Compared with σBC, near-surface σ weakly correlates in space with MODIS N d , accounting for a 16% variance. It is concluded that the linear regression calculated from ln(N d )–ln(σBC) (the standard method for quantifying ACIs) is more physically meaningful than that derived from the Nd–AOD pair.
Thomson scattering measurements in high energy density experiments are often recorded using optical streak cameras. In the low-signal regime, noise introduced by the streak camera can become an important and sometimes the dominant source of measurement uncertainty. In this paper, we present a formal method of accounting for the presence of streak camera noise in our measurements. We present a phenomenological description of the noise generation mechanisms and present a statistical model that may be used to construct the covariance matrix associated with a given measurement. This model is benchmarked against simulations of streak camera images. We demonstrate how this covariance may then be used to weight fitting of the data and provide quantitative assessments of the uncertainty in the fitting parameters determined by the best fit to the data and build confidence in the ability to make statistically significant measurements in the low-signal regime, where spatial correlations in the noise become apparent. These methods will have general applicability to other measurements made using optical streak cameras.
We study, for the first time, the cross correlation between the angular distribution of radial peculiar velocities (PV) and the lensing convergence of cosmic microwave background (CMB) photons. We derive theoretical expectations for the signal and its covariance and assess its detectability with existing and forthcoming surveys. We find that such cross-correlations are expected to improve constraints on different gravitational models by partially breaking degeneracies with the matter density. We identify in the distance-scaling dispersion of the peculiar velocities the most relevant source of noise in the cross correlation. For this reason, we also study how the above picture changes assuming a redshift-independent scatter for the PV, obtained for example using a reconstruction technique. Our results show that the cross correlation might be detected in the near future combining PV measurements from DESI and the convergence map from CMB-S4. Using realistic direct PV measurements we predict a cumulative signal-to-noise ratio of approximately 3.8σ using data on angular scales 3 ≤ ℓ ≤ 200. For an idealized reconstructed peculiar velocity map extending up to redshift z = 0.15 and a smoothing scale of 4 Mpc h -1 we predict a cumulative signal-to-noise ratio of approximately 27σ from angular scales 3 ≤ ℓ ≤ 200. We conclude that currently reconstructed peculiar velocities have more constraining power than directly observed ones, even though they are more cosmological-model dependent.
m, n = 2, 1 tearing mode onset empirical probability and machine learning analyses of a multiscenario DIII-D database of over 14 000 H-mode discharges show that the normalized plasma beta, the rotation profile, and the magnetic equilibrium shape have the strongest impact on the 2,1 tearing mode stability, in qualitative agreement with neoclassical tearing modes (m and n are the poloidal and toroidal mode numbers, respectively). In addition, 2,1 tearing modes are most likely to destabilize when n > 1 tearing modes are already present in the core plasma. The covariance matrix of tearing sensitive plasma parameters takes a nearly block-diagonal form, with the blocks incorporating thermodynamic, current and safety factor profile, separatrix shape, and plasma flow parameters, respectively. This suggests a number of paths to improved stability at fixed pressure and edge safety factor primarily by preserving a minimum of 1 kHz differential rotation, increasing the minimum safety factor above unity, using upper single null magnetic configuration, and reducing the core impurity radiation. In addition, lower triangularity, lower elongation, and lower pedestal pressure may also help to improve stability. Furthermore, the electron and ion temperature, collisionality, resistivity, internal inductance, and the parallel current gradient appear to only weakly correlate with the 2,1 tearing mode onsets in this database.
Nuclear data are necessary for reliable modeling and simulation of the next generation of nuclear reactors. However, the variation in the ratio of the computed to experimental values (C/E) for certain types of nuclear systems is much less than predicted by evaluated nuclear data file (ENDF)/B covariances. Figure 1 provides an example of a set of metal-plutonium-fueled, fast-spectrum (PU-MET-FAST) integral experiments from the International Criticality Safety Benchmark Evaluation Project (ICSBEP) in the Oak Ridge National Laboratory (ORNL) VALID database. The variation in the C/E values, shown with one standard deviation error bars, is 100% covered by both the SCALE and ENDB/VIII.0 covariance data. This is due to the comparisons to integral data that are essential during the evaluation process. However, the ENDF evaluations represent uncertainties and correlations in differential data only; they do not reflect the impact of the comparison to integral data in the covariance evaluations.
Knowledge of the distribution of soil thermal properties is important for understanding subsurface hydrological and biogeochemical processes. This study describes and evaluates quick thermal profiling (QTP), a new measurement technique aimed at providing rapid, depth-resolved measurements of soil thermal inertia at numerous locations across the landscape. A cylindrical probe with temperature sensors at multiple depths is quickly inserted into the ground, and soil thermal inertia is estimated from how quickly the probe temperature equilibrates with the soil. To this end, a finite volume heat transfer model is used to generate temperature equilibration time series across combinations of controlling factors, and a gridded search inversion approach is applied to infer soil thermal inertia. Field tests in the Arctic indicate that QTP measurements have a minimum uncertainty of 0.14 J m −2 K −1 s −1/2 and covary with dual-probe heat pulse thermal analyzer measurements (concordance correlation coefficient = 0.56) with a root-mean-square error of 0.40 J m −2 K −1 s −1/2 . Besides demonstrating the value of QTP for estimating thermal inertia, this study identifies various sources of measurement uncertainty, particularly probe-soil contact resistance and frictional heating. Further, analysis of soil samples indicates that thermal inertia can be used to estimate thermal conductivity and dry bulk density in the studied area, although such inferences are highly site-specific. Overall, the QTP method holds promise to generate thermal inertia data products and to complement other characterization approaches for advancing understanding of soil properties across far more locations than is currently possible.
We present a Gibbs sampling solution to the mapmaking problem for cosmic microwave background (CMB) measurements that builds on existing destriping methodology. Gibbs sampling breaks the computationally heavy destriping problem into two separate steps: noise filtering and map binning. Considered as two separate steps, both are computationally much cheaper than solving the combined problem. This provides a huge performance benefit as compared to traditional methods and it allows us, for the first time, to bring the destriping baseline length to a single sample. Here, we applied the Gibbs procedure to simulated Planck 30 GHz data. We find that gaps in the time-ordered data are handled efficiently by filling them in with simulated noise as part of the Gibbs process. The Gibbs procedure yields a chain of map samples, from which we are able to compute the posterior mean as a best-estimate map. The variation in the chain provides information on the correlated residual noise, without the need to construct a full noise covariance matrix. However, if only a single maximum-likelihood frequency map estimate is required, we find that traditional conjugate gradient solvers converge much faster than a Gibbs sampler in terms of the total number of iterations. The conceptual advantages of the Gibbs sampling approach lies in statistically well-defined error propagation and systematic error correction. This methodology thus forms the conceptual basis for the mapmaking algorithm employed in the BEYONDPLANCK framework, which implements the first end-to-end Bayesian analysis pipeline for CMB observations.
There are new evaluation methods being explored for combined RRR + URR. Simple example demonstrated in SAMMY code using. Potential to be applied to future LoFi, Medium-Fi covariance efforts. Much more work needs to be done to fully correlate regions.
The generality of the worldwide leaf economics spectrum (LES) has made it a pillar of trait-based ecological research. Yet, few studies have examined the processes shaping the evolution of the LES within species, in part, because most species occupy only a small portion of the LES. Here, to address this gap, we took advantage of the distinct leaf economics strategies present in different ecotypes of the phenotypically diverse perennial grass Panicum virgatum (switchgrass) to generate a genetic mapping population, which we planted in common gardens at three sites spanning 12 degrees of latitude in the central United States. With this genetic mapping population, we evaluated two potentially interacting causes of LES evolution: 1) genetic architecture, where multiple traits are influenced by either the same gene (pleiotropy) or by genes in close physical proximity (genetic linkage), and 2) correlational selection, where selection acts on traits in combination rather than in isolation. We found that shared genetic architecture influenced covariation between photosynthetic rate (A MASS ) and leaf nitrogen (N MASS ) and between A MASS and leaf mass per area (LMA). We also found that correlational selection favored the trait combinations predicted by the LES (e.g., high LMA with low N MASS or low LMA with high N MASS ) and disfavored other, mismatched trait combinations at two of the three sites. Together, these results demonstrate how the evolution of an integrated LES within species can arise from multiple evolutionary causes.
Abstract The increasing statistical power of cosmic microwave background (CMB) datasets requires a commensurate effort in understanding their noise properties. The noise in maps from ground-based instruments is dominated by large-scale correlations, which poses a modeling challenge. This paper develops novel models of the complex noise covariance structure in the Atacama Cosmology Telescope Data Release 6 (ACT DR6) maps. We first enumerate the noise properties that arise from the combination of the atmosphere and the ACT scan strategy. We then prescribe a class of Gaussian, map-based noise models, including a new wavelet-based approach that uses directional wavelet kernels for modeling correlated instrumental noise. The models are empirical, whose only inputs are a small number of independent realizations of the same region of sky. We evaluate the performance of these models against the ACT DR6 data by drawing ensembles of noise realizations. Applying these simulations to the ACT DR6 power spectrum pipeline reveals a ∼ 20% excess in the covariance matrix diagonal when compared to an analytic expression that assumes noise properties are uniquely described by their power spectrum. Along with our public code,mnms, this work establishes a necessary element in the science pipelines of both ACT DR6 and future ground-based CMB experiments such as the Simons Observatory (SO).
Abstract When the vortex center location is estimated from a radar-scanned tornadic mesocyclone, the estimated location is not error-free. This raises an important issue concerning the sensitivities of analyzed vortex flow (VF) fields by the VF-Var (formulated in Part I of this paper series and tested in Part II) to vortex center location errors, denoted by Δ x c . Numerical experiments are performed to address this issue with the following findings: The increase of |Δ x c | from zero to a half of vortex core radius causes large analysis error increases in the vortex core but the increased analysis errors decrease rapidly away from the vortex core especially for dual-Doppler analyses. The increased horizontal-velocity errors in the vortex core are mainly in the Δ x c -normal component, because this component varies much more rapidly than the other component along the Δ x c direction in the vortex core. The vertical variations of Δ x c distort the vertical correlation structure of Δ x c -dislocated VF-dependent background error covariance, which can increase the analysis errors in the vortex core. The dual-Doppler analyses have adequate accuracies outside the vortex core even when |Δ x c | increases to a half of vortex core radius, while single-Doppler analyses can also have adequate accuracies outside the vortex core mainly for the single-Doppler-observed velocity component. The sensitivities to Δ x c are largely unaffected by the vortex slanting. The above findings are important and useful for assessing the accuracies of analyzed VFs for real radar-observed tornadic mesocyclones. Significance Statement When the vortex center location is estimated from a radar-scanned tornadic mesocyclone, the estimated location is not error-free. This raises an issue concerning the sensitivity of analyzed vortex flow (VF) by the VF-Var (formulated in Part I of this paper series and tested with simulated radar observations in Part II) to vortex center location error. This issue and its required investigations are very important for the VF-Var to be applied to real radar-observed tornadic mesocyclones, especially in an operational setting with the WSR-88Ds. Numerical experiments are performed to address this issue. The findings from these experiments are important and useful for assessing the accuracies of VF-Var analyzed VF fields for real radar-observed tornadic mesocyclones.