Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Gaussian distribution”

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 505 records · Page 28

Alleviating prior dependencies for DESI DR1 clustering fits through reparameterization

Bayesian analyses of the full-shape clustering of Dark Energy Spectroscopic Instrument (DESI) Data Release 1 (DR1) exhibit prior-volume projection effects, whereby weakly constrained nuisance parameters of the Effective Field Theory of Large Scale Structure (EFTofLSS) shift marginalized cosmological posteriors away from the posterior maximum. We reanalyze DESI DR1 power spectrum multipoles using two complementary mitigation strategies: (i) nonlinear orthogonalization to decorrelate nuisance and cosmological parameter priors, and (ii) a fully reparameterization-invariant Jeffreys prior over all EFTofLSS coefficients, evaluated on-the-fly via closed-form Jacobians. Including data from DESI, Big-Bang Nuclesynthesis and a constraint on $n_{\mathrm{s}}$, baseline priors lead to multi-$σ$ projection in the Hubble parameter $H_{0}$ and dark energy equation of state parameters $w_{0}$ and $w_{a}$; the Jeffreys prior successfully recenters these posteriors to enclose the maximum a posteriori estimate within the 68% credible regions, demonstrating clear mitigation of projection effects for these late-time expansion parameters. A hybrid Jeffreys+baseline-Gaussian configuration controls residual over-broad tails in the physical cold dark matter density $ω_{\mathrm{c}}$ while preserving the volume correction, and is our favoured approach. We compare the credible intervals derived using our methodology to those obtained using Halo Occupation Distribution (HOD)-informed priors and to confidence intervals derived using frequentist profile likelihood analyses, finding agreement in both central values and degeneracy directions in the $w_{0}$--$w_{a}$ plane. This demonstrates that, once projection effects are properly controlled, we can make robust inferences about the late-time cosmological expansion independent of the statistical framework adopted.

