Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distribution parameters”

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 127 records · Page 7

Uncertainty quantification of an empirical shell-model interaction using principal component analysis

Recent investigations have emphasized the importance of uncertainty quantification (UQ) in nuclear theory. Here, we carry out UQ for configuration-interaction shell-model calculations in the 1$\textit{s}$–0$\textit{d}$ valence space, investigating the sensitivity of observables to perturbations in the 66 parameters (matrix elements) of a high-quality empirical interaction. The large parameter space makes computing the corresponding Hessian numerically costly, so we compare a cost-effective approximation, using the Feynman-Hellmann theorem, to the full Hessian and find it works well. Diagonalizing the Hessian yields the principal components of the interaction: linear combinations of parameters ordered by sensitivity. This approximately decoupled distribution of parameters facilitates theoretical uncertainty propagation onto structure observables: electromagnetic transitions, Gamow-Teller decays, and dark matter-nucleus scattering matrix elements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Advancing earth system model calibration: a diffusion-based method

Understanding of climate impact on ecosystems globally requires site-specific model calibration. Here we introduce a novel diffusion-based uncertainty quantification (DBUQ) method for efficient model calibration. DBUQ is a score-based diffusion model that leverages Monte Carlo simulation to estimate the score function and evaluates a simple neural network to quickly generate samples for approximating parameter posterior distributions. DBUQ is stable, efficient, and can effectively calibrate the model given diverse observations, thereby enabling rapid and site-specific model calibration on a global scale. This capability significantly advances Earth system modeling and our understanding of climate impacts on Earth systems. We demonstrate DBUQ's capability in E3SM land model calibration at the Missouri Ozark AmeriFlux forest site. Both synthetic and real-data applications indicate that DBUQ produces accurate parameter posterior distributions similar to those generated by Markov Chain Monte Carlo sampling but with 30X less computing time. This efficiency marks a significant stride in model calibration, paving the way for more effective and timely climate impact analyses.

Liu, Yanfang↗

Drop Size Distribution Variability in Central Argentina during RELAMPAGO-CACTI

The Remote sensing of Electrification, Lightning, And Meso-scale/micro-scale Processes with Adaptive Ground Observations (RELAMPAGO) and the Cloud, Aerosol, and Complex Terrain Interactions Experiment Proposal (CACTI) field campaigns provided an unprecedented thirteen-disdrometer dataset in Central Argentina during the Intensive (IOP, 15 November to 15 December 2018) and Extended (EOP, 15 October 2018 to 30 April 2019) Observational Periods. The drop size distribution (DSD) parameters and their variability were analyzed across the region of interest, which was divided into three subregions characterized by the differing proximity to the Sierras de Córdoba (SDC), in order to assess the impact of complex terrain on the DSD parameters. A rigorous quality control of the data was first performed. The frequency distributions of DSD-derived parameters were analyzed, including the normalized intercept parameter (logN w ), the mean volume diameter (D 0 ), the mean mass diameter (D m ), the shape parameter (μ), the liquid water content (LWC), and the rain rate (R). The region closest to the SDC presented higher values of logN w , lower D 0 , and higher μ, while the opposite occurred in the farthest region, i.e., the concentration of small drops decreased while the concentration of bigger drops increased with the distance to the east of the SDC. Furthermore, the region closest to the SDC showed a bimodal distribution of D 0 : the lower values of D 0 were associated with higher values of logN w and were found more frequently during the afternoon, while the higher D 0 were associated with lower logN w and occurred more frequently during the night. The data were analyzed in comparison to the statistical analysis of Dolan et al. 2018 and sorted according to the classification proposed in the cited study. The logN w -D 0 and LWC-D 0 two-dimensional distributions allowed further discussion around the applicability of other mid-latitude and global precipitation classification schemes (startiform/convection) in the region of interest. Finally, three precipitation case studies were analyzed with supporting polarimetric radar data in order to relate the DSD characteristics to the precipitation type and the microphysical processes involved in each case.

54 ENVIRONMENTAL SCIENCES↗

A multiscale approach to analyze the effect of radial void fraction distributions on BWR neutronic lattice parameters

