Search NASA⌕ Search

SEARCH · Search NASA

Results for “Particle identication methods”

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.

137 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Establishing metrics to quantify spatial similarity in spherical and red blood cell distributions

As computational power increases and systems with millions of red blood cells can be simulated, it is important to note that varying spatial distributions of cells may affect simulation outcomes. Since a single simulation may not represent the ensemble behavior, many different configurations may need to be sampled to adequately assess the entire collection of potential cell arrangements. In order to determine both the number of distributions needed and which ones to run, we must first establish methods to identify well-generated, randomly placed cell distributions and to quantify distinct cell configurations. We utilize metrics to assess (1) the presence of any underlying structure to the initial cell distribution and (2) similarity between cell configurations. We propose the use of the radial distribution function to identify long-range structure in a cell configuration and apply it to a randomly distributed and structured set of red blood cells. To quantify spatial similarity between two configurations, we make use of the Jaccard index, and characterize sets of red blood cell and sphere initializations. As an extension to our work submitted to the International Conference on Computational Science, we significantly increase our data set size from 72 to 1048 cells, include a similar set of studies using spheres, compare the effects of varying sphere size, and utilize the Jaccard index distribution to probe sets of extremely similar configurations. Our results show that the radial distribution function can be used as a metric to determine long-range structure in both distributions of spheres and RBCs. We determine that the ideal case of spheres within a cube versus bi-concave shaped cells within a cylinder affects the shape of the Jaccard index distributions, as well as the range of Jaccard values, showing that both the shape of particle and the domain may play a role. Furthermore, we also find that the distribution is able to capture very similar configurations through Jaccard index values greater than 95% when appending several nearly identical configurations into the data set.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Solitons in lattice field theories via tight-binding supersymmetry

Reflectionless potentials play an important role in constructing exact solutions to classical dynamical systems (such as the Korteweg-de Vries equation), non-perturbative solutions of various large-N field theories (such as the Gross-Neveu model), and closely related solitonic solutions to the Bogoliubov-de Gennes equations in the theory of superconductivity. These solutions rely on the inverse scattering method, which reduces these seemingly unrelated problems to identifying reflectionless potentials of an auxiliary one-dimensional quantum scattering problem. There are several ways of constructing these potentials, one of which is quantum mechanical supersymmetry (SUSY). In this paper, motivated by recent experimental platforms, we generalize this framework to develop a theory of lattice solitons. We first briefly review the classical inverse scattering method in the continuum limit, focusing on the Korteweg-de Vries (KdV) equation and SU(N) Gross-Neveu model in the large N limit. We then generalize this methodology to lattice versions of interacting field theories. Our analysis hinges on the use of trace identities, which are relations connecting the potential of an equation of motion to the scattering data. For a discrete Schrödinger operator, such trace identities had been known as far back as Toda; however, we derive a new set of identities for the discrete Dirac operator. We then use these identities in a lattice Gross-Neveu and chiral Gross-Neveu (Nambu-Jona-Lasinio) model to show that lattice solitons correspond to reflectionless potentials associated with the discrete scattering problem. These models are of significance as they are equivalent to a mean-field theory of a lattice superconductor. To explicitly construct these solitons, we generalize supersymmetric quantum mechanics to tight-binding models. We show that a matrix transformation exists that maps a tight-binding model to an isospectral one which shares the same structure and scattering properties. The corresponding soliton solutions have both modulated hopping and onsite potential, the former of which has no analogue in the continuum limit. We explicitly compute both topological and non-topological soliton solutions as well as bound state spectra in the aforementioned models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1993

LDEF: Dosimetric measurement results (AO 138-7 experiment)

One of the objectives of the AO 138-7 experiment on board the Long Duration Exposure Facility (LDEF) was a total dose measurement with Thermo Luminescent Detectors (TLD 100). Two identical packages, both of them including five TLD's inside various aluminum shields, are exposed to the space environment in order to obtain the absorbed dose profile. Radiation fluence received during the total mission length was computed, taking into account the trapped particles (AE8 and AP8 models during solar maximum and minimum periods) and the cosmic rays; due to the magnetospheric shielding the solar proton fluences are negligible on the LDEF orbit. The total dose induced by these radiations inside a semi infinite plane shield of aluminum are computed with the radiation transport codes available at DERTS. The dose profile obtained is in good agreement with the evaluation by E.V. Benton. TLD readings are performed after flight; due to the mission duration increase a post flight calibration was necessary in order to cover the range of the in flight induced dose. The results obtained, similar (plus or minus 30 percent) for both packages, are compared with the dose profile computation. For thick shields it seems that the measurements exceed the forecast (about 40 percent). That can be due to a cosmic ray and trapped proton contributions coming from the backside (assumed as perfectly shielded by the LDEF structure in the computation), or to an underestimate of the proton or cosmic ray fluences. A fine structural shielding analysis should be necessary in order to determine the origin of this slight discrepancy between forecast and in flight measurements. For the less shielded dosimeters, mainly exposed to the trapped electron flux, a slight overestimation of the dose (less than 40 percent) appears. Due to the dispersion of the TLD's response, this cannot be confirmed. In practice these results obtained on board LDEF, with less than a factor 1.4 between measurements and forecast, reinforce the validity of the computation methods and models used for the long term evaluation of the radiation levels (flux and dose) encountered in space on low inclination and altitude Earth orbits.

Bourrieau, J.↗
NASA NTRS1990

Noncontact measurement of high temperature using optical fiber sensors

The primary goal of this research program was the investigation and application of noncontact temperature measurement techniques using optical techniques and optical fiber methods. In particular, a pyrometer utilizing an infrared optical light pipe and a multiwavelength filtering approach was designed, revised, and tested. This work was motivated by the need to measure the temperatures of small metallic pellets (approximately 3 mm diameter) in free fall at the Microgravity Materials Processing Drop Tube at NASA Marshall Space Flight Center. In addition, research under this program investigated the adaptation of holography technology to optical fiber sensors, and also examined the use of rare-earth dopants in optical fibers for use in measuring temperature. The pyrometer development effort involved both theoretical analysis and experimental tests. For the analysis, a mathematical model based on radiative transfer principles was derived. Key parameter values representative of the drop tube system, such as particle size, tube diameter and length, and particle temperature, were used to determine an estimate of the radiant flux that will be incident on the face of an optical fiber or light pipe used to collect radiation from the incandescent falling particle. An extension of this work examined the advantage of inclining or tilting the collecting fiber to increase the time that the falling particle remains in the fiber field-of-view. Those results indicate that increases in total power collected of about 15 percent may be realized by tilting the fiber. In order to determine the suitability of alternative light pipes and optical fibers, and experimental set-up for measuring the transmittance and insertion loss of infrared fibers considered for use in the pyrometer was assembled. A zirconium fluoride optical fiber and several bundles of hollow core fiber of varying diameters were tested. A prototype two-color pyrometer was assembled and tested at Virginia Tech, and then tested on the Drop Tube at Marshall Space Flight Center. Radiation from 5 mm diameter niobium drops falling in the Drop Tube was successfully detected, and recorded for later analysis. Subsequent analysis indicated that the imaging of light output from the light pipe onto the detector active areas was not identical for both detectors.

Claus, R. O.↗
DOE OSTI2018

Neutrino Oscillations in Supernovae

A (core-collapse) supernova is the final death call of a massive star during which a gigantic amount of a special kind of elementary particles known as neutrinos are released for a brief period of time. An earlier work by the PI and his collaborators has shown that the dense neutrino gas deep inside the supernova can experience "flavor transformation" or "oscillations" in a collective fashion, a macroscopic quantum phenomenon spanning over distances of 100-1000 kilometers. During this process neutrinos of different energies and propagation directions change their identities simultaneously. This phenomenon may be a key to understanding the mechanism of supernova explosion and the origin of many of the chemical elements found on the earth which were produced in the supernovae exploded in the past. The previous understanding of neutrino oscillations in supernovae was based on a stationary 3-dimensional model with 1 spatial dimension and 2 momentum dimensions. In this "neutrino bulb model" neutrinos are emitted from a spherically symmetric surface and they are decoupled from the motion of matter and the production of chemical elements. In order to accurately model the neutrino oscillations in supernovae and determine its physical consequences, sophisticated 7-dimensional (1 temporal dimension + 3 spatial dimensions + 3 momentum dimensions) modelling of neutrino transport must be performed self-consistently together with the hydrodynamic simulation of the flow of energy and matter inside the supernova and with the calculation of the production and destruction of chemical elements through a large nucleosynthesis network. Although it will be long journey to achieve the lofty goal, our project has made important progress in three different fronts towards this goal. In the first front we explored new computing hardware on which the next generation simulation codes may be run. We ported the neutrino bulb code to the platform of Graphics Processing Unit (GPU) and then later developed a new bulb code that can be run on supercomputers equipped with the Intel Xeon CPUs and/or the Xeon Phi coprocessors. In the second front we developed new algorithms that can be used towards computing neutrino oscillations in multi-dimensional models. We have developed a moment expansion method which is approximately 10 times more efficient than the old algorithm. We have also developed numerical codes to solve neutrino oscillations in toy models with two spatial dimensions. In the last front we used analytic methods to study neutrino oscillations in more sophisticated models, e.g. of multiple spatial dimensions and/or with time dependence. By analyzing the neutrino oscillations in toy models with two spacial dimensions we found that the approximate spacial translation symmetry, which is similar to the spherical symmetry in the neutrino bulb model, can be broken spontaneously as neutrino oscillations begin. By analyzing a time-dependent neutrino bulb model we found that certain modes of neutrino oscillations that are suppressed in the stationary model can become active in the time-dependent model. We have also developed theories which give insights to how ordinary matter may affect neutrino oscillations in supernovae and other astrophysical environments.

Duan, Huaiyu↗
DOE OSTI2024

A prototype scintillator real‐time beam monitor for ultra‐high dose rate radiotherapy

Background: FLASH Radiotherapy (RT) is an emergent cancer RT modality where an entire therapeutic dose is delivered at more than 1000 times higher dose rate than conventional RT. For clinical trials to be conducted safely, a precise and fast beam monitor that can generate out-of-tolerance beam interrupts is required. This paper describes the overall concept and provides results from a prototype ultra-fast, scintillator-based beam monitor for both proton and electron beam FLASH applications. Purpose: A FLASH Beam Scintillator Monitor (FBSM) is being developed that employs a novel proprietary scintillator material. The FBSM has capabilities that conventional RT detector technologies are unable to simultaneously provide: (1) large area coverage; (2) a low mass profile; (3) a linear response over a broad dynamic range; (4) radiation hardness; (5) real-time analysis to provide an IEC-compliant fast beam-interrupt signal based on true two-dimensional beam imaging, radiation dosimetry and excellent spatial resolution. Methods: The FBSM uses a proprietary low mass, less than 0.5 mm water equivalent, non-hygroscopic, radiation tolerant scintillator material (designated HM: hybrid material) that is viewed by high frame rate CMOS cameras. Folded optics using mirrors enable a thin monitor profile of ∼10 cm. A field programmable gate array (FPGA) data acquisition system generates real-time analysis on a time scale appropriate to the FLASH RT beam modality: 100–1000 Hz for pulsed electrons and 10–20 kHz for quasi-continuous scanning proton pencil beams. An ion beam monitor served as the initial development platform for this work and was tested in low energy heavy-ion beams ( 86 Kr +26 and protons). A prototype FBSM was fabricated and then tested in various radiation beams that included FLASH level dose per pulse electron beams, and a hospital RT clinic with electron beams. Results: Results presented in this report include image quality, response linearity, radiation hardness, spatial resolution, and real-time data processing. Furthermore, the HM scintillator was found to be highly radiation damage resistant. It exhibited a small 0.025%/kGy signal decrease from a 216 kGy cumulative dose resulting from continuous exposure for 15 min at a FLASH compatible dose rate of 237 Gy/s. Measurements of the signal amplitude versus beam fluence demonstrate linear response of the FBSM at FLASH compatible dose rates of >40 Gy/s. Comparison with commercial Gafchromic film indicates that the FBSM produces a high resolution 2D beam image and can reproduce a nearly identical beam profile, including primary beam tails. The spatial resolution was measured at 35–40 µm. Tests of the firmware beta version show successful operation at 20 000 Hz frame rate or 50 µs/frame, where the real-time analysis of the beam parameters is achieved in less than 1 µs. Conclusions: The FBSM is designed to provide real-time beam profile monitoring over a large active area without significantly degrading the beam quality. A prototype device has been staged in particle beams at currents of single particles up to FLASH level dose rates, using both continuous ion beams and pulsed electron beams. Using a novel scintillator, beam profiling has been demonstrated for currents extending from single particles to 10 nA currents. Radiation damage is minimal and even under FLASH conditions would require ≥50 kGy of accumulated exposure in a single spot to result in a 1% decrease in signal output. Beam imaging is comparable to radiochromic films, and provides immediate images without hours of processing. Real-time data processing, taking less than 50 µs (combined data transfer and analysis times), has been implemented in firmware for 20 kHz frame rates for continuous proton beams.

2D beam imaging↗
NASA NTRS2020

Analysis of Sample Acquisition Dynamics Using Discrete Element Method

The analysis presented in this paper is conducted in the framework of the Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) project, currently under development at NASA Ames Research Center. OceanWATERS aims at designing a simulation environment which allows for testing autonomy of scientific lander missions to the icy moons of our solar system. Mainly focused on reproducing the end effector interaction with the inherent terrain, this paper introduces a novel discrete element method (DEM)-based approach to determine forces and torques acting on the lander’s scoop during the sample acquisition process. An accurate force feedback from the terrain on the scoop is required by fault-detection and autonomous decision-making algorithms to identify when the requested torque on the robotic arm’s joints exceeds the maximum available torque. Knowledge of the terrain force feedback significantly helps evaluating the arm’s links structural properties and properly selecting actuators for the joints. Models available in literature constitute a partial representation of the dynamics of the interaction. As an example, Balovnev derived an analytical expression of the vertical and horizontal force acting on a bucket while collecting a sample as a function of its geometry and velocity, soil parameters and reached depth. Although the model represents an adequate approximation of the two force components, it ignores the direction orthogonal to the scoop motion and neglects the torque. This work relies on DEM analysis to compensate for analytical models’ deficiencies and inaccuracies, i. e. provide force and torque 3D vectors, defined in the moving reference (body) frame attached to the scoop, at each instant of the sample collection process. Results from the first presented analysis relate to the specific OceanWATERS sampling strategy, which consists of collecting the sample through five consecutive passes with increasing depth, each pass following the same circularlinear- circular trajectory. Data is collected given a specific scoop design interacting with two types of bulk materials, which may characterize the surface of icy planetary bodies: snow and ice. Although specifically concerned with the OceanWATERS design, this first analysis provides the expected force trends for similar sampling strategies and allows to deduce phenomenological information about the general scooping process. In order to further instruct the community on the use of DEM tools as a solution to the sampling collection problem, two more analyses have been carried out, mainly focused on reducing the DEM computation time, which increases with a decrease in particle size. After running a set of identical simulations, where the only changing parameter is the size of the spherical particle, it is observed that the resulting force trajectories, starting from a given particle size, converge to the true trend. It is deducible that a further decrease in size yields negligible improvements in the accuracy, while it sensibly increases computation time. A final analysis aims at discussing limitations of approximating bulk material particles having a complex shape, e. g. ice fragments, with spheres, by comparing force trends resulting in the two cases for the same simulation scenario.

Catanoso, Damiana↗
DOE OSTI2022

Nb3Sn Coating of Twin Axis Cavity for SRF Applications

The twin axis cavity with two identical accelerating beams has been proposed for Energy recovery linac (ERL) applications. Nb3Sn is a superconducting material with a higher critical temperature and a higher critical field as compared to Nb, which promises a lower operating cost due to higher quality factors. Two niobium twin axis cavities fabricated at JLab and were proposed to be coated with Nb3Sn. Due to their more complex geometry, the typical coating process used for basic elliptical cavities needs to be improved to coat these cavities. This devel-opment advances the current coating system at JLab for coating complex cavities. Two twin axis cavities were coated recently for the first time. This contribution dis-cusses initial results from coating of twin axis cavities, RF testing and witness sample analysis with an overview of the current challenges towards high performance Nb3Sn coated twin axis cavities.

Tiskumara, J.↗
NASA NTRS1999

Search for Pulsations from a Nearby Millisecond Pulsar and Wasilewski 49: Mirror for a Hidden Seyfert 1 Nucleus

Five studies are reported in this final report. The recently discovered 5.3 ms pulsar J1012+5307 at a distance of 520 pc is in an area of the sky which is particularly deficient in absorbing gas. The column density along the line of sight is less than 7.5 x 10(exp 19) CM(exp -2) which facilitates soft X-ray observations. Halpern reported a possible ROSAT Position Sensitive Proportional Counter (PSPC) detection of the pulsar in a serendipitous, off-axis observation. We have now confirmed the X-ray emission of PSR J1012+,5307 in a 23 ksec observation with the ROSAT High Resolution Imager (HRI). A point source is detected within 3" of the radio position. Its count rate of 1.6 +/- 0.3 x 10(exp -3) s(exp -1) corresponds to an unabsorbed 0. 1-2.4 keV flux of 6.4 x 10(exp -14) ergs cm(exp -2) s(exp -1), similar to that reported previously. This counts-to-flux conversion is valid for N(sub H) = 5 x 10(exp 19) cm(exp -2), and either a power-law spectrum of photon index 2.5 or a blackbody of kT = 0.1 keV. The implied X-ray luminosity of 2.0 x 10(exp 30) ergs s(exp -1) is 5 X 10(exp -4) of the pulsar's spin-down power dot-E, and similar to that of the nearest millisecond pulsar J0437-4715, which is nearly a twin of J1012+5307 in P dot-E. We subjected the 37 photons (and 13 background counts) within the source region to a pulsar search, but no evidence for pulsation was found. The pulsar apparently emits over a large fraction of its rotation cycle, and the absence of sharp modulation can be taken as evidence for surface thermal emission, as is favored for PSR J0437-4715, rather than magnetospheric X-ray emission which is apparent in the sharp pulses of the much more energetic millisecond pulsar B1821-24. A further test of this interpretation will be made with a longer ROSAT observation, which will increase the number of photons collected by a factor of 5, and permit a more sensitive examination of the light curve for modulation due to emission from heated polar caps. If found, such modulation will be further evidence that surface reheating by the impact of particles accelerated along open fiel;d lines operates in these approximately 10(exp -9) year old pulsars. In a second study, a new AM Her star serendipitously in a 25 day observation was detected with the EUVE satellite. A coherent period of 85.82 min is present in the EUVE Deep Survey imager light curve of this source. A spectroscopic optical identification is made with a 19th magnitude blue star that has H and He emission lines, and broad cyclotron humps typical of a magnetic cataclysmic variable. A lower limit to the polar magnetic field of 50 MG is estimated from the spacing of the cyclotron harmonics. EUVE J0425.6-5714 is also detected in archival ROSAT HRI observations spanning two months, and its stable and highly structured light curve permits us to fit a coherent ephemeris linking the ROSAT and EUVE data over a 1.3 yr gap. The derived period is 85.82107 +/- 0.00020 min, and the ephemeris should be accurate to 0.1 cycles until the year 2005. A narrow but partial X-ray eclipse suggests that this object belongs to the group of AM Her stars whose viewing geometry is such that the accretion stream periodically occults the soft X-ray emitting accretion spot on the surface of the white dwarf. A non-detection of hard X-rays from ASCA observations that are contemporaneous with the ROSAT HRI shows that the soft X-rays must dominate by at least an order of magnitude, which is consistent with a known trend among AM Her stars with large magnetic field. This object should not be confused with the Seyfert galaxy IH 0419-577 (= LB 1727), another X-ray/EUV source which lies only 3'95 away, and was the principal target of these monitoring observations. In a third report; the identity of the persistent high-energy (> 100 MeV) gamma-ray sources in the Galaxy, still largely a mystery is investigated. The second installment of the EGRET (2EG) lists a total of 128 sources, of which 51 are likely or possibly identified with AGNs, five with rotation-powered pulsars, and one is the LMC. There are 71 unidentified sources, of which 33, or almost half, lie in the narrow band of absolute value of b < 10 degrees along the Galactic plane. This excess of low-latitude sources must, therefore, constitute a Galactic population that is either similar to the already identified gamma-ray pulsars, or an entirely new class of gamma-ray emitters associated with the disk population. We are continuing our program, begun in A06, that is aimed at intermediate-latitude sources, arguing that X-ray detection of them is the most plausible method of identifying the Galactic population. The sources at high latitude must statistically be mostly AGNS, and are more straightforwardly identified through radio and optical means. For the fourth report, identification of the EGRET Source the identity of the persistent, high-energy gamma-ray sources in the Galactic plane is a mystery. The most likely scenario is a population of middle-aged pulsars, many of which could be radio quiet like Geminga. We have an ongoing program of ROSAT, VLA, and optical observations of selected EGRET error circles at intermediate Galactic latitude. For one of these fields, at (l, b) = (106 degrees, +3), our complete census of X-ray and radio sources, reveals a remarkable association between a radio shell with unique properties, and a complete census of X-ray and radio sources reveals a remarkable association between a radio shell with unique properties, and a compact X-ray source. Further observations are needed to determine whether or not this source has a hard X-ray spectrum like that of other gamma ray pulsars and, ideally, to find its pulsations. A fifth report reviews, Wasilewski 49, an interacting pair of Seyfert galaxies at z = 0.063, one of which contains a hidden Seyfert I nucleus as evidenced by broad wings on its Balmer lines.

Halpern, Jules P.↗
DOE OSTI2022

MECHANISTIC MULTIPHYSICS MODELING OF CLADDING RUPTURE IN NUCLEAR FUEL RODS DURING LOSS-OF-COOLANT ACCIDENT CONDITIONS

The Loss of Coolant Accident (LOCA) is a design basis accident that is included as part of the safety analysis of nuclear power plants. As the nuclear industry desires to increase the cycle length and discharge burnup of existing nuclear power plants they must demonstrate safe operation during a LOCA on high burnup fuel. During a LOCA transient on high burnup fuel rods, the rods may undergo a process known as fuel fragmentation, relocation, and dispersal (FFRD). To permit dispersal, the cladding encapsulating the fuel must rupture with an opening size large enough to allow the fragmented fuel particles to release. Current licensing tools used by industry and the United States Nuclear Regulatory Commission are limited in geometric fidelity and materials that can be analyzed. These simulation tools generally employ a quasi-two-dimensional (1.5D or Layered1D) or 2D-RZ axisymmetric geometric representations exclusively. While a valid approach under some instances, there are times when important physics have an asymmetric behavior in the fuel rod. Examples include fuel fragmentation, thermal-hydraulic boundary conditions, and cladding rupture, all of which are important for LOCA analysis. As industry pursues burnup extensions it must be demonstrated that fuel dispersal can be mitigated or eliminated. To do this, an understanding of the rupture opening after cladding failure is required. This work presents the development of a model for predicting the size and location of the rupture opening in failed fuel rods during LOCA conditions using advanced modeling and simulation tools. In order to supply the rupture model with appropriate boundary conditions, improvements to fuel fragmentation, axial relocation and oxidation modeling were required. First, the eXtended Finite Element Method (XFEM) is used to mechanistically predict the number of fuel fragments that form due to material strength randomization, criteria for strength randomization, mesh density, power ramping rates and irradiation effects. These predictions with associated uncertainty were compared to empirical correlations developed for UO2 verifying that they can be used with increased confidence in subsequent axial relocation analyses. Secondly, a new first-of-its-kind Layered2D computational framework was developed that provides the ability to apply azimuthally varying boundary conditions while providing discrete layers to track fuel movement during the LOCA. An existing fuel axial relocation model developed for Layered1D was extended to work within the Layered2D framework. A large sensitivity study was performed on the initial version of the model to identify modeling parameters of particular importance, with the emissivity used for radiation after blowdown being the primary source of uncertainty. Then, a simplistic approach to incorporate mechanical degradation of the cladding due to oxidation was also developed to investigate the impact of reduced cladding thickness on predictions of the time to failure of cladding tubes. It was found that the cladding will typically rupture prior to a reduction in thickness that is sufficient to impact the rupture behavior. A model was then developed for predicting cladding rupture that transfers the cladding surface temperatures, rod internal and external pressures, fast neutron fluence, and fast neutron flux from a more detailed Layered1D, Layered2D, or 2D-RZ analysis to a 3D cladding only analysis. Comparisons of the rupture model to a few experiments indicate reasonable predictions. The rupture model was then applied to two accident tolerant fuel concepts (FeCrAl and Cr-coated Zircaloy) where it predicted that both ATF concepts would have smaller rupture openings and delayed rupture times than the standard Zircaloy-4 cladding material under identical loading conditions.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

Development of a Comprehensive Two-Phase Flow Database for the Validation of NEK-2P

Three-dimensional (3-D) two-phase Computational Fluid Dynamics (CFD) codes are emerging as a powerful and potentially practical tool for applications in which detailed local flow information is needed. However, two-phase flow models and associated closure relations are not well established for CFD applications, which is partly due to the lack of high-quality validation data. The main objective of this work is to develop a comprehensive database of two-phase flows that can be used to validate two-phase CFD codes such as NEK-2P. In this project, four advanced local measurement systems, including Particle Image Velocimetry and Planar Laser-Induced Fluorescence (PIV-PLIF), high-speed imaging, x-ray densitometry, and multi-sensor conductivity probe are employed to measure the local two-phase flow parameters of both gas and liquid phases. By combining these techniques, the local void fraction, bubble velocity, interfacial area concentration, bubble frequency, liquid velocity, turbulence intensity, etc., in various two-phase flow regimes can be obtained. These local measurement techniques are first used in a 25.4 mm circular pipe test section. Seven air-water two-phase flow conditions spanning the bubbly, slug, churn-turbulent, and annular flow regimes are measured in this facility. The obtained database contains the radial profiles of both gas- and liquid-phase parameters at three axial locations along the test section. The second facility used in this work contains a 30 mm × 10 mm rectangular test section, in which three two-phase flow and two single-phase flow conditions are measured. Two-dimensional distributions of local two-phase flow parameters in the cross-sectional plane are measured in this facility at three axial locations as well. A facility featuring a 3×3 electrically heated rod bundle is also designed and being constructed in this project. This facility is specially designed for optical measurements and is expected to provide high-quality boiling data in the future. Preliminary analyses have been performed for the data taken in the 25.4 mm circular pipe. Both center-peaked and wall-peaked void fraction profiles have been observed in the data depending on the two-phase flow conditions and/or developing lengths. The 1-D drift-flux model was evaluated with the newly obtained datasets, in which both gas- and liquid-phases data were directly measured. The distribution parameter model has been optimized based on a new void-profile classification method proposed in this study. The optimized drift-flux model shows a significant improvement in predicting both gas velocity and void fraction. The measured liquid-phase turbulence was used to benchmark Sato’s turbulence model considering the bubble-induced shear stress for the three tested bubbly flows. The benchmark results showed good agreement between the PIV measurements and model predictions. In the bubbly flows tested that have low void fractions less than 3%, the effect of the bubble-induced turbulence was found not significant. However, the bubble-induced shear stress becomes important with the increase of the void fraction. The Conjugate Heat Transfer (CHT) model was developed and implemented in NEK-2P. This model allows the coupled simulation of the solid domain and two-phase fluid domain, allowing the specification of realistic boundary conditions. The CHT implementation was verified first with non-boiling simulations. The predicted temperatures in the fluid domain were shown to be identical in simulations with or without the CHT model. The CHT model was then validated through simulations of three Becker benchmark CHF tests performed under both Dryout (DO) and Departure from Nucleate Boiling (DNB) conditions. Reasonably good agreement was obtained between calculated wall temperatures and corresponding experimental data.

42 ENGINEERING↗