Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “ordered statistics”

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,243 records · Page 69

Building a 2.5D Digital Elevation Model from 2D Imagery

When projecting imagery into a georeferenced coordinate frame, one needs to have some model of the geographical region that is being projected to. This model can sometimes be a simple geometrical curve, such as an ellipse or even a plane. However, to obtain accurate projections, one needs to have a more sophisticated model that encodes the undulations in the terrain including things like mountains, valleys, and even manmade structures. The product that is often used for this purpose is a Digital Elevation Model (DEM). The technology presented here generates a high-quality DEM from a collection of 2D images taken from multiple viewpoints, plus pose data for each of the images and a camera model for the sensor. The technology assumes that the images are all of the same region of the environment. The pose data for each image is used as an initial estimate of the geometric relationship between the images, but the pose data is often noisy and not of sufficient quality to build a high-quality DEM. Therefore, the source imagery is passed through a feature-tracking algorithm and multi-plane-homography algorithm, which refine the geometric transforms between images. The images and their refined poses are then passed to a stereo algorithm, which generates dense 3D data for each image in the sequence. The 3D data from each image is then placed into a consistent coordinate frame and passed to a routine that divides the coordinate frame into a number of cells. The 3D points that fall into each cell are collected, and basic statistics are applied to determine the elevation of that cell. The result of this step is a DEM that is in an arbitrary coordinate frame. This DEM is then filtered and smoothed in order to remove small artifacts. The final step in the algorithm is to take the initial DEM and rotate and translate it to be in the world coordinate frame [such as UTM (Universal Transverse Mercator), MGRS (Military Grid Reference System), or geodetic] such that it can be saved in a standard DEM format and used for projection.

Padgett, Curtis W.↗

High-Resolution Mesoscale Model Setup for the Eastern Range and Wallops Flight Facility

Mesoscale weather conditions can have an adverse effect on space launch, landing, ground processing, and weather advisories, watches, and warnings at the Eastern Range (ER) in Florida and Wallops Flight Facility (WFF) in Virginia. During summer, land-sea interactions across Kennedy Space Center (KSC) and Cape Canaveral Air Force Station (CCAFS) lead to sea breeze front formation, which can spawn deep convection that can hinder operations and endanger personnel and resources. Many other weak locally-driven low-level boundaries and their interactions with the sea breeze front and each other can also initiate deep convection in the KSC/CCAFS area. These convective processes often last 60 minutes or less and pose a significant challenge to the local forecasters. Surface winds during the transition seasons (spring and fall) pose the most difficulties for the forecasters at WFF. They also encounter problems forecasting convective activity and temperature during those seasons. Therefore, accurate mesoscale model forecasts are needed to better forecast a variety of unique weather phenomena. Global and national scale models cannot properly resolve important local-scale weather features at each location due to their horizontal resolutions being much too coarse. Therefore, a properly tuned local data assimilation (DA) and forecast model at a high resolution is needed to provide improved capability. To accomplish this, a number of sensitivity tests were performed using the Weather Research and Forecasting (WRF) model in order to determine the best DA/model configuration for operational use at each of the space launch ranges to best predict winds, precipitation, and temperature. A set of Perl scripts to run the Gridpoint Statistical Interpolation (GSI)/WRF in real-time were provided by NASA's Short-term Prediction Research and Transition Center (SPoRT). The GSI can analyze many types of observational data including satellite, radar, and conventional data. The GSI/WRF scripts use a cycled GSI system similar to the operational North American Mesoscale (NAM) model. The scripts run a 12-hour pre-cycle in which data are assimilated from 12 hours prior up to the model initialization time. A number of different model configurations were tested for both the ER and WFF by varying the horizontal resolution on which the data assimilation was done. Three different grid configurations were run for the ER and two configurations were run for WFF for archive cases from 27 Aug 2013 through 10 Nov 2013. To quantify model performance, standard model output will be compared to the Meteorological Assimilation Data Ingest System (MADIS) data. The MADIS observation data will be compared to the WRF forecasts using the Model Evaluation Tools (MET) verification package. In addition, the National Centers for Environmental Prediction's Stage IV precipitation data will be used to validate the WRF precipitation forecasts. The author will summarize the relative skill of the various WRF configurations and how each configuration behaves relative to the others, as well as determine the best model configuration for each space launch range.