The goal of this study is to analyze the effect of radial void fraction distribution, considered at three degrees of fidelity, on lattice and 2-group homogenized cross sections. The novelty of this work lies in its use of experimentally-determined radial void fraction distributions which allows for fine-scale coolant density data with no predictive uncertainty and which could be prohibitively expensive to simulate over the range of operational parameters used in the experiments. The experimentally-determined radial void fraction distributions are also characterized using a statistical tool called a semi-variogram in order to determine the scale of continuity. Void fraction distributions are averaged at the lattice and subchannel level for comparison with the fine-scale model. Models averaged at the subchannel and lattice level show a difference in about 68% of the time. Also, by comparing 2-group homogenized cross-sections, the differences between the subchannel and homogeneous cases are measurable.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Multi-Parameter Optical Fiber for Distributed Sensing of Humidity, CH4, CO2, and Corrosion

This work describes the use of the optical fiber sensor (OFS) for successful monitoring of humidity, CH4, and CO2 based on the strain produced along the single-mode fiber (SMF) sensor. This is enabled by absorption of H2O/gases onto the commercially available polyacrylate coated jacketed portion of the fiber resulting in a change in strain. Under equilibrium, a differential microstrain was observed along the jacketed portion of the SMF with N2, CH4, and CO2 at different relative humidity (RH) conditions. The strain response of the SMF under various mixed gas composition and different RH conditions were also measured and made calibration curves accordingly. Linear regression and principal component analysis of the strain datasets provided deconvolution of the impact of strain from H2O, N2, CH4, and CO2. Additionally, modified OFS comprised of the Fe coated fiber section was employed to monitor corrosion based on the increase in backscattered intensity amplitude of the light being passed once corrosion of Fe occurs. Also, corrosion of Fe was studied under soil by installing the fiber in soil with protective measures which would prevent the sensor being mechanically disturbed or broken during installation. The corrosion rates were studied by monitoring the rate at which the intensity of backscattered light amplitude attains a steady state value when complete corrosion of Fe with a specific coating thickness occurs.

Mainali, Badri↗

Two-particle Bose-Einstein correlations and their Lévy parameters in PbPb collisions at s N N = 5.02 TeV

Two-particle Bose–Einstein momentum correlation functions are studied for charged-hadron pairs in lead-lead collisions at a center-of-mass energy per nucleon pair of s NN = 5.02 TeV . The data sample, containing 4.27 × 10 9 minimum bias events corresponding to an integrated luminosity of 0.607 nb − 1 , was collected by the CMS experiment in 2018. The experimental results are discussed in terms of a Lévy-type source distribution. The parameters of this distribution are extracted as functions of particle pair average transverse mass and collision centrality. These parameters include the Lévy index or shape parameter α , the Lévy scale parameter R , and the correlation strength parameter λ . The source shape, characterized by α , is found to be neither Cauchy nor Gaussian, implying the need for a full Lévy analysis. Similarly to what was previously found for systems characterized by Gaussian source radii, a hydrodynamical scaling is observed for the Lévy R parameter. The λ parameter is studied in terms of the core-halo model. ©2024 CERN, for the CMS Collaboration 2024 CERN

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Uniqueness and global optimality of the maximum likelihood estimator for the generalized extreme value distribution

The three-parameter generalized extreme value distribution arises from classical univariate extreme value theory and is in common use for analysing the far tail of observed phenomena, yet important asymptotic properties of likelihood-based estimation under this standard model have not been established. In this paper, we prove that the maximum likelihood estimator is global and unique. An interesting secondary result entails the uniform consistency of a class of limit relations in a tight neighbourhood of the true shape parameter.

54 ENVIRONMENTAL SCIENCES↗

Residual mean field model of valence quarks in the nucleon