Bonici, M. [Waterloo U.; Perimeter Inst. Theor. Ph↗

A Bimodal Diagnostic Cloud Fraction Parameterization. Part I: Motivating Analysis and Scheme Description

Cloud fraction parameterizations are beneficial to regional, convection-permitting numerical weather prediction. For its operational regional midlatitude forecasts, the Met Office uses a diagnostic cloud fraction scheme that relies on a unimodal, symmetric subgrid saturation-departure distribution. This scheme has been shown before to underestimate cloud cover and hence an empirically based bias correction is used operationally to improve performance. This first of a series of two papers proposes a new diagnostic cloud scheme as a more physically based alternative to the operational bias correction. The new cloud scheme identifies entrainment zones associated with strong temperature inversions. Additionally, for model grid boxes located in this entrainment zone, collocated moist and dry Gaussian modes are used to represent the subgrid conditions. The mean and width of the Gaussian modes, inferred from the turbulent characteristics, are then used to diagnose cloud water content and cloud fraction. It is shown that the new scheme diagnoses enhanced cloud cover for a given gridbox mean humidity, similar to the current operational approach. It does so, however, in a physically meaningful way. Using observed aircraft data and ground-based retrievals over the southern Great Plains in the United States, it is shown that the new scheme improves the relation between cloud fraction, relative humidity, and liquid water content. An emergent property of the scheme is its ability to infer skewed and bimodal distributions from the large-scale state that qualitatively compare well against observations. A detailed evaluation and resolution sensitivity study will follow in Part II.

54 ENVIRONMENTAL SCIENCES↗

Search for new phenomena in multi-body invariant masses in events with at least one isolated lepton and two jets using $\sqrt{s}$ = 13 TeV proton–proton collision data collected by the ATLAS detector

A search for resonances in events with at least one isolated lepton (e or μ ) and two jets is performed using 139 fb -1 of $\sqrt{s}$ = 13 TeV proton–proton collision data recorded by the ATLAS detector at the LHC. Deviations from a smoothly falling background hypothesis are tested in three- and four-body invariant mass distributions constructed from leptons and jets, including jets identified as originating from bottom quarks. Model-independent limits on generic resonances characterised by cascade decays of particles leading to multiple jets and leptons in the final state are presented. The limits are calculated using Gaussian shapes with different widths for the invariant masses. The multi-body invariant masses are also used to set 95% confidence level upper limits on the cross-section times branching ratios for the production and subsequent decay of resonances predicted by several new physics scenarios.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Cosmology with persistent homology: a Fisher forecast

Abstract Persistent homology naturally addresses the multi-scale topological characteristics of the large-scale structure as a distribution of clusters, loops, and voids. We apply this tool to the dark matter halo catalogs from theQuijotesimulations, and build a summary statistic for comparison with the joint power spectrum and bispectrum statistic regarding their information content on cosmological parameters and primordial non-Gaussianity. Through a Fisher analysis, we find that constraints from persistent homology are tighter for 8 out of the 10 parameters by margins of 13–50%. The complementarity of the two statistics breaks parameter degeneracies, allowing for a further gain in constraining power when combined. We run a series of consistency checks to consolidate our results, and conclude that our findings motivate incorporating persistent homology into inference pipelines for cosmological survey data.

Astronomy & Astrophysics↗

Two-pion femtoscopic correlations in Be+Be collisions at $\sqrt{s_{\text {NN}}} = 16.84$ GeV measured by the NA61/SHINE at CERN

This paper reports measurements of two-pion femtoscopic correlations in Be+Be collisions at a beam momentum of 150$A\,\hbox {GeV}\!/\!c$ (energy available in the center-of-mass system for nucleon pair $\sqrt{s_{\text {NN}}} = 16.84$ GeV) by the NA61/SHINE experiment at the CERN SPS accelerator. The obtained momentum space correlation functions can be well described by a Lévy distributed source model. The transverse mass dependence of the Lévy source parameters is presented, and their possible theoretical interpretations are discussed. The results show that the Lévy exponent $\alpha $ is approximately constant as a function of $m_{\text {T}}$ , and far from both the Gaussian case of $\alpha = 2$ or the conjectured value at the critical endpoint, $\alpha = 0.5$. The radius scale parameter R shows a slight decrease in $m_{\text {T}}$ , which can be explained as a signature of transverse flow. Finally, an approximately constant trend of the intercept parameter $\lambda $ as a function of $m_{\text {T}}$ was observed, similar to previous NA44 S + Pb results (obtained with a Gaussian approximation, but unlike RHIC results).

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Improving qubit readout with hidden Markov models

We demonstrate the application of pattern recognition algorithms via hidden Markov models (HMM) for qubit readout. This scheme provides a state-path trajectory approach capable of detecting qubit-state transitions and makes for a robust classification scheme with higher starting-state assignment fidelity than when compared to a multivariate Gaussian or a support vector machine scheme. Therefore, the method also eliminates the qubit-dependent readout time optimization requirement in current schemes. Using a HMM state discriminator we estimate fidelities reaching the ideal limit. Unsupervised learning gives access to transition matrix, priors, and IQ distributions, providing a toolbox for studying qubit-state dynamics during strong projective readout.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Maximizing the impact of the Vera Rubin Observatory LSST, and DESI and CMB-S4 experiments in constraining gravity and dark sector physics

The work in this award has addressed our objectives on a number of fronts: Current graduate student Chris Wilson and Bean are determining how cosmic voids can be used to test the properties of gravity, leveraging simulations from DESI collaborator Baojiu Li and refining theoretical Gaussian Streaming Model approaches to predict galaxy dynamics in voids when the cosmology is not GR. This work is still underway, including analyses of the efficacy of void quadrupole correlation statistics and void size functions when studying dark matter halo statistics or the expected galaxy distributions from upcoming galaxy surveys; Current graduate student Yulin Gong, Bean with collaborators are refining techniques to extract out the kinetic Sunyaev Zel’dovich effect from galaxy cluster survey data and CMB temperature anisotropies, with the intent to further ready this pipeline development in readiness for spectroscopic data from DESI and LSST along with CMB-S4 measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

The Impacts of Single-Scattering and Microphysical Properties of Ice Particles Smaller Than 100 µm on the Bulk Radiative Properties of Tropical Cirrus

There are large uncertainties in the single-scattering (i.e., morphologies) and microphysical (i.e., concentrations) properties of ice particles whose size are less than ~100 µm. Insufficient resolutions of the most advanced cloud probes (e.g., cloud particle imager) cannot resolve the micrometer-scale morphologies of small ice particles. Further, the shattering of large ice particles on probes’ inlets or tips causes uncertainties in the measurement of the concentrations of small ice particles. These uncertainties have large impacts on the single-scattering and microphysical properties of small ice particles that are utilized to quantify the bulk radiative properties of cirrus. In this study, the impacts of uncertainties in the morphologies and concentrations of small ice particles on the bulk radiative properties of tropical cirrus were calculated using measurements acquired during the Tropical Warm Pool-International Cloud Experiment. Five different models (i.e., budding Buckyball, Chebyshev particle, droxtal, Gaussian random sphere, and sphere) that represent the shapes of small ice particles were used to calculate the single-scattering properties. The bulk radiative properties, average phase-function ($\overline{P_{11}}$), and average asymmetry parameter ($\overline{g}$) were computed by combining the measured size/habit distributions and the calculated single-scattering properties of ice particles. The impacts of the selection of varying morphologies of small particles on the bulk radiative properties were quantified. For these calculations, the possible range of the concentrations of small ice particles which depend on the degree of shattered large particles were also used. The impacts of varying the single-scattering properties of small ice particles on the bulk radiative properties were the largest in the upper parts of cirrus (T < –60 °C), while they were the smallest in the lower parts of cirrus (–45 < T < –30 °C). The impacts of uncertainties in the concentrations of small ice particles on the bulk radiative properties were largest in the lower parts of cirrus (–45 < T < –30 °C), whereas they were smallest in the upper parts of cirrus (T < –60 °C). The effect of shattering was maximum in the lower parts of cirrus, whilst it was minimum in the upper parts of cirrus. The combined impacts of uncertainties in the single-scattering (i.e., morphologies) and microphysical (i.e., concentrations) properties of small ice particles revealed variations of up to 11.2% (127.1%; 67.3%) of the integrated intensity in the forward (sideward; backward) angles in $\overline{P_{11}}$ and a corresponding change in $\overline{g}$ by up to 12.61%.

54 ENVIRONMENTAL SCIENCES↗

Tomography of the rho meson in the QCD instanton vacuum: Transverse momentum dependent parton distribution functions

We analyze the rho meson unpolarized and polarized transverse momentum dependent parton distribution functions (TMDPDFs) in the instanton liquid model (ILM). The corresponding TMDs in ILM are approximated by a constituent quark beam function in the leading Fock state multiplied by a rapidity-dependent factor resulting from the staple-shaped Wilson lines, for fixed longitudinal momentum, transverse separation, and rapidity. At the resolution of the ILM, all of the rho meson TMDs are symmetric in parton 𝑥 for fixed transverse momentum, and Gaussian-like in the transverse momentum for fixed parton 𝑥. The latter is a direct consequence of the profiling of the quark zero modes in the ILM. The evolved TMDs at higher rapidity using the Collins-Soper kernel, and higher resolution using the renormalization group, show substantial skewness towards low parton 𝑥.

Gluons↗

Transverse momentum structure of proton within the basis light-front quantization framework

We obtain the leading-twist valence quark transverse-momentum-dependent parton distribution functions (TMD PDFs) for the proton within the basis light-front quantization (BLFQ) framework. Our results are consistent with lattice QCD calculations and our previous results for the collinear limit. We also obtain consistency with the Soffer-type bounds. Within our approach, we find that six T-even TMDs in the leading twist are all independent of each other, and previously found model-dependent relations do not hold. This is a promising sign that our results are representative of future, more extensive treatments of QCD. Furthermore, we obtain a non-trivial x-dependence of the <(p ⊥ ) 2 > and some consistency with the Gaussian ansatz but only in the small <(p ⊥ ) 2 > region. Those features suggest our results may be a useful alternative in future experimental extractions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Scalable statistical inference of photometric redshift via data subsampling

Handling big data has largely been a major bottleneck in traditional statistical models. Consequently, when accurate point prediction is the primary target, machine learning models are often preferred over their statistical counterparts for bigger problems. But full probabilistic statistical models often outperform other models in quantifying uncertainties associated with model predictions. We develop a data-driven statistical modeling framework that combines the uncertainties from an ensemble of statistical models learned on smaller subsets of data carefully chosen to account for imbalances in the input space. We demonstrate this method on a photometric redshift estimation problem in cosmology, which seeks to infer a distribution of the redshift—the stretching effect in observing the light of far-away galaxies—given multivariate color information observed for an object in the sky. Our proposed method performs balanced partitioning, graph-based data subsampling across the partitions, and training of an ensemble of Gaussian process models.

data subsampling↗

Multiple emission components in the Cygnus cocoon detected from Fermi -LAT observations

Star-forming regions may play an important role in the life cycle of Galactic cosmic rays (CRs), notably as home to specific acceleration mechanisms and transport conditions. Gamma-ray observations of Cygnus X have revealed the presence of an excess of hard-spectrum gamma-ray emission, possibly related to a cocoon of freshly accelerated particles. We seek an improved description of the gamma-ray emission from the cocoon using ~13 yr of observations with the Fermi-Large Area Telescope (LAT) and use it to further constrain the processes and objects responsible for the young CR population. We developed an emission model for a large region of interest, including a description of interstellar emission from the background population of CRs and recent models for other gamma-ray sources in the field. Thus, we performed an improved spectro-morphological characterisation of the residual emission including the cocoon. The best-fit model for the cocoon includes two main emission components: an extended component FCES G78.74+1.56, described by a 2D Gaussian of extension r 68 = 4.4° ± 0.1° -0.1° +0.1° and a smooth broken power law spectrum with spectral indices 1.67 ± 0.05 -0.01 +0.02 and 2.12 ± 0.02 -0.01 +0.00 below and above 3.0 ± 0.6 -0.2 +0.0 GeV, respectively; and a central component FCES G80.00+0.50, traced by the distribution of ionised gas within the borders of the photo-dissociation regions and with a power law spectrum of index 2.19 ± 0.03 -0.01 +0.00 that is significantly different from the spectrum of FCES G78.74+1.56. An additional extended emission component FCES G78.83+3.57, located on the edge of the central cavities in Cygnus X and with a spectrum compatible with that of FCES G80.00+0.50, is likely related to the cocoon. For the two brightest components FCES G80.00+0.50 and FCES G78.74+1.56, spectra and radial-azimuthal profiles of the emission can be accounted for in a diffusion-loss framework involving one single population of non-thermal particles with a flat injection spectrum. Particles span the full extent of FCES G78.74+1.56 as a result of diffusion from a central source, and give rise to source FCES G80.00+0.50 by interacting with ionised gas in the innermost region. For this simple diffusion-loss model, viable setups can be very different in terms of energetics, transport conditions, and timescales involved, and both hadronic and leptonic scenarios are possible. The solutions range from long-lasting particle acceleration, possibly in prominent star clusters such as Cyg OB2 and NGC 6910, to a more recent and short-lived release of particles within the last 10–100 kyr, likely from a supernova remnant. The observables extracted from our analysis can be used to perform detailed comparisons with advanced models of particle acceleration and transport in star-forming regions.

79 ASTRONOMY AND ASTROPHYSICS↗

An experimental and computational analysis of combustion heat release transformation in dual fuel combustion

Dual fuel (DF) diesel-methane combustion, which employs a high-reactivity fuel (diesel) to ignite a low-reactivity fuel (methane), is a widely studied combustion strategy for internal combustion engines, with significant potential for engine-out emissions reductions without the need for major hardware modifications. A phenomenon, which has been reported in the DF literature, but not explained fully, is the transformation of the shape of the apparent heat release rate (AHRR) curve as the start of injection (SOI) of diesel is advanced beyond a certain threshold; coincidentally, this AHRR transformation is usually accompanied by a sharp decrease in engine-out emissions of oxides of nitrogen (NOx). The goal of the present work is to establish the underlying physical reason(s) that cause the AHRR transformation. The AHRR transformation was observed on a single cylinder research engine (SCRE) at an indicated mean effective pressure (IMEP) of 5 bar at a speed of 1500 rev/min. The transformation occurred over a range of SOIs from 330 to 320 crank angle degrees (CAD). While the 330 CAD SOI exhibited a typical two-stage AHRR curve, with a clearly definable first-stage peak followed by a second-stage AHRR with little-to-no low temperature heat release (LTHR) present and high engine-out NOx, the 320 CAD SOI exhibited a single-stage, Gaussian-like AHRR curve, with noticeable LTHR and at least one order-of-magnitude lower NOx emissions. Here, leveraging analysis of experimental data and three-dimensional computational fluid dynamic simulations, the authors show that the AHRR transformation is impacted mainly by differences in local equivalence ratio distributions within the cylinder at ignition onset for different diesel SOIs.

33 ADVANCED PROPULSION SYSTEMS↗

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↗

Atomistic Simulations of Polydisperse Lignin Melts Using Simple Polydisperse Residue Input Generator

Understanding the physics of lignin will help rationalize its function in plant cell walls as well as aiding practical applications such as deriving biofuels and bioproducts. Here, in this work, we present SPRIG (Simple Polydisperse Residue Input Generator), a program for generating atomic-detail models of random polydisperse lignin copolymer melts i.e., the state most commonly found in nature. Using these models, we use all-atom molecular dynamics (MD) simulations to investigate the conformational and dynamic properties of polydisperse melts representative of switchgrass (Panicum virgatum L.) lignin. Polydispersity, branching and monolignol sequence are found to not affect the calculated glass transition temperature, T g . The Flory–Huggins scaling parameter for the segmental radius of gyration is 0.42 ± 0.02, indicating that the chains exhibit statistics that lie between a globular chain and an ideal Gaussian chain. Below T g the atomic mean squared displacements are independent of molecular weight. In contrast, above T g , they decrease with increasing molecular weight. Therefore, a monodisperse lignin melt is a good approximation to this polydisperse lignin when only static properties are probed, whereas the molecular weight distribution needs to be considered while analyzing lignin dynamics.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Directivity Modes of Earthquake Populations with Unsupervised Learning

We present a novel approach for resolving modes of rupture directivity in large populations of earthquakes. A seismic spectral decomposition technique is used to first produce relative measurements of radiated energy for earthquakes in a spatially compact cluster. The azimuthal distribution of energy for each earthquake is then assumed to result from one of several distinct modes of rupture propagation. Rather than fitting a kinematic rupture model to determine the most likely mode of rupture propagation, we instead treat the modes as latent variables and learn them with a Gaussian mixture model. The mixture model simultaneously determines the number of events that best identify with each mode. The technique is demonstrated on four datasets in California, each with compact clusters of several thousand earthquakes with comparable slip mechanisms. We show that the datasets naturally decompose into distinct rupture propagation modes that correspond to different rupture directions, and the fault plane is unambiguously identified for all cases. We find that these small earthquakes exhibit unilateral ruptures 63–73% of the time on average. Here, the results provide important observational constraints on the physics of earthquakes and faults.

58 GEOSCIENCES↗

A Data-Driven Nonparametric Approach for Probabilistic Load-Margin Assessment Considering Wind Power Penetration

A modern power system is characterized by an increasing penetration of wind power, which results in large uncertainties in its states. These uncertainties must be quantified properly; otherwise, the system security may be threatened. Facing this challenge, here we propose a cost-effective, data-driven approach to assessing a power system's load margin probabilistically. Using actual wind data, a kernel density estimator is applied to infer the nonparametric wind speed distributions, which are further merged into the framework of a vine copula. The latter enables us to simulate complex multivariate and highly dependent model inputs with a variety of bivariate copulae that precisely represent the tail dependence in the correlated samples. Furthermore, to reduce the prohibitive computational time of traditional Monte-Carlo simulations that process a large amount of samples, we propose to use a nonparametric, Gaussian-process-emulator-based reduced-order model to replace the original complicated continuation power-flow model through a Bayesian-learning framework. To accelerate the convergence rate of this Bayesian algorithm, a truncated polynomial chaos surrogate, which serves as a highly efficient, parametric Bayesian prior, is developed. This emulator allows us to execute the time-consuming continuation power-flow solver at the sampled values with a negligible computational cost. Results of simulations that are performed on several test systems reveal the impressive performance of the proposed method in the probabilistic load-margin assessment.

17 WIND ENERGY↗

Optimization of the number and locations of the calibration stations needed to monitor soil moisture using distributed temperature sensing systems: A proof-of-concept study

The single-probe heat-pulse (SPHP) technique combined with the Fiber-optic Distributed Temperature Sensing (DTS) technology can offer novel high-resolution measurements of soil moisture (θ) over spatial scales ranging from several centimeters to several kilometers. However, the key limitation of this method is in obtaining the calibration relationship between θ and soil thermal conductivity (λ) across a specific field. In a previous study, a new methodology using a Gaussian processes model was presented to account for the spatial variability in the λ-θ relationship. The model aggregated θ measurements from soil moisture sensors scattered over the SPHP transect with the corresponding DTS λ measurements at their locations. In this study, a novel methodology is tested to optimize the number and locations of soil moisture sensors required to account for the spatial variability of the λ - θ relationship to achieve higher accuracy from the SPHP technique. The proposed methodology utilizes hierarchical clustering to analyze the information contained in the spatial structure of the SPHP measurements as the soil dries from a nearly-saturated condition. The proposed methodology was tested using data from a field in Oklahoma. Monte-Carlo simulation was performed to validate the performance of the proposed methodology. The predictions obtained from the proposed methodology resulted in θ measurements accuracy comparable to those obtained from the 10% best Monte-Carlo iterations of randomly assigned soil moisture locations. Further, this study demonstrates that the proposed methodology is more efficient than the traditional practice of randomly spreading calibration soil moisture sensors along the SPHP transect.

54 ENVIRONMENTAL SCIENCES↗