Search NASA⌕ Search

SEARCH · Search NASA

Results for “Probability Distribution Function”

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 469 records · Page 26
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

Derivation of a Multiparameter Gamma Model for Analyzing the Residence-Time Distribution Function for Nonideal Flow Systems as an Alternative to the Advection-Dispersion Equation

A new residence-time distribution (RTD) function has been developed and applied to quantitative dye studies as an alternative to the traditional advection-dispersion equation (AdDE). The new method is based on a jointly combined four-parameter gamma probability density function (PDF). The gamma residence-time distribution (RTD) function and its first and second moments are derived from the individual two-parameter gamma distributions of randomly distributed variables, tracer travel distance, and linear velocity, which are based on their relationship with time. The gamma RTD function was used on a steady-state, nonideal system modeled as a plug-flow reactor (PFR) in the laboratory to validate the effectiveness of the model. The normalized forms of the gamma RTD and the advection-dispersion equation RTD were compared with the normalized tracer RTD. The normalized gamma RTD had a lower mean-absolute deviation (MAD) (0.16) than the normalized form of the advection-dispersion equation (0.26) when compared to the normalized tracer RTD. The gamma RTD function is tied back to the actual physical site due to its randomly distributed variables. The results validate using the gamma RTD as a suitable alternative to the advection-dispersion equation for quantitative tracer studies of non-ideal flow systems.

Embry, Irucka↗
DOE OSTI2016

Universal scaling in the aging of the strong glass former SiO{sub 2}

We show that the aging dynamics of a strong glass former displays a strikingly simple scaling behavior, connecting the average dynamics with its fluctuations, namely, the dynamical heterogeneities. We perform molecular dynamics simulations of SiO{sub 2} with van Beest-Kramer-van Santen interactions, quenching the system from high to low temperature, and study the evolution of the system as a function of the waiting time t{sub w} measured from the instant of the quench. We find that both the aging behavior of the dynamic susceptibility χ{sub 4} and the aging behavior of the probability distribution P(f{sub s,r}) of the local incoherent intermediate scattering function f{sub s,r} can be described by simple scaling forms in terms of the global incoherent intermediate scattering function C. The scaling forms are the same that have been found to describe the aging of several fragile glass formers and that, in the case of P(f{sub s,r}), have been also predicted theoretically. A thorough study of the length scales involved highlights the importance of intermediate length scales. We also analyze directly the scaling dependence on particle type and on wavevector q and find that both the average and the fluctuations of the slow aging dynamics are controlled by a unique aging clock, which is not only independent of the wavevector q, but is also the same for O and Si atoms.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Probability of Pinching Failure in PWR SNF after Long-Term Storage

A probabilistic assessment was performed on pinching failure of PWR SNF rod cladding under accidental SNF cask transport conditions. a set of 3,000 rod pinching scenarios was analyzed, with different geometrical, material, and loading parameters, using a Latin Hypercube Sampling (LHS). The investigated accidents include train derailments leading to cask collisions or cask drops. Pinching failure was related to cladding offset strain, defined as the plastic deformation interval in radial cladding direction normalized with respect to the outer cladding diameter. Finite element (FE) models were used to compute SNF rod deformations and the offset strain demand OSD in each of the assessed pinching scenarios. The numerical models account for the contribution of the degraded fuel pellets to rod pinching stiffness. The offset strain capacity OSC was computed using regression and interpolation models based on ring compression (RCT) test results, and account for hydride-related effects on cladding ductility. The results show a low probability of pinching failure. Hydride-related cladding embrittlement affects the structural integrity of PWR SNF rods after long-term dry storage of up to 300 years. Radial hydride patterns or large amounts of circumferentially oriented hydrides reduce the cladding ductile deformation capability when pinched at low material temperatures. The probability of hydride-related pinching failure of PWR cladding is estimated after long-term dry storage to address: How do PWR SNF rod conditions like hydrogen content, vacuum drying peak cladding hoop stress σ{sub θ}, and cladding material temperature affect cladding ductility? Is hydride-related cladding embrittlement a material degradation mechanism that affects the structural integrity of PWR SNF rods under pinching loading? What are the controlling parameters for deformation demands in PWR SNF rods under pinching loading? Is the pinching load amplitude an important parameter that controls the probability of pinching failure? How do parameters like pellet-cladding gap or fuel pellet degradation affect the pinching loading capacity of PWR SNF rods? Statistical modelling of PWR SNF rod conditions in U.S. inventory: A pseudo rod database was generated based on EIA database. Analysis of train accident scenarios: a probability density function (pdf) of a pinching load distribution P{sub 9m} of a 9-m cask drop is linearly scaled to generate a more representative: distribution P{sub AS}. The new distribution considers train collision velocities, cask drops for on- and off-bridge accidents, and cask impact orientation. LHS generates 3,000 rod pinching scenarios based on nine independent geometrical, material, and loading parameters. Probabilistic FE Analysis: The LHS realizations are used to generate 3,000 FE models. The FE analysis is conducted for P{sub 9m}, and P{sub AS}. The workflow is automatized. The OSC is computed with regression and interpolation models, as a function of the hydrogen content C{sub H}, vacuum drying peak σ{sub θ}, and SNF temperature T{sub FR}: The regression data was obtained in RCTs on PWR cladding. The OSD is derived from the cladding deformations in the FE models. Cladding failure was identified via a comparison of OSC and OSD. A 300 and a 30 years-of-dry-storage duration is investigated: The dry-storage period affects the expected SNF temperature and consequently, the OSC and the OSD. The material temperature LHS population is adapted to account for a SNF temperature variation. The probability of hydride-related, ductile cladding failure during a severe transportation accident is not significantly affected by the pinching load amplitude: The probability is slightly reduced from a range of less than 0.0333% to 0.1667% for P{sub 9m}, to a range of less than 0.0333% to 0.1333% for P{sub AS}, depending on the vacuum drying peak cladding hoop stress modelling assumptions. Further, the SNF dry-storage duration does not significantly affect the probability of ductile pinching failure. Based on available train accidents data, the pinching failure probability for a 5,000 km long transportation route ranges from less than 1.62 E-7 to a value of 6.47 E-7, depending on the vacuum drying peak cladding hoop stress modelling assumptions. The pinching load amplitude is not a controlling parameter for the probability of ductile cladding pinching failure, and does not control the cladding deformations once pellet-cladding contact is established, The pellet-cladding gap width, however, plays a key role. Fuel pellet degradation only affects the deformation capacity, if the pellet stiffness is significantly reduced (more than 90%). Only highly embrittled cladding is likely to fail. The deformation demands on the cladding under pinching loading are small.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2010

A UNIFIED THEORY FOR THE EFFECTS OF STELLAR PERTURBATIONS AND GALACTIC TIDES ON OORT CLOUD COMETS

We examine the effects of passing field stars on the angular momentum of a nearly radial orbit of an Oort cloud comet bound to the Sun. We derive the probability density function of the change in angular momentum from one stellar encounter, assuming a uniform and isotropic field of perturbers. We show that the total angular momentum follows a Levy flight, and determine its distribution function. If there is an asymmetry in the directional distribution of perturber velocities, the marginal probability distribution of each component of the angular momentum vector can be different. The constant torque attributed to Galactic tides arises from a non-cancellation of perturbations with an impact parameter of order the semimajor axis of the comet. When the close encounters are rare, the angular momentum is best modeled by the stochastic growth of stellar encounters. If trajectories passing between the comet and the Sun occur frequently, the angular momentum exhibits the coherent growth attributed to the Galactic tides.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

A behaviorally-explicit approach for delivering vaccine baits to mesopredators to control epizootics in fragmented landscapes