Abstract We develop a non-perturbative model for valence parton distribution functions (PDFs) based on the mean field interactions of valence quarks in the nucleonic interior. The main motivation for the model is to obtain a mean field description of the valence quarks as a baseline to study the short range quark–quark interactions that generate the high x tail of PDFs. The model is based on the separation of the valence three-quark cluster and residual system in the nucleon. Then the nucleon structure function is calculated within the effective light-front diagrammatic approach introducing nonperturbative light-front valence quark and residual wave functions. Within the model a new relation is obtained between the position, $$x_p$$ x p , of the peak of $$xq_V(x)$$ x q V ( x ) distribution of the valence quark and the effective mass of the residual system, $$m_R$$ m R , in the form: $$x_{p} \approx {1\over 4} (1-{m_R\over m_N})$$ x p ≈ 1 4 ( 1 - m R m N ) at starting $$Q^2$$ Q 2 . This relation explains the difference in the peak positions for d- and u-quarks through the expected difference of residual masses for valence d- and u-quark distributions. The parameters of the model are fixed by fitting the calculated valence quark distributions to the phenomenological PDFs. This allowed us to estimate the overall mean field contribution in baryonic and momentum sum rules for valence d- and u-quarks. Finally, the evaluated parameters of the non-perturbative wave functions of valence 3q-cluster and residual system can be used in calculation of other quantities such as nucleon form factors, generalized partonic and transverse momentum distributions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Complete quantification of parametric uncertainties in ( d , p ) transfer reactions

Deuteron-induced transfer reactions are a popular probe in nuclear structure and nuclear astro physics studies. The interpretation of these transfer measurements relies on reaction theory that takes as input effective interactions between the nucleons and the target nucleus. Previous work quantified the uncertainty associated with the optical potentials between the nucleons and the target. In this study, we extend that work by also including the parameters of the mean field associated with the overlap function of the final bound state, thus obtaining the full parametric uncertainty on transfer observables. We use Bayesian Markov Chain Monte Carlo simulations to obtain parameter posterior distributions. We use elastic-scattering cross sections to constrain the optical potential parameters and use the asymptotic normalization coefficient of the final state to constrain the bound state interaction. We then propagate these posteriors to the transfer angular distributions and obtain confidence intervals for this observable. We study (d,p) reactions on 14 C, 16 O, and 48 Ca at energies in the range $E$ $d$ = 7 - 24 MeV. Our results show a strong reduction in uncertainty by using the asymptotic normalization coefficient as a constraint, particularly for those reactions most sensitive to ambiguities in the mean field. For those reactions, the importance of constraining the bound state interaction is equal to that of constraining the optical potentials. Furthermore, the case of 14 C is an outlier because the cross section is less sensitive to the nuclear interior. When minimal constraints are used on the parameters of the nucleon-target interaction, the 1σ uncertainties on the differential cross sections are large (~ 500-100%). However, if elastic-scattering data and the asymptotic normalization coefficient are used in the analysis, with an error of 10% (5%), this uncertainty reduces to ~ 30% (~ 15%).

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Quantifying the Turbulent Entrainment‐Mixing Processes Based on Z ‐LWC Relationships of Cloud Droplets

Abstract Turbulent entrainment‐mixing processes profoundly influence the relationship between radar reflectivity factor and liquid water content ( Z ‐LWC) of cloud droplets. However, quantification of the entrainment‐mixing mechanisms based on the Z ‐LWC relationship is still lacking. To address this gap, 12,218 entrainment‐mixing cases are simulated using the Explicit Mixing Parcel Model. We examine the variations of the parameters in the power‐law relationship Z = a LWC b , and the relationship between parameter b and homogeneous mixing degree ( ψ ), a measure quantifying entrainment‐mixing processes. The results indicate that parameter b distributes within the range of 1–2, with a positive correlation between parameter b and ψ . The b ‐ ψ relationship is fitted, which connects the Z ‐LWC relationship for various entrainment‐mixing types. The results suggest the potential for employing a remote sensing approach to investigate the entrainment‐mixing mechanisms of non‐precipitating small cumulus/stratocumulus clouds, thereby overcoming the limitations of traditional observational studies that rely solely on aircraft observations.

54 ENVIRONMENTAL SCIENCES↗

Understanding nanoscale structural distortions in Pb(Zr 0.2 Ti 0.8 )O 3 by utilizing X-ray nanodiffraction and clustering algorithm analysis