Watson, Leela R.↗

Quantitative Approach for the Accurate CFD Simulation of Hover in Turbulent Flow

Time-dependent Navier-Stokes simulations have been carried out for a V22 rotor in hover using an improved, low-diffusion, HLLE++ upwind algorithm in the OVERFLOW code. Emphasis is placed on lessons learned over the past decade regarding the effects of high-order spatial accuracy, grid resolution, and the use of detached eddy simulation in predicting the rotor performance, i.e., figure-of-merit (FM). A quick-start procedure is described together with a statistical measure of FM convergence that reduces hover computations by 5 times, similar in computational work for forward flight. Furthermore, Cartesian adaptive mesh refinement in the rotor wake revealed a complex turbulent flow with worm-like structures of various scales. These turbulent worms, found computationally more than a decade ago, have recently been observed in a separate German Aerospace Center (DLR) experiment. Moreover, adaptive mesh refinement has been used to resolve the tip-vortex to its correct physical size.

Neal M. Chaderjian↗

Dark Energy Survey Year 3 results: Cosmology from combined galaxy clustering and lensing validation on cosmological simulations

In this work, we present a validation of the Dark Energy Survey Year 3 (DES Y3) 3×2-point analysis choices by testing them on Buzzard2.0, a new suite of cosmological simulations that is tailored for the testing and validation of combined galaxy clustering and weak-lensing analyses. We show that the buzzard2.0 simulations accurately reproduce many important aspects of the DES Y3 data, including photometric redshift and magnitude distributions, and the relevant set of two-point clustering and weak-lensing statistics. We then show that our model for the 3×2-point data vector is accurate enough to recover the true cosmology in simulated surveys assuming the true redshift distributions for our source and lens samples, demonstrating robustness to uncertainties in the modeling of the nonlinear matter power spectrum, nonlinear galaxy bias, and higher-order lensing corrections. Additionally, we demonstrate for the first time that our photometric redshift calibration methodology, including information from photometry, spectroscopy, clustering cross-correlations, and galaxy–galaxy lensing ratios, is accurate enough to recover the true cosmology in simulated surveys in the presence of realistic photometric redshift uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Hemispherical power asymmetry in intensity and polarization for Planck PR4 data

Abstract One of the foundations of the Standard Model of Cosmology is statistical isotropy, which can be tested, among other probes, through the study of the Cosmic Microwave Background (CMB). However, a hemispherical power asymmetry on large scales has been reported for WMAP and Planck data by different works.The statistical significance is above 3σfor temperature, suggesting a directional dependence of the local power spectrum, and thus a feature beyond the ΛCDM model. With the third release of the Planck data (PR3), a new analysis was performed including the E-mode polarization maps, finding an asymmetry at a modest level of significance. In this work, we perform an asymmetry analysis in intensity and polarization maps for the latest Planck processing pipeline (PR4). We obtain similar results to those obtained with PR3, with a slightly lower significance (2.8% for the Sevem method)for the amplitude of the E-mode local variance dipole as well as a significant variability with the considered mask.In addition, a hint of a possible T-E alignment between the asymmetry axes is found at the level of ∼ 5%. For the analysis, we have implemented an alternative inpainting approach in order to get an accurate reconstruction of the E-modes. More sensitive all-sky CMB polarization data, such as those expected from the future LiteBIRD experiment, are needed to reach a more robust conclusion on the possible existence of deviations from statistical isotropy in the form of a hemispherical power asymmetry.

Astronomy & Astrophysics↗

Statistical Short-Range Guidance for Peak Wind Speed Forecasts at Edwards Air Force Base, CA

The peak winds near the surface are an important forecast element for space shuttle landings. As defined in the Flight Rules (FR), there are peak wind thresholds that cannot be exceeded in order to ensure the safety of the shuttle during landing operations. The National Weather Service Spaceflight Meteorology Group (SMG) is responsible for weather forecasts for all shuttle landings, and is required to issue surface average and 10-minute peak wind speed forecasts. They indicate peak winds are a challenging parameter to forecast. To alleviate the difficulty in making such wind forecasts, the Applied Meteorology Unit (AMU) developed a PC-based graphical user interface (GUI) for displaying peak wind climatology and probabilities of exceeding peak wind thresholds for the Shuttle Landing Facility (SLF) at Kennedy Space Center (KSC; Lambert 2003). However, the shuttle occasionally may land at Edwards Air Force Base (EAFB) in southern California when weather conditions at KSC in Florida are not acceptable, so SMG forecasters requested a similar tool be developed for EAFB.

Dreher, Joseph G.↗

Proposed Mars Surveyor Landing Sites in Northern Meridiani Sinus, Southern Elysium Planitia, and Argyre Planitia

Our objective is to propose two landing sites that the Mars Surveyor 2001 Lander and Athena Rover could go to on Mars that should meet the safety requirements of the spacecraft landing system and optimize surface operations (chiefly driven by power and communications requirements). An additional site within Argyre Planitia, initially proposed by Parker to the Mars Surveyor Landing Site program, is also proposed for potential consideration for post-2001 missions to Mars, as it is well outside the current latitude limits for the Athena Rover. All three sites are designed to be situated as close to a diversity of geologic units within a few kilometers of the landing site so that diversity can be placed in a geologic context. This objective is very different from the Mars Pathfinder requirement to land at a site with a maximum chance for containing a diversity of rocks within a few tens of meters of the lander. That requirement was driven by the Sojourner mobility limit of a few tens of meters. It can be argued that the Athena project, with its much larger mobility capability, might actually want to avoid such a site, because placing collected samples in geologic context would be difficult. While it has been argued, both before and after the Mars Pathfinder landing, that the provenance for local blocks may be determined by orbiter spectra, primarily from the MGS TES instrument, our ability to do so has yet to be demonstrated. Indeed, several months after conclusion of the Pathfinder mission, we have yet to reach a consensus on the composition of local materials. Our primary data set for selecting a landing site within the latitude and elevation constraints of the 2001 mission is the Viking Orbiter image archive. The site must be selected to place the landing ellipse so as to avoid obvious hazards, such as steep slopes, large or numerous craters, or abundant large knobs. For this purpose, we chose a resolution limit of better than 50 m/pixel. This necessarily excludes from the present study images from current and future orbiter spacecraft, until such data does become readily available. Within each proposed region, it may be possible to identify additional sites once these data become available. Second, the fine-component thermal inertia data, should be greater than about 5 or 6 cgs Units (10(exp -3) cal/sq cm s(exp -0.5)/K). Low thermal inertias imply dusty environments, which could pose a mobility hazard. Similarly, the albedo of the site should not be particularly high, which would also suggest dusty surfaces. Low albedos are preferred, as they often coincide with low Viking red:violet ratios and indicate less dusty surfaces. Next, the Modeled Block Abundance should also not be too high or too low. Based on the Viking Lander and Mars Pathfinder experiences, percentages of blocks should be on the order of 5-25%. Too many blocks could pose a hazard to the landing and mobility. Too few blocks could also indicate a dusty surface. Primary Landing Site: Northern Meridiani Sinus (Proposed by T. J. Parker and K., S. Edgett) Vital Statistics: (1) Latitude, Longitude: 0-3 N, 350-2 W. *Elevation (Viking): about0.5-1.5 Ian. (2) Viking Orbiter Image coverage: Excellent coverage by 15 - 25 m/pixel images (orbits 709A and 410B). Possible stereo coverage in region where two orbits overlap (probably small parallax angle, as these orbits are not listed in NASA Contractor Report 3501) (3) Albedo: about .18 -.26 (4) Block Abundance: 5-26% (5)Fine-Component Thermal Inertia: 5-9 cgs units This region consists of bright deposits similar to those described by Edgett et al, that also lie within a prominent dark albedo region. These deposits are flat-lying, to such a degree that they ramp against topography rather than draping over it. This led Edgett and Parker to suggest that they may be subaqueous sediments, possibly lacustrine or marine evaporites, laid down sometime from the late Noachian to middle Hesperian (age determination pending crater counts). A contact between this material and elevated, dissected highlands to the south was identified , and is described by Edgett et al. Our desire in proposing this landing site is to sample the edge of this deposit where it has been exposed through etching, presumably eolian deflation (the deposit, though in the highlands, is itself only lightly to moderately cratered). This should enable access to in situ stratigraphy. The actual landing site will be selected where slopes are not expected to be steep, such that the rover itself should be able to traverse them and sample layered materials on the way, either up or down the slope. Perhaps due to uncertainties at this time as to the friability or meter-scale roughness of the deposit, it might make sense to place the landing ellipse on the exhumed highland surface adjacent to the deflated margin of the deposit and plan on driving to the deposit rather than landing on it and driving downslope. This should also enable imaging the margin for evidence of layering should it prove too difficult to climb. A target ellipse on the highland surface should also allow Athena access to ancient Noachian highland materials, particularly if placed near crater ejecta or an inlier of knobby material. Secondary Landing Site: Southern Elysium Planitia (Proposed by T. J. Parker) Vital Statistics: (1) Latitude, Longitude: 1.5-3.5 S, 195-198 W. (2) Elevation (Viking): -1.0 km. (3) Viking Orbiter Image coverage: Excellent coverage by 15 - 25 m/pixel images (orbit 725). Possible stereo coverage between images from beginning and end of orbit that overlap (probably small parallax angle) (4) Albedo: about .27-.28 (5) Block Abundance: 4-7% (6) Fine-Component Thermal Inertia: about 3 cgs units This region consists of eroded knobby material, probably of Noachian age, though much of the crater population has been destroyed, that is onlapped at a sharp contact by an extensive plains unit in southern Elysium Planitia that is Amazonian in age. The plains materials have been attributed to unusually low-viscosity flood lavas from fissures south of the Elysium volcanic rise, or to lacustrine materials associated with a large, Amazonian lake at the source of Marte Vallis. Parker and Schenk presented evidence in support of the latter interpretation, though they attributed the putative shore morphology to an embayment of a northern plains ocean into the southern Elysium region. Detailed examination of the margin of the deposit, showing erosion, not simply burial, of small crater rims and fluidized ejecta blankets, also points to lacustrine or marine sedimentation rather than volcanic plains burial. The plains surface exhibits a "crusty" appearance that many researchers have attributed to pressure ridges in lava flows. In a lacustrine context, they also resemble pressure ridges in desiccated evaporite deposits and salt-rimmed pools (now dry) similar in scale and morphology to spectacular, hundred meter-scale pool rims in alkaline Lake Natron, East African Rift. The eroded highland margin surface adjacent to these plains appears to be fairly smooth, even at 15 m/pixel. Isolated knob inliers are scattered from a few kilometers to several tens of "kilometers apart. Heights of the knobs have not been measured yet but, based on experience with similar features in the Pathfinder landing ellipse, are probably typically on the order of several tens of meters high and smaller, though some of the largest knobs in the region are probably up to a few hundred meters high. Two craters larger than a kilometer in diameter, with fluidized deposits, lie nearby the proposed landing site. Very high-resolution images from MOC should help to determine whether a landing site navigable by the Athena rover could be placed in this region. The space between knobs and craters is large enough to enable placement of a target landing ellipse between them but still provide access to one or more of them and to the margin of the Elysium plains material. Post-2001 Mars Surveyor Landing Site: Argyre Planitia (Proposed by T. J. Parker) Vital Statistics: (1) Latitude, Longitude: 55-56 S, 41-43 W. (2) Elevation (Viking): 1.0 km. (3) Viking Orbiter Image coverage: Excellent coverage by 40 m/pixel images (orbits 567B, 568B, and 569B). Excellent stereo coverage with large parallax angles over the entire landing site region, and much of central and southern Argyre. (4) Albedo: about .23-.24 (5) Block Abundance: No data (6) Fine-Component Thermal Inertia: No data The floors of both the Argyre and Hellas basins contain etched layered materials that are probably thick accumulations of channel or lacustrine sediments. The deposits in Hellas are much more eroded than those in Argyre, and Hellas lacks a channel outlet. Argyre is unique in that Uzboi Vallis flowed out of the basin, requiring overflow of a standing body of water within Argyre. This makes it the largest impact basin on Mars with channels both draining into it and flowing out from it. Hellas' channels may be catastrophic flood channels, whereas Argyre was fed by modest-scale valley networks, though the outlet at Uzboi Vallis was a catastrophic flood Highland craters and basins of this kind should be high-priority landing targets for missions intended to focus on the search for either prebiotic organic materials or even simple fossil microorganisms. Basins with internally-draining valley networks should be preferred over flood channels, as they could have provided the long-term influx of water favorable to the origin of life. (Catastrophic floods are not conducive to fossil preservation, due to their very short durations and high transportation energies). They also afford an opportunity to study the evolution of the planet's climate and volatiles during the period of time between the late Noachian and early Hesperian, when a drastic change from a proposed early warm, wet climate to one more closely resembling the modern environment is thought to have occurred. Large basin

Parker, T. J.↗

Linear growth of circuit complexity from Brownian dynamics

How rapidly can a many-body quantum system generate randomness? Using path integral methods, we demonstrate that Brownian quantum systems have circuit complexity that grows linearly with time. In particular, we study Brownian clusters of N spins or fermions with time-dependent all-to-all interactions, and calculate the Frame Potential to characterize complexity growth in these models. In both cases the problem can be mapped to an effective statistical mechanics problem which we study using path integral methods. Within this framework it is straightforward to show that the kth Frame Potential comes within ϵ of the Haar value after a time of order t ~ kN + k log k + log ϵ –1 . Using a bound on the diamond norm, this implies that such circuits are capable of coming very close to a unitary k-design after a time of order t ~ kN. We also consider the same question for systems with a time-independent Hamiltonian and argue that a small amount of time-dependent randomness is sufficient to generate a k-design in linear time provided the underlying Hamiltonian is quantum chaotic. These models provide explicit examples of linear complexity growth that are analytically tractable and are directly applicable to practical applications calling for unitary k-designs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Robust verification of stochastic simulation codes

We introduce a robust verification tool for computational codes, which we call Stochastic Robust Extrapolation based Error Quantification (StREEQ). Unlike the prevalent Grid Convergence Index (GCI) [1] method, our approach is suitable for both stochastic and deterministic computational codes and is generalizable to any number of discretization variables. Building on ideas introduced in the Robust Verification [2] approach, we estimate the converged solution and orders of convergence with uncertainty using multiple fits of a discretization error model. In contrast to Robust Verification, we perform these fits to many bootstrap samples yielding a larger set of predictions with smoother statistics. Here, bootstrap resampling is performed on the lack-of-fit errors for deterministic code responses, and directly on the noisy data set for stochastic responses. This approach lends a degree of robustness to the overall results, capable of yielding precise verification results for sufficiently resolved data sets, and appropriately expanding the uncertainty when the data set does not support a precise result. For stochastic responses, a credibility assessment is also performed to give the analyst an indication of the trustworthiness of the results. Furthermore, this approach is suitable for both code and solution verification, and is particularly useful for solution verification of high-consequence simulations..

97 MATHEMATICS AND COMPUTING↗

The Atacama Cosmology Telescope: Large-scale velocity reconstruction with the kinematic Sunyaev-Zel'dovich effect and DESI LRGs

The kinematic Sunyaev-Zel'dovich (kSZ) effect induces a non-zero density-density-temperature bispectrum, which we can use to reconstruct the large-scale velocity field from a combination of cosmic microwave background (CMB) and galaxy density measurements, in a procedure known as “kSZ velocity reconstruction”. This method has been forecast to constrain large-scale modes with future galaxy and CMB surveys, improving their measurement beyond what is possible with the galaxy surveys alone. Such measurements will enable tighter constraints on large-scale signals such as primordial non-Gaussianity, deviations from homogeneity, and modified gravity. In this work, we demonstrate a statistically significant measurement of kSZ velocity reconstruction for the first time, by applying quadratic estimators to the combination of the ACT DR6 CMB+kSZ map and the DESI LRG galaxies (with photometric redshifts) in order to reconstruct the velocity field. We do so using a formalism appropriate for the 2-dimensional projected galaxy fields that we use, which naturally incorporates the curved-sky effects important on the largest scales. We find evidence for the signal by cross-correlating with an external estimate of the velocity field from the spectroscopic BOSS survey and rejecting the null (no-kSZ) hypothesis at 3.8σ. Our work presents a first step towards the use of this observable for cosmological analyses.

Sunyaev-Zeldovich effect↗

Effect of the volume fraction gradient on the phase interaction force model for disperse two-phase flows

In this work, the effects of the particle volume fraction gradient on fluid-particle interactions are studied. The phase interaction force is decomposed into three terms. For the first term, namely the symmetrized force density, we present theoretical reasoning and numerical evidence to assume that it is independent of the particle volume fraction gradient. The second term is the particle volume fraction gradient times a newly introduced diffusion stress. The third term is the divergence of the particle-fluid-particle (PFP) stress. If this assumption of independence of the particle volume fraction gradient for the first term can be verified, to the first order of the ratio of the mean distance between particles to the macroscopic lengthscale, all three terms can be studied and modeled in flows with uniform particle distributions. Models thus obtained are applicable to statistically inhomogeneous flows, with the second and third terms accounting for statistical inhomogeneity. To verify this assumption, numerical simulations of flows passing fixed arrays of particles are performed. Both uniform and nonuniform particle volume fractions are studied and compared for disperse multiphase flows with the particle Reynolds numbers ranging from 1 to 100, and particle volume fraction ranging from 1% to 26% in statistically steady states. It is found that the symmetrized force (first) term can be well approximated by the drag force obtained from studies of uniform flows. The diffusion stress is positive along the flow direction and negative in the directions perpendicular to the flow. In the case of moving particles, this stress could potentially cause particle aggregation in the flow direction and dispersion in the directions perpendicular to the flow. Finally, the diffusion stress is only important when there is a volume fraction gradient, while the PFP stress can be important in inhomogeneous flows with either nonuniform particle concentrations or nonuniform average relative velocities between the phases.

42 ENGINEERING↗

Moments of parton distribution functions of the pion from lattice QCD using gradient flow

We present a nonperturbative determination of the pion valence parton distribution function (PDF) moment ratios ⟨𝑥 𝑛−1 ⟩/⟨𝑥⟩ up to 𝑛 = 6, using the gradient flow in lattice quantum chromodynamics (QCD). As a testing ground, we employ SU(3) isosymmetric gauge configurations generated by the OpenLat initiative with a pseudoscalar mass of 𝑚 𝜋 ≃ 411 MeV. Our analysis uses four lattice spacings and a nonperturbatively improved action, enabling full control over the continuum extrapolation, and the limit of vanishing flow time, 𝑡 →0. The flowed ratios exhibit O(𝑎 2 ) scaling across the ensembles, and the continuum-extrapolated results, matched to the $\overline{MS}$ scheme at 𝜇 = 2 GeV using next-to-next-to-leading order matching coefficients, show only mild residual flow-time dependence. The resulting ratios, computed with a relatively small number of configurations, are consistent with phenomenological expectations for the pion’s valence distribution, with statistical uncertainties that are competitive with modern global fits. These findings demonstrate that the gradient flow provides an efficient and systematically improvable method to access partonic quantities from first principles. Future extensions of this work will target lighter pion masses toward the physical point, and applications to nucleon structure such as the proton PDFs and the gluon and sea-quark distributions.

lattice QCD↗

A consistent, Bayesian, approach to the cross section probability distribution in the unresolved resonance region

The cross sections of neutron-induced reactions can be divided into three energy ranges: the resolved resonance region (RRR), the unresolved resonance region (URR), and the fast region. In general, the cross sections in the URR show significant fluctuations that cannot be predicted and cannot be experimentally resolved, thus, it is commonly assumed that the cross section at a specific energy is given by a probability distribution function (PDF) over a range of values that can span several orders of magnitude. The current methodology used to describe such behavior is to construct the PDF by stochastically generating resonance ladders and numerically measuring the PDF. The resonance ladders are sampled using known resonance statistical properties and average resonance widths and spacings extrapolated from the RRR. Although this is a standard and widely used technique, it is computationally very expensive, therefore, an alternative, analytical, approach would be preferable due to the considerable speed up of the computational time in real life applications. Moreover, the current methodology does not take into account existing experimental data, such for total and capture cross sections, that are available for many nuclei. Finally, this approach was developed to be used in reactor-scale applications and it is not suited for use in single-event applications. In this work we will rethink the entire approach to the PDF construction using a Bayesian mindset. This will allow us to provide a different definition of the PDF that allows a much faster calculation of the higher-temperature PDFs and a proper combination of theoretical and experimental PDFs following the probability theory. We will also show that our definition is well suited for single-event applications and we will make an explicit connection between our method and the standard approach. We do this by showing that the central limit theorem applies and our method leads to the same PDF obtained with the standard methodology, for a large number of events per history.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Systematic Study of the Self-Renormalized Nucleon Gluon PDF in Large-Momentum Effective Theory

We present a systematic study of the nucleon gluon parton distribution function (PDF) using the self-renormalized large-momentum effective theory (LaMET) approach in lattice QCD. This work extends previous gluon-PDF extractions by performing a detailed analysis of key systematic effects, including gauge-link smearing, lattice spacing, pion mass, and nucleon boost momentum. The self-renormalization framework mitigates ultraviolet divergences associated with Wilson-line self-energy and renormalon contributions by combining lattice matrix elements with perturbative short-distance information, thereby preserving the correct infrared structure. Calculations are performed on $N_f=2+1+1$ HISQ ensembles generated by the MILC Collaboration at three lattice spacings and two pion masses, with boosted nucleon states reaching momenta up to 2.2~GeV. We determine renormalization factors from zero-momentum matrix elements and apply hybrid renormalization to suppress discretization artifacts. After extrapolating large-separation behavior and performing Fourier transforms, we reconstruct quasi-PDFs and match them to lightcone PDFs using next-to-leading order Wilson coefficients. Our results demonstrate that smearing and lattice-spacing effects are under control, and pion-mass and lattice-spacing dependence is mild relative to the current $O(10^6)$ statistics; however, momentum dependence remains a significant source of uncertainty. Future work including even larger boost momenta will be essential to reduce systematics in lattice determinations of the gluon PDF and to advance toward precision QCD phenomenology at the LHC and the future Electron-Ion Collider.

FOS: Physical sciences↗

Quantifying experimental edge plasma evolution via multidimensional adaptive Gaussian process regression

The edge density and temperature of tokamak plasmas are strongly correlated with energy and particle confinement and their quantification is fundamental to understanding edge dynamics. These quantities exhibit behaviours ranging from sharp plasma gradients and fast transient phenomena (e.g. transitions between low and high confinement regimes) to nominal stationary phases. Analysis of experimental edge measurements therefore require robust fitting techniques to capture potentially stiff spatiotemporal evolution. Additionally, fusion plasma diagnostics inevitably involve measurement errors and data analysis requires a statistical framework to accurately quantify uncertainties. This paper outlines a generalized multidimensional adaptive Gaussian process routine capable of automatically handling noisy data and spatiotemporal correlations. We focus on the edge-pedestal region in order to underline advancements in quantifying time-dependent plasma profiles including transport barrier formation on the Alcator C-Mod tokamak.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Solar cosmic ray hazard to interplanetary and earth-orbital space travel

A statistical treatment of the radiation hazards to astronauts due to solar cosmic ray protons is reported to determine shielding requirements for solar proton events. More recent data are incorporated into the present analysis in order to improve the accuracy of the predicted mission fluence and dose. The effects of the finite data sample are discussed. Mission fluence and dose versus shield thickness data are presented for mission lengths up to 3 years during periods of maximum and minimum solar activity; these correspond to various levels of confidence that the predicted hazard will not be exceeded.

Yucker, W. R.↗

Simulation of the velocity diffusion of charged particles in turbulent magnetic fields

Numerical simulation experiments have been performed in order to study the velocity diffusion of charged particles in a static turbulent magnetic field. By following orbits of particles moving in a large ensemble of random magnetic field realizations with suitably-chosen statistical properties, a pitch-angle diffusion coefficient is derived. Results are presented for a variety of particle rigidities and r.m.s. random field strengths and a power spectrum typical of conditions found in interplanetary space.

Kaiser, T. B.↗

Trends in total column ozone measurements

It is important to ensure the best available data are used in any determination of possible trends in total ozone in order to have the most accurate estimates of any trends and the associated uncertainties. Accordingly, the existing total ozone records were examined in considerable detail. Once the best data set has been produced, the statistical analysis must examine the data for any effects that might indicate changes in the behavior of global total ozone. The changes at any individual measuring station could be local in nature, and herein, particular attention was paid to the seasonal and latitudinal variations of total ozone, because two dimensional photochemical models indicate that any changes in total ozone would be most pronounced at high latitudes during the winter months. The conclusions derived from this detailed examination of available total ozone can be split into two categories, one concerning the quality and the other the statistical analysis of the total ozone record.

Rowland, F. S.↗