Despite the widespread use of aerial baiting to manage epizootics among free-ranging populations,particularly in rabies management, bait acceptance and seroconversion rates often are lower than required to eliminate spread of disease. Our objectives in this study, therefore, were to evaluate the performance of stratified bait distribution models derived from resource selection functions (RSF) on uptake of placebo rabies baits by raccoons ( Procyon lotor ) and Virginia opossums ( Didelphis virginiana ), as well as the probability of bait uptake as a function of proximity to bait distribution areas in fragmented agricultural ecosystems. Among 478 raccoons and 108 opossums evaluated for presence of Rhodamine B (RB) across 8 sites, only 26% of raccoons and 20% of opossums exhibited marking consistent with bait consumption 14–24 days post-baiting. The effective area treated, based on 90% kernel density estimators of marked individuals, ranged from 99–240 ha larger than bait distribution zones, with RB marked individuals captured up to 753m beyond the bait zone. Despite incorporation of RSF data into bait distribution models, no differences in uptake rates were observed between treatment and control sites. These data likely reflect the underlying constraints imposed by the loss and fragmentation of habitat on animal movement in heterogeneous landscapes, forcing individuals to optimize movements at coarse (i.e., patch-level) rather than fine spatial scales in highly fragmented environments. Our data also confirm that the probability of bait acceptance decreases with increasing distance from bait zone interiors, even within the zone itself. Thus, although bait acceptance was confirmed beyond bait zone boundaries, the proportion of vaccinated individuals may comprise a small minority of the population at increasing distances from baiting interiors. These data suggest focal baiting creates a buffered area of treated individuals around bait zones or bait stations, but repeated treatments may be needed to achieve sufficient uptake to eradicate disease.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2015

Note: Precise radial distribution of charged particles in a magnetic guiding field

Current high precision beta decay experiments of polarized neutrons, employing magnetic guiding fields in combination with position sensitive and energy dispersive detectors, resulted in a detailed study of the mono-energetic point spread function (PSF) for a homogeneous magnetic field. A PSF describes the radial probability distribution of mono-energetic electrons at the detector plane emitted from a point-like source. With regard to accuracy considerations, unwanted singularities occur as a function of the radial detector coordinate which have recently been investigated by subdividing the radial coordinate into small bins or employing analytical approximations. In this note, a series expansion of the PSF is presented which can numerically be evaluated with arbitrary precision.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2024

A Eulerian three-fluid flow framework for the study of fuel dispersal behavior under loss-of-coolant accident conditions

As the nuclear industry explores new fuel designs to accommodate increased burnup, studying fuel behavior during loss-of-coolant accidents is essential for ensuring the safe operation of light-water reactors. Here, this study presents a computational model to simulate the intricate three-phase flow dynamics of fuel dispersal following cladding breach, where solid fuel particles and carrier fission gases interact with the surrounding liquid or gas medium. The model utilizes a Eulerian three-fluid framework to depict the phase interaction within interpenetrating continua, treating fragmented fuel as granular material. Closure relationships for bulk-flow properties of the solid phase and detailed interfacial exchange terms in high-temperature, high-speed gas-solid-liquid flow scenarios are proposed. The model was validated using adiabatic experimental results in the literature, focusing on characterizing fuel dispersal behavior during loss-of-coolant accidents conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched probability density functions reported in experimental studies. Key highlights of this study include the theoretical description of a coupled three-phase flow with phase change, interpretation of the high-pressure boundary conditions and insights into the transient behavior of gas-solid-liquid phase dynamics during a simulated high-speed dispersal event.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2015

Fickian dispersion is anomalous

The thesis put forward here is that the occurrence of Fickian dispersion in geophysical settings is a rare event and consequently should be labeled as anomalous. What people classically call anomalous is really the norm. In a Lagrangian setting, a process with mean square displacement which is proportional to time is generally labeled as Fickian dispersion. With a number of counter examples we show why this definition is fraught with difficulty. In a related discussion, we show an infinite second moment does not necessarily imply the process is super dispersive. By employing a rigorous mathematical definition of Fickian dispersion we illustrate why it is so hard to find a Fickian process. We go on to employ a number of renormalization group approaches to classify non-Fickian dispersive behavior. Scaling laws for the probability density function for a dispersive process, the distribution for the first passage times, the mean first passage time, and the finite-size Lyapunov exponent are presented for fixed points of both deterministic and stochastic renormalization group operators. The fixed points of the renormalization group operators are p-self-similar processes. A generalized renormalization group operator is introduced whose fixed points form a set of generalized self-similar processes. Finally, power-law clocks are introduced to examine multi-scaling behavior. Several examples of these ideas are presented and discussed.

58 GEOSCIENCES↗
DOE OSTI2015

Thermal Recovery of Colloidal Quantum Dot Ensembles Following Photoinduced Dimming

Colloidal CdSe quantum dot (QD) core ensembles were photodimmed and allowed to recover in the dark using ambient thermal energy at a range of temperatures. Nonlinear thermal recovery is well described by a stretched exponential function, and further analysis yields an underlying probability distribution of rate constants. Casting the rate constants as a collection of first-order activated processes provides an activation barrier probability distribution with significant density at room-temperature thermal energy that peaks at 200 meV before decaying to zero. This treatment for the recovery transition intuitively describes the distributed kinetics observed and complements commonly proposed blinking mechanisms.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Convective shells in the interior of Cepheid variable stars: Overshooting models based on hydrodynamic simulations

Context. Because Cepheid variable stars have long been used as a cosmic benchmark for scaling distances in our Galaxy and beyond, the accuracy of stellar evolution models for Cepheids have wide-reaching effects. However, our understanding of the dynamics in the interiors of these physically complex stars is limited. Aims. Our goal is to provide a detailed multi-dimensional picture of hydrodynamic convection and convective boundary mixing in the interior of Cepheids. Methods. Using the Modules for Experiments in Stellar Astrophysics (MESA), we studied the structure of intermediate-mass stars that cross the instability strip. Then, we performed two-dimensional hydrodynamic simulations of six stars with the fully compressible Multidimensional Stellar Implicit Code (MUSIC). Our simulations did not model the radial pulsations but focused on the interior structure of this family of stars. We developed and applied a new statistical analysis to examine convection and convective boundary mixing in the interior of these stellar simulations. Results. Based on a grid of MESA models, we demonstrated that a common structure for intermediate mass Cepheids includes an interior convective shell as well as a thin outer convective envelope. Using the extreme value theory approach to analyze our MUSIC simulation data, we found that overshooting above the convective shell fills the space between these convectively unstable layers. We developed a new statistical analysis that provides a clearer picture of how overshooting fills this layer; it also allowed us to formulate a detailed comparison between overshooting above and below the convective shell. Our analysis effectively decomposes the overshooting layer into two layers: a weak overshooting layer and a strong overshooting layer. Statistically, this is accomplished by decomposing the strongly non-Gaussian probability density function into a mixture of gamma distributions. Using our mixture model, we showed that the ratio of overshooting lengths above and below the convective shell depends directly on the radial extent of the convective shell as well as its depth in the star. We proposed a new form for the diffusion coefficient that addresses the need for overlapping overshooting layers between convective shells. We introduced the idea of a “super-mixing layer” where overshooting from both the convective shell and the convective envelope results in efficient mixing and could be viewed as merging the two adjacent convective zones.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Data-Driven Benchmarking of Building Energy Efficiency Utilizing Statistical Frontier Models

Frontier methods quantify the energy efficiency of buildings by forming an efficient frontier (best-practice technology) and by comparing all buildings against that frontier. Because energy consumption fluctuates over time, the efficiency scores are stochastic random variables. Existing applications of frontier methods in energy efficiency either treat efficiency scores as deterministic values or estimate their uncertainty by resampling from one set of measurements. Availability of smart meter data (repeated measurements of energy consumption of buildings) enables using actual data to estimate the uncertainty in efficiency scores. Additionally, existing applications assume a linear form for an efficient frontier; i.e.,they assume that the best-practice technology scales up and down proportionally with building characteristics. However, previous research shows that buildings are nonlinear systems. This paper proposes a statistical method called stochastic energy efficiency frontier (SEEF) to estimate a bias-corrected efficiency score and its confidence intervals from measured data. The paper proposes an algorithm to specify the functional form of the frontier, identify the probability distribution of the efficiency score of each building using measured data, and rank buildings based on their energy efficiency. To illustrate the power of SEEF, this paper presents the results from applying SEEF on a smart meter data set of 307 residential buildings in the United States. SEEF efficiency scores are used to rank individual buildings based on energy efficiency, to compare subpopulations of buildings, and to identify irregular behavior of buildings across different time-of-use periods. SEEF is an improvement to the energy-intensity method (comparing kWh/sq.ft.): whereas SEEF identifies efficient buildings across the entire spectrum of building sizes, the energy-intensity method showed bias toward smaller buildings. The results of this research are expected to assist researchers and practitioners compare and rank (i.e.,benchmark) buildings more robustly and over a wider range of building types and sizes. Eventually, doing so is expected to result in improved resource allocation in energy-efficiency programs.

Kavousian, A↗
DOE OSTI2014

Transit and radial velocity survey efficiency comparison for a habitable zone Earth

Transit and radial velocity searches are two techniques for identifying nearby extrasolar planets to Earth that transit bright stars. Identifying a robust sample of these exoplanets around bright stars for detailed atmospheric characterization is a major observational undertaking. In this study we describe a framework that answers the question of whether a transit or radial velocity survey is more efficient at finding transiting exoplanets given the same amount of observing time. Within the framework we show that a transit survey's window function can be approximated using the hypergeometric probability distribution. We estimate the observing time required for a transit survey to find a transiting Earth-sized exoplanet in the habitable zone (HZ) with an emphasis on late-type stars. We also estimate the radial velocity precision necessary to detect the equivalent HZ Earth-mass exoplanet that also transits when using an equal amount of observing time as the transit survey. We find that a radial velocity survey with σ{sub rv} ∼ 0.6 m s{sup –1} precision has comparable efficiency in terms of observing time to a transit survey with the requisite photometric precision σ{sub phot} ∼ 300 ppm to find a transiting Earth-sized exoplanet in the HZ of late M dwarfs. For super-Earths, a σ{sub rv} ∼ 2.0 m s{sup –1} precision radial velocity survey has comparable efficiency to a transit survey with σ{sub phot} ∼ 2300 ppm.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

STAR CLUSTER PROPERTIES IN TWO LEGUS GALAXIES COMPUTED WITH STOCHASTIC STELLAR POPULATION SYNTHESIS MODELS

We investigate a novel Bayesian analysis method, based on the Stochastically Lighting Up Galaxies (slug) code, to derive the masses, ages, and extinctions of star clusters from integrated light photometry. Unlike many analysis methods, slug correctly accounts for incomplete initial mass function (IMF) sampling, and returns full posterior probability distributions rather than simply probability maxima. We apply our technique to 621 visually confirmed clusters in two nearby galaxies, NGC 628 and NGC 7793, that are part of the Legacy Extragalactic UV Survey (LEGUS). LEGUS provides Hubble Space Telescope photometry in the NUV, U, B, V, and I bands. We analyze the sensitivity of the derived cluster properties to choices of prior probability distribution, evolutionary tracks, IMF, metallicity, treatment of nebular emission, and extinction curve. We find that slug's results for individual clusters are insensitive to most of these choices, but that the posterior probability distributions we derive are often quite broad, and sometimes multi-peaked and quite sensitive to the choice of priors. In contrast, the properties of the cluster population as a whole are relatively robust against all of these choices. We also compare our results from slug to those derived with a conventional non-stochastic fitting code, Yggdrasil. We show that slug's stochastic models are generally a better fit to the observations than the deterministic ones used by Yggdrasil. However, the overall properties of the cluster populations recovered by both codes are qualitatively similar.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Uncertainty Visualization of Critical Points of 2D Scalar Fields for Parametric and Nonparametric Probabilistic Models

This paper presents a novel end-to-end framework for closed-form computation and visualization of critical point uncertainty in 2D uncertain scalar fields. Critical points are fundamental topological descriptors used in the visualization and analysis of scalar fields. The uncertainty inherent in data (e.g., observational and experimental data, approximations in simulations, and compression), however, creates uncertainty regarding critical point positions. Uncertainty in critical point positions, therefore, cannot be ignored, given their impact on downstream data analysis tasks. Here, in this work, we study uncertainty in critical points as a function of uncertainty in data modeled with probability distributions. Although Monte Carlo (MC) sampling techniques have been used in prior studies to quantify critical point uncertainty, they are often expensive and are infrequently used in production-quality visualization software. We, therefore, propose a new end-to-end framework to address these challenges that comprises a threefold contribution. First, we derive the critical point uncertainty in closed form, which is more accurate and efficient than the conventional MC sampling methods. Specifically, we provide the closed-form and semianalytical (a mix of closed-form and MC methods) solutions for parametric (e.g., uniform, Epanechnikov) and nonparametric models (e.g., histograms) with finite support. Second, we accelerate critical point probability computations using a parallel implementation with the VTK-m library, which is platform portable. Finally, we demonstrate the integration of our implementation with the ParaView software system to demonstrate near-real-time results for real datasets.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2014

Physical Constraints on Geologic CO2 Sequestration in Low-Volume Basalt Formations

Deep basalt formations within large igneous provinces have been proposed as target reservoirs for carbon capture and sequestration on the basis of favorable CO2-water-rock reaction kinetics that suggest carbonate mineralization rates on the order of 102–103 d. Although these results are encouraging, there exists much uncertainty surrounding the influence of fracture-controlled reservoir heterogeneity on commercial-scale CO2 injections in basalt formations. This work investigates the physical response of a low-volume basalt reservoir to commercial-scale CO2 injections using a Monte Carlo numerical modeling experiment such that model variability is solely a function of spatially distributed reservoir heterogeneity. Fifty equally probable reservoirs are simulated using properties inferred from the deep eastern Snake River Plain aquifer in southeast Idaho, and CO2 injections are modeled within each reservoir for 20 yr at a constant mass rate of 21.6 kg s–1. Results from this work suggest that (1) formation injectivity is generally favorable, although injection pressures in excess of the fracture gradient were observed in 4% of the simulations; (2) for an extensional stress regime (as exists within the eastern Snake River Plain), shear failure is theoretically possible for optimally oriented fractures if Sh is less than or equal to 0.70SV; and (3) low-volume basalt reservoirs exhibit sufficient CO2 confinement potential over a 20 yr injection program to accommodate mineral trapping rates suggested in the literature.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Development and Integration of a Stochastic Clad Damage Propagation Model into PRONGHORN-SC Subchannel Analysis Code

The failure of fuel pins in nuclear reactors is intrinsically stochastic. Typically, a combination of variation in manufacturing that affects the material characteristics and the fuel assembly dimensions, variation in operating conditions, such as local power, coolant flow rate, and irradiation induced changes in material properties lead to a large uncertainty in failure margin of the fuel pins. Failure, therefore, may occur in exceptional pins with adverse combinations of these variations. Upon a metal fuel pin (U-Pu-Zr/HT9) failure, depressurization of the fuel pin takes place by release of fission gas, liquid sodium bond, and potentially solid fuel particles or molten/eutectic fuel droplets through the hole in cladding. The effect of a fission gas jet on neighbor fuel pins and possible propagation of a clad damage during normal operation was studied experimentally in 1970s and it was found that the post-failure fission gas jet insulates the jet impingement area of the target fuel pin surface and could increase the target pin’s surface temperature by as much as 100 – 200 K during the failed pin depressurization. It was concluded that the effect should not lead to fuel pin failure propagation during normal operation. In accident scenarios of sodium and lead fast reactors such as Unprotected Loss-Of-Flow (ULOF) or Unprotected Transient Over Power (UTOP), the fuel pins can be subjected to higher clad temperatures and fuel pin pressures or fuel clad mechanical/chemical interaction where thermal creep margin becomes significantly lower compared to the normal operation conditions. Therefore, possible stochastic failure and the post-failure fission gas/fuel jet impingement could be critical in order to predict fuel pin failure propagation. Pin depressurization due to fission gas release may degrade the heat transfer by formation of a gas blanket on a neighboring pin surface, which is a local phenomenon, and by causing coolant flow deceleration and starvation, which could affect a surrounding region as well. Furthermore, the potential presence of solid fuel particles or molten fuel at the time of clad failure could boost post-failure jet induced degradation even further. The present study models the U-Pu-Zr/HT9 metal fuel pin failure and stochastic clad damage propagation by biased sampling based on a Cumulative Damage Fraction (CDF) type clad failure criterion and the normal distribution of fuel failure probability density as a function of logarithm of Cumulative Damage Fraction. In addition, the effect of post-failure fission gas jet on heat transfer degradation is modeled for the target pins. This model is called stochastic Clad Damage Propagation (CDAP). The CDAP model is now fully integrated into developmental version of PRONGHORN-SC subchannel analysis code, allowing for modeling local failures and its propagation potential. Section 2 describes the components of the CDAP models. Section 3 describes the model implementation to PRONGHORN-SC and input specifications. Section 4 describes the CDAP model validation coupled to PRONGHORN-SC.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2011

Wave functions of the proton ground state in the presence of a uniform background magnetic field in lattice QCD

We calculate the probability distributions of quarks in the ground state of the proton, and how they are affected in the presence of a constant background magnetic field. We focus on wave functions in the Landau and Coulomb gauges. We observe the formation of a scalar u-d diquark clustering. The overall distortion of the quark probability distribution under a very large magnetic field, as demanded by the quantization conditions on the field, is quite small. The effect is to elongate the distributions along the external field axis while localizing the remainder of the distribution.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2012

Respiratory motion sampling in 4DCT reconstruction for radiotherapy

Purpose: Phase-based and amplitude-based sorting techniques are commonly used in four-dimensional CT (4DCT) reconstruction. However, effect of these sorting techniques on 4D dose calculation has not been explored. In this study, the authors investigated a candidate 4DCT sorting technique by comparing its 4D dose calculation accuracy with that for phase-based and amplitude-based sorting techniques.Method: An optimization model was formed using organ motion probability density function (PDF) in the 4D dose convolution. The objective function for optimization was defined as the maximum difference between the expected 4D dose in organ of interest and the 4D dose calculated using a 4DCT sorted by a candidate sampling method. Sorting samples, as optimization variables, were selected on the respiratory motion PDF assessed during the CT scanning. Breathing curves obtained from patients' 4DCT scanning, as well as 3D dose distribution from treatment planning, were used in the study. Given the objective function, a residual error analysis was performed, and k-means clustering was found to be an effective sampling scheme to improve the 4D dose calculation accuracy and independent with the patient-specific dose distribution. Results: Patient data analysis demonstrated that the k-means sampling was superior to the conventional phase-based and amplitude-based sorting and comparable to the optimal sampling results. For phase-based sorting, the residual error in 4D dose calculations may not be further reduced to an acceptable accuracy after a certain number of phases, while for amplitude-based sorting, k-means sampling, and the optimal sampling, the residual error in 4D dose calculations decreased rapidly as the number of 4DCT phases increased to 6.Conclusion: An innovative phase sorting method (k-means method) is presented in this study. The method is dependent only on tumor motion PDF. It could provide a way to refine the phase sorting in 4DCT reconstruction and is effective for 4D dose accumulation. Optimized sorting techniques could achieve acceptable residuals (less than 0.5% of the prescription dose) using 6 sorting samples, which is much better than amplitude-based or phase-based sorting. Further increase in sorting phase number exceeding 6 or more may not be necessary when using the k-means sampling or optimal sampling points.

62 RADIOLOGY AND NUCLEAR MEDICINE↗