Hard X-ray nanodiffraction provides a unique nondestructive technique to quantify local strain and structural inhomogeneities at nanometer length scales. However, sample mosaicity and phase separation can result in a complex diffraction pattern that can make it challenging to quantify nanoscale structural distortions. In this work, a k-means clustering algorithm was utilized to identify local maxima of intensity by partitioning diffraction data in a three-dimensional feature space of detector coordinates and intensity. This technique has been applied to X-ray nanodiffraction measurements of a patterned ferroelectric PbZr 0.2 Ti 0.8 O 3 sample. The analysis reveals the presence of two phases in the sample with different lattice parameters. A highly heterogeneous distribution of lattice parameters with a variation of 0.02 Å was also observed within one ferroelectric domain. This approach provides a nanoscale survey of subtle structural distortions as well as phase separation in ferroelectric domains in a patterned sample.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Implicit neural representations for experimental steering of advanced experiments

Scattering measurements using electrons, neutrons, or photons are essential for obtaining microscopic insights into materials. However, limited facility availability and high-dimensional scattering data necessitate more efficient experimental steering techniques. Here, we report a machine learning method that guides scattering data collection and facilitates real-time estimation of model parameters, given a reliable forward model to simulate experimental signals. We employ implicit neural representations as efficient surrogates that link model parameters with simulated spectroscopies. This enables a Bayesian optimal experimental design framework to estimate the probability distributions of parameters from high-dimensional scattering data. We demonstrate the proposed method using inelastic neutron scattering with simulated and real experimental data, highlighting the method’s ability to provide real-time parameter estimation with quantified uncertainties and to deliver informed experimental guidance that reduces experimental time while maximizing scientific output. This approach paves the way for accelerated discoveries in condensed matter through scattering measurements.

36 MATERIALS SCIENCE↗

Predicting Li-Ion Battery Capacity Fade Using Early-Life Data and a Hybrid Data-Driven Gaussian Process-Bayesian Regression Approach

Accurately predicting Li-ion battery capacity trajectories using early-life data can dramatically improve battery-life understandings and be used to rapidly evaluate design/cost/performance trade-offs when developing new battery materials. Accurate early-life predictions enable researchers to quickly iterate over cell designs and material precursor properties without consistently cycling cells to failure. To this end, we present a toolbox that uses a combined Gaussian Process and Bayesian regression approach that capitalizes on signals other than just capacity (e.g., dQ/dV, voltage drops) to rapidly predict capacity-fade trajectories. The prediction tool uses Bayesian regression to fit functional forms, e.g., power law, sigmoids, etc., to predict capacity-fade dynamics. By fitting functional forms, the capacity fade can be interrogated at any point in the future, allowing for early cell-failure prediction. Additionally, Bayesian regression allows for accurate uncertainty estimates that account for cell-to-cell variability (aleatoric uncertainty) and the lack of observation data (epistemic uncertainty). By only using early cycle data to predict the capacity fade trajectory, uncertainty bounds at end-of-life can be extremely large. The large uncertainty bounds are further exacerbated because there is no systematic way to define the prior distribution of the functional forms' parameters. We improve our the predicted trajectory confidence interval of our predicted trajectory using two methods. First, we shows that a small amount of held-out cycling data is sufficientuse some train cells, that have been cycled to failure to derive information regarding the appropriate prior distributions for the functional forms' parameters of the functional form, effectively leading to data-driven priors.. We propose constructing the data-driven priors by first running a Bayesian regression starting with uninformed priors to generate intermediate cell-specific posterior parameter distributions. These posterior distributions are combined using a Ggaussian mixture model for each parameter to create the data-driven priors. These mixture models serve as the data-driven prior distributions for the parameters for. Second, we derive multiple features, e.g., C_dchg 0.5 DoD 0.5, log (|mean(dQ/dV_(w_3-w_0 ) (V)|), etc., from the train cellsheld-out cycling data, identify which the features are that best predicting capacity at early/mid-life cycles, and then create Ggaussian process regression models that are used for predicting capacity at early/mid-life cycles for the test cells (see blue dots with error bars in Fig 1b). Finally, these predicted data-points are used in addition to the actual early cycle data capacity fade to construct the Bayesian regression trajectory for the test cell s. Notably. We note that these two methods are complementary and can be combined with each other. We evaluate the performance of our proposed method on an testing open-source dataset from Iowa State University and Iowa Lakes Community College (ISU-ILCC). This dataset comprises of 251 nickel-manganese-cobalt/graphite Lithium-ion cells that are cycled under 63 different conditions. We compute the mean average percentage error (MAPE) and negative log predictive density (NLPD) to quantify the efficacy of our method. Our initial findings suggest that, when only few observations are available, for test cells, when using only Bayesian regression with uninformed priors, a power law functional provides the most accurate predictions. with very few data points. However, asHowever, a the number of data points increases, a twin sigmoidal function becomes more accurate as the number of observations further increases. We also find that using as little as 10% of the data set towards generating data-driven priors can lead to significant improvement in prediction accuracy when using early cycle data. Lastly, we found that augmenting early-cycle data with Gaussian process-predicted capacity data for Bayesian regression greatly improves the prediction accuracy. We will present a comprehensive comparison of our methods to other methods available in the literature and apply this method to additional battery datasets.

42 ENGINEERING↗

Cosmologically Coupled Compact Objects: A Single-parameter Model for LIGO–Virgo Mass and Redshift Distributions

Abstract We demonstrate a single-parameter route for reproducing higher mass objects as observed in the LIGO–Virgo mass distribution, using only the isolated-binary stellar evolution channel. This single parameter encodes the cosmological mass growth of compact stellar remnants that exceed the Tolman–Oppenheimer–Volkoff limit. Cosmological mass growth appears in known solutions to General Relativity with cosmological boundary conditions. We consider the possibility of solutions with cosmological boundary conditions, which reduce to Kerr on timescales short compared to the Hubble time. We discuss complementary observational signatures of these solutions that can confirm or invalidate their astrophysical relevance.

79 ASTRONOMY AND ASTROPHYSICS↗

Analysis of multipacting threshold sensitivity to the random distributions of the secondary electron yield parameters

The way multipacting develops, depends strongly on the secondary emission property of the surface material. The knowledge of secondary electron yield is crucial for accurate prediction of the multipacting threshold. Variations in secondary electron yield parameters from experimental measurements create uncertainty, stemming from handling and surface preparation, and these uncertainties significantly affect multipacting threshold predictions. Despite their significance, the previous studies on the multipacting phenomenon did not adequately address the effect of an assumed random distribution of the secondary emission parameters on the multipacting threshold. Therefore, this paper aims to provide a comprehensive statistical study on how the different random distributions of the secondary emission parameters and, as a result, the uncertainty in the secondary electron yield affect multipacting thresholds. We focus on three commonly used distributions, namely uniform, normal, and truncated normal distributions, to define the uncertainty of random inputs. We use the chaos polynomial expansion method to determine how much each of the random parameters contributes to the multipacting threshold uncertainty. Additionally, we calculate Sobol sensitivity indices to evaluate the impact of the individual parameters or groups of parameters on the model outputs and study how different random distributions of these parameters affected the Sobol index results.

physics↗

Nucleon helicity generalized parton distribution at physical pion mass from lattice QCD

The generalized parton distributions (GPDs) offer a window on three-dimensional imaging of the nucleon, providing understanding of how the fundamental properties of the nucleon, such as its mass and spin, arise from the underlying quark and gluon degrees of freedom. In this work, we present the first lattice calculation of the nucleon isovector helicity GPD at physical pion mass, using an a≈0.09 fm lattice ensemble with 2+1+1 flavors of highly improved staggered quarks generated by MILC Collaboration. We perform the GPD calculation in Breit frame using averaged nucleon boost momentum P z ≈2.2 GeV with nonzero momentum transfers in [0.2,1.0] GeV 2 . Nonperturbative renormalization in RI/MOM scheme is used to obtain the quasi-distribution before matching to the lightcone GPDs. The three-dimensional distribution $\tilde{H}$(x, Q 2 ) is presented, along with the three-dimensional nucleon tomography and impact-parameter-dependent distribution for selected Bjorken x at μ=3 GeV in $\overline{MS}$ scheme.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