Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Re”

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,405 records · Page 78

Combustion of 10–100 μm aluminum droplets in detonation products gases

Here we describe a two-phase model of combustion effects in aluminized high explosive (HE) charges. It is based on: (i) a Gas Dynamic Model of the expansion of the detonation product gases and their turbulent combustion with air; and (ii) a Heterogeneous Continuum Model of aluminum (Al) droplets and their combustion with the detonation product gases. Initial conditions are based on an analytical similarity solution for a cylindrical Chapman-Jouguet (CJ) detonation propagating at the CJ detonation velocity. The CJ jump conditions are computed at the thermodynamic equilibrium state by the Cheetah code, assuming the Al droplets are inert. We assume that the Al is 10 % of the charge mass and occurs as droplets at the CJ state. Different initial droplet diameters, ranging from 10 to 100 microns, are studied. A hydrodynamic combustion model based on large Damköhler numbers is employed in this study, -3σ(1 + 0.276 $\sqrt{Re}$)/($Kd^2_w$). It has a square-root dependence on the Reynolds number (Re) and inversed-squared dependence on the droplet diameter (d w ). The burnout time (t B ) of the Al droplets has a three-halves dependence on the droplet diameter, t B ~ ($d^0_w$) 3/2 . After burnout, the detonation products act as detonation products of the HE charge with active Al. They turbulently mix with air and form a combustion layer on the outer edge of the fireball. Details of the two-phase model, initial conditions and evolution of the flow field will be described.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Effects of sub‐mm cylindrical voids on detonation performance in PBX 9501

Abstract Internal features of varying scale and geometry are always present in explosives systems. Below a critical length, dependent on the explosive, these features can operate as a driving force for energy concentration and reaction, known as hotspots. At larger length scales, internal features result in jetting and detonation wave shaping, allowing for bulk work to be done by the explosive as is seen in shape charges. To date, a large volume of work has been performed to simulate hotspot ignition and large‐scale wave shaping. However, little experimental data exists on the effects of features in intermediate length scales, 0.1 to 1 millimeter. It has been observed in many tests that these small‐scale features can influence the high explosive (HE) performance and in some cases cause substantial damage to adjacent systems. This work provides quantitative data on the effects of machined voids moderately above more typical hotspot lengths, 0.3–0.8 mm in diameter, in PBX 9501 pressed to 1.785 g cm −3 ±2.5 mg cm −3 . Streak imaging was used to visualize void collapse, jet velocity, re‐initiation time, and wave shape evolution. Delay in detonation front propagation time was found to be linearly dependent on the void diameter and jet velocity was found to be independent within the tested range and resolution. Cut‐back experiments were used to investigate wave shape distortion and evolution downstream of the void, showing consistent growth and decay shapes across all void sizes. Simulations using CTH, a hydrocode by Sandia National Laboratory, were used to investigate void collapse showing agreement with the trends of experimental results but yielded inaccurate wave shape development delay values. Both simulation and experimental results identified several re‐initiation mechanisms with jetting and subsequent double shocking of localized HE being the dominant mechanisms for the void sizes that were studied.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

$\mathrm{RADICAL}$-Pilot and $\mathrm{PMIx}$/$\mathrm{PRRTE}$: Executing Heterogeneous Workloads at Large Scale on Partitioned $\mathrm{HPC}$ Resources

Execution of heterogeneous workflows on high-performance computing (HPC) platforms present unprecedented resource management and execution coordination challenges for runtime systems. Task heterogeneity increases the complexity of resource and execution management, limiting the scalability and efficiency of workflow execution. Re-source partitioning and distribution of tasks execution over portioned re-sources promises to address those problems but we lack an experimental evaluation of its performance at scale. Here this paper provides a performance evaluation of the Process Management Interface for Exascale (PMIx) and its reference implementation PRRTE on the leadership-class HPC plat-form Summit, when integrated into a pilot-based runtime system called RADICAL-Pilot. We partition resources across multiple PRRTE Distributed Virtual Machine (DVM) environments, responsible for launching tasks via the PMIx interface. We experimentally measure the work-load execution performance in terms of task scheduling/launching rate and distribution of DVM task placement times, DVM startup and termination overheads on the Summit leadership-class HPC platform. Integrated solution with PMIx/PRRTE enables using an abstracted, standardized set of interfaces for orchestrating the launch process, dynamic process management and monitoring capabilities. It extends scaling capabilities allowing to overcome a limitation of other launching mechanisms (e.g., JSM/LSF). Explored different DVM setup configurations provide insights on DVM performance and a layout to leverage it. Our experimental results show that heterogeneous workload of 65,500 tasks on 2048 nodes, and partitioned across 32 DVMs, runs steady with resource utilization not lower than 52%. While having less concurrently executed tasks resource utilization is able to reach up to 85%, based on results of heterogeneous workload of 8200 tasks on 256 nodes and 2 DVMs.

97 MATHEMATICS AND COMPUTING↗

Coarse-grained fixed-point tensor networks and holographic reflected entropy in 3D gravity

We use the framework of fixed-point BCFT tensor networks to present a microscopic CFT derivation of the correspondence between reflected entropy (RE) and entanglement wedge cross section (EW) in AdS 3 /CFT 2 , for both bipartite and multipartite settings. These fixed-point tensor networks, obtained by triangulating Euclidean CFT path integrals, allow us to explicitly construct the canonical purification via cutting-and-gluing CFT path integrals. Employing modular flow in the large-c limit, we demonstrate that these intrinsic CFT manipulations reproduce bulk geometric prescriptions, without assuming the AdS/CFT dictionary. The emergence of bulk geometry is traced to coarse-graining over heavy states in the large-c limit. Universal coarse-grained BCFT data for compact 2D CFTs, through the relation to Liouville theory with ZZ boundary conditions, yields hyperbolic geometry on the Cauchy slice. The corresponding averaged replica partition functions reproduce all candidate EWs, arising from different averaging patterns, with the dominant one providing the correct RE and EW. In this way, many heuristic tensor-network intuitions in toy models are made precise and established directly from intrinsic CFT data.

AdS-CFT correspondence↗

Quark sivers function at small x: spin-dependent odderon and the sub-eikonal evolution

We apply the formalism developed earlier for studying transverse momentum dependent parton distribution functions (TMDs) at small Bjorken x to construct the small- x asymptotics of the quark Sivers function. First, we explicitly construct the complete fundamental “polarized Wilson line” operator to sub-sub-eikonal order: this object can be used to study a variety of quark TMDs at small x . We then express the quark Sivers function in terms of dipole scattering amplitudes containing various components of the “polarized Wilson line” and show that the dominant (eikonal) term which contributes to the quark Sivers function at small x is the spin-dependent odderon, confirming the re- cent results of Dong, Zheng and Zhou. Our conclusion is also similar to the case of the gluon Sivers function derived by Boer, Echevarria, Mulders and Zhou. We also analyze the sub-eikonal corrections to the quark Sivers function using the constructed “polarized Wilson line” operator. We derive new small- x evolution equations re-summing double-logarithmic powers of α s ln 2 (1 /x ) with α s the strong coupling constant. We solve the corresponding novel evolution equations in the large- N c limit, obtaining a sub-eikonal correction to the spin-dependent odderon contribution. We conclude that the quark Sivers function at small x receives contributions from two terms and is given by ${f}_{1T}^{\perp q}\left(x,{k}_T^2\right)={C}_O\left(x,{k}_T^2\right)\frac{1}{x}+{C}_1\left({k}_T^2\right){\left(\frac{1}{x}\right)}^0+\cdots$ with the function C O ( x, ${k}_T^2$) varying slowly with x and the ellipsis denoting the subasymptotic and sub-sub-eikonal (order- x ) corrections.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Strong constraints on the b → sγ photon polarisation from B0 → K*0e+e-decays

An angular analysis of the B 0 → K* 0 e + e - decay is performed using a data sample corresponding to an integrated luminosity of 9 fb - 1 of pp collisions collected with the LHCb experiment. The analysis is conducted in the very low dielectron mass squared ( q 2 ) interval between 0 . 0008 and 0 . 257 GeV 2 , where the rate is dominated by the B 0 → K * 0 γ transition with a virtual photon. The fraction of longitudinal polarisation of the K * 0 meson, F L , is measured to be F L = (4 . 4 ± 2 . 6 ± 1 . 4)%, where the first uncertainty is statistical and the second systematic. The \( {A}_{\mathrm{T}}^{\mathrm{Re}} \) observable, which is related to the lepton forward-backward asymmetry, is measured to be \( {A}_{\mathrm{T}}^{\mathrm{Re}} \) = - 0 . 06 ± 0 . 08 ± 0 . 02. The \( {A}_{\mathrm{T}}^{(2)} \) and \( {A}_{\mathrm{T}}^{\mathrm{Im}} \) transverse asymmetries, which are sensitive to the virtual photon polarisation, are found to be \( {A}_{\mathrm{T}}^{(2)} \) = 0 . 11 ± 0 . 10 ± 0 . 02 and \( {A}_{\mathrm{T}}^{\mathrm{Im}} \) = 0 . 02 ± 0 . 10 ± 0 . 01. The results are consistent with Standard Model predictions and provide the world’s best constraint on the b → sγ photon polarisation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Growth of organized flow coherent motions within a single-stream shear layer: 4D-PTV measurements

Abstract This study investigates the evolution of a single-stream shear layer (SSSL) originating from a wall boundary layer past a backward-facing step. Utilizing a time-resolved 3D-Particle Tracking Velocimetry (4D-PTV) technique, we track the trajectories of fluorescent particles to gain insight into the flow characteristics of the SSSL. A compact water tunnel facility ( $$\textrm{Re}_\tau =1\,240$$ Re τ = 1 240 ) is fabricated to obtain an SSSL with a perpendicular slow entrainment stream past the separation edge. A hybrid interpolation approach that combines ensemble binning and Gaussian weighting is implemented to derive minimally filtered mean and instantaneous lower- and higher-order flow field parameters. Spanwise-dominant coherent motion accompanied by finer flow scales is observed to grow due to flow entrainment through “nibbling” actions of small-scale vortices, “engulfing” by large-scale vortices, and vortex pairing events. Furthermore, the non-zero-speed stream edge grows relatively faster than the zero-speed stream edge, showing a strong asymmetry in mixing composition across a mixing layer. The SSSL reaches self-similarity at a streamwise distance of $$\approx 55\,\theta _{0}$$ ≈ 55 θ 0 , where $$\theta _0$$ θ 0 is the initial momentum thickness from the separation edge, i.e., considerably shorter than reported in previous studies. A literature comparison of growth rate parameters raises intriguing questions regarding a potential inclusive growth scaling unifying the free shear layers. A turbulent kinetic energy (TKE) budget analysis reveals a negative production region immediately downstream of the separation edge attributed to a large positive streamwise gradient of streamwise velocity. In the self-similar region, the phase-averaged flow mapping demonstrates a larger concentration of turbulence production rate around the outer edges of spanwise vortices, specifically at the intersection of braids and vortices. Furthermore, a spatial separation exists in the regions of peak production and dissipation rates within the vortex core region favoring dissipation. The braids exhibit a larger concentration of turbulence diffusion rates, indicating their function as a conduit for exchanging turbulence between neighboring coherent motions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Revisiting a Drag Partition Model For Canopy-Like Roughness Elements

Turbulent flows over a large surface area (S) covered by n obstacles experience an overall drag due to the presence of the ground and the protruding obstacles into the flow. The drag partition between the roughness obstacles and the ground is analyzed using an analytical model proposed by Raupach (Boundary-Layer Meteorol 60:375-395, 1992) and is hereafter referred to as R92. The R92 is based on the premise that the wake behind an isolated roughness element can be described by a shelter area A and a shelter volume V. The individual sizes of A and V without any interference from other obstacles can be determined from scaling analysis for the spread of wakes. To upscale from an individual roughness element to n/S elements where wakes may interact, R92 adopted a background stress re-normalizing instead of reducing A or V with each element addition. This work demonstrates that R92’s approach results in a linear background stress reduction in A and V only when the ratio of n/S is small, due to a low probability of wake interactions. This probabilistic nature suggests that up-scaling from individual to multiple roughness elements can be re-formulated using stochastic averaging methods proposed here. The two approaches are shown to recover R92 under plausible conditions. An alternative scaling for the shelter volume is also proposed here using thermodynamic arguments of work and dissipation though the final outcome remains similar to R92. Comparisons between R92 and available data spanning more than two decades after R92 on blocks and vegetation-like roughness elements confirm the practical utility of R92. The agreement between R92 and this updated databases of experiments and simulations confirm the potential use of R92 in large-scale models provided that the relevant parameters accommodate certain features of the roughness element type (cube versus vegetation-like) and, to a lesser extent, their configuration throughout S. Last, a comparison between R92 and models based on first-order closure principles with constant mixing length suggests that R92 can outperform such models when evaluated across a wide range of roughness densities.

canopy turbulence↗

Optimizing power system restoration with damaged communications

Utility procedures for power system blackstart and restoration typically assume that energization decisions can be reliably communicated across the grid. In reality, the communications and control network would likely also be affected in power outages, such as those caused by extreme weather events or cyber-attacks. This paper studies the effect of damage to the power system communications and control infrastructure on restoration operations following a blackout. We model the communications infrastructure as a graph, overlaying the power grid, and imposing the requirement that every energized element in the power grid be observable from a control center. We expand on a specialized branch-and-bound algorithm from the literature to optimize the restoration process and devise an initialization heuristic and a rounding heuristic to improve solution speed. We perform numerical experiments on synthetic systems for Illinois and Texas with outages based on a solar flare or hurricane. We compare the results of our specialized branch-and-bound algorithm to the results from (i) the initialization heuristic alone, (ii) a variation of this heuristic that we use as a baseline, and (iii) the restoration optimization for the power system without communications constraints. Here, we find that damage to the communications infrastructure significantly increases the time required to re-energize the grid. Moreover, by simultaneously optimizing communications repairs and grid energization decisions, we are able to re-energize the grid significantly faster than if communications repairs and energization decisions were made independently or with partial coordination, motivating improvements to current industry practice.

97 MATHEMATICS AND COMPUTING↗

Slip localization behavior at triple junctions in nickel-base superalloys

Incipient slip localization in the vicinity of hundreds of grain boundary triple junctions (TJs) in a lightly deformed nickel-base superalloy IN718 is studied using a combination of three-dimensional (3D) crystal plasticity finite element (CPFE) modeling, high resolution digital image correlation (HR-DIC) and 3D electron back-scatter diffraction tomography (3D EBSD). A 3D reconstruction method enables identification of thousands of TJs and correspondence of any observed slip bands with their originating TJ lines below the specimen surface. Here, we present a large-scale CPFE model of the experimental 3D microstructure composed of high-fidelity representation of the TJ lines and the boundaries and interiors of the parent grains and use it to calculate the local micromechanical response and slip activity of all TJs at the onset of macroscopic yielding. Statistical analysis of the calculated quantities reveal TJs develop larger stress concentration and grain-average re-orientation than grain interiors and grain boundaries, however no substantial differences in cumulative slip were found among these microstructural regions. We find that TJs with observed slip bands generate lower grain-average re-orientation, fewer active slip systems, and more localized slip on a single system than those without. The distinctions in the reorientation and slip activity are stronger in TJs that experience more intense slip.

36 MATERIALS SCIENCE↗

Degradation of electrical resistivity of tungsten following shielded neutron irradiation

A major challenge for heat transfer in nuclear materials is to ensure thermal mobility after high amounts of neutron irradiation. Tungsten is widely selected as a heat transfer material in fusion reactors. In metals, thermal conductivity is dominated by electrons’ ability to transfer energy. Neutron irradiation generates point defects, clusters, and solid transmutation (e.g.rhenium and osmium in tungsten), which inhibit electron motion. The purpose of this work is to quantify the irradiation-induced change in electron mobility and deconvolute transmutation and microstructural effects on observed changes to electron mobility. Single and polycrystalline tungsten were fast neutron irradiated in the High Flux Isotope Reactor at Oak Ridge National Laboratory to doses between 0.2 and 0.7 displacements per atom (dpa) and temperatures from 500 °C to 1000 °C. Grain growth was observed in all samples. Microstructure and transmutation were quantified. The geometric orientation of samples with elongated grains has been shown to affect electrical resistivity. A mathematical model was developed and used to deconvolute solid-solution transmutation, grain, and temperature-dependent lattice effects on resistivity. At ~0.4 dpa at ~590 °C, the combined resistivity degradation due to voids, vacancies, interstitials, and dislocations is estimated to be greater than the contribution from solid solution Re transmutation, which is greater than the contribution from grain boundaries. At doses of ~0.7 dpa at ~750 °C, solid solution Re contributions are greater than all other effects combined. As a result, this work establishes a basis to predict the effects of irradiation temperature and transmutation on thermal properties of tungsten and highlights the importance of irradiation temperature.

36 MATERIALS SCIENCE↗

On dysprosium utilisation in multi-main-phase Nd–Dy–Fe–B magnets with core–shell microstructures

The development of high-performance Nd–Dy–Fe–B magnets that minimise the consumption of the scarce rare earth (RE) element Dy remains a major global scientific and technological quest. Here, we designed an alloy microstructure comprising of a uniform Dy-lean core–Dy-rich shell in a series of multi-main-phase (MMP) Nd–Dy–Fe–B magnets. The resulting MMP Dy1 and Dy3 magnets with an overall Dy level of 1 and 3 wt.% possessed values of 0.48 and 0.29 T/wt.% of coercivity increment per unit weight percentage of the Dy addition, respectively. Most importantly, the resulting MMP Dy3 magnet exhibited a high coercivity (2.38 T), an excellent thermal stability of the coercivity (|β| = 0.531%/°C), a high squareness factor (> 95%), all with little diminishment in the remanent magnetisation (1.35 T) and maximum energy product (43.6 MGOe). These properties are superior to the currently available sintered Nd–Dy–Fe–B magnets which utilise higher levels of Dy of 5 wt.%. Via magnetic and multi-scale microstructural characterisation experiments and micromagnetic simulations, the formation of the Dy-lean core–Dy-rich shell microstructure is rationalised via solid-state-diffusion and solution reprecipitation during liquid-phase sintering. The Dy-lean core–Dy-rich shell microstructure and the non-ferromagnetic low-Fe RE-rich grain boundary phase led to the synergistic magnetic performance. This is significant in the context of the MMP Nd–Dy–Fe–B magnets being applied to large-scale production. The present work establishes a pathway for the more sustainable utilisation of Dy in permanent magnets via formation of a uniform core–shell microstructure.

36 MATERIALS SCIENCE↗

Distinguishing localization from chaos: Challenges in finite-size systems

Highlights: • Provides an overview of the current understanding of the many-body localization phase transition. • Discusses the subtleties of finite-size scaling near this transition • Assesses the implications of these subtleties for numerical studies of spin chains. • Explores scaling of diagnostics in models with known localization transitions. • Presents suggestions for future numerical work. We re-examine attempts to study the many-body localization transition using measures that are physically natural on the ergodic/quantum chaotic regime of the phase diagram. Using simple scaling arguments and an analysis of various models for which rigorous results are available, we find that these measures can be particularly adversely affected by the strong finite-size effects observed in nearly all numerical studies of many-body localization. This severely impacts their utility in probing the transition and the localized phase. In light of this analysis, we discuss a recent study (Šuntajs et al., 2020) of the behaviour of the Thouless energy and level repulsion in disordered spin chains, and its implications for the question of whether MBL is a true phase of matter.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Multiscale Effects Masked the Impact of the COVID-19 Pandemic on Electricity Demand in the United States

Shelter-in-place orders and business closures related to COVID-19 changed the hourly profile of electricity demand and created an unprecedented source of uncertainty for the grid. The potential for continued shifts in electricity profiles has implications for electricity sector investment and operating decisions that maintain reserve margins and provide grid reliability. This study reveals that understanding this uncertainty requires an understanding of the underlying drivers at the customer-class scale. This paper utilizes three datasets to compare the impacts of COVID-19 on electricity consumption across a range of spatiotemporal and customer scales. At the utility/customer-class scale, COVID-19-induced shutdowns in the spring of 2020 shifted weekday residential load profiles to resemble weekend profiles from previous years. Total commercial loads declined, but the commercial diurnal load profile was unchanged. With only total loads available at the balancing authority scale, the apparent impact of COVID-19 was smaller during the summer due in part to phased re-opening and spatial variability in re-opening, but there were still clear variations once total loads were broken down zonally. Monthly data at the state scale showed an increase in state-level residential electricity sales, a decrease in commercial sales, and a small net decrease in total sales in most states from April-August 2020. Analyses that focus on total load or a single scale may miss important changes that become apparent when the load is broken down regionally or by customer class.

COVID-19, electricity demand, multiscale, Commonwe↗

Experimental demonstration of a 5 kW th granular-flow reactor for solar thermochemical energy storage with aluminum-doped calcium manganite particles

A two-step cycle was discussed for solar thermochemical energy storage based on aluminum-doped calcium manganite reduction/oxidation reactions for direct integration into Air Brayton cycles. The two steps encompassed 1) the storage of concentrated solar direct irradiation via the thermal reduction of aluminum-doped calcium manganite and 2) the delivery of heat to an Air-Brayton cycle via re-oxidation of oxygen-deficient aluminum-doped calcium manganite. The re-oxidized aluminum-doped calcium manganite was fed back to the first step to complete the cycle. A 5 kWth solar thermochemical reactor operating under vacuum was fabricated and tested to examine the first cycle reduction step. Reactor operating conditions and high-flux solar simulator control were tuned for continuous reactor operation with particle temperatures > 1073 K. Continuous operation was achieved using intermittent, dense granular flows. A maximum absorption efficiency of 64.7% was shown, accounting for both sensible and chemical heat storage.

25 ENERGY STORAGE↗

Heat Transfer Coefficients of Additively Manufactured Tubes with Internal Pin Fins for Supercritical Carbon Dioxide Cycle Recuperators

This paper describes the measurement of convective heat transfer coefficients and friction factors for sCO 2 flowing in additively manufactured tubes with internal pin fins at the US DOE’s National Energy Technology Laboratory in Morgantown, WV. The measurement procedures were validated by conducting benchmark tests with smooth stainless-steel tube and comparing the results with published correlations for Nusselt number (Nu) and friction factor. Over Reynolds numbers (ReD) ranging from 5 × 10 4 to 2.5 × 10 5 , measured Nu was within 5% of the Dittus-Boelter correlation and measured friction factors were within 5% of the McAdams correlation for smooth tube flow. The candidate pin fin patterned pipes were additively manufactured (AM) at the Oak Ridge National Laboratory. The pins were circular or elliptical in cross-section and were printed at a 30⁰ angle relative to the inner wall (to meet AM constraints). The pin arrangement was helical to promote enhanced heat transfer due to swirl flow. Pin length to diameter aspect ratio was 1.33, 2, and 8, while the pin diameter to tube diameter ratio was 0.188, 0.125, and 0.063. Tests were performed for ReD varying from 6.9 × 10 4 to 2.2 × 10 5 and at conditions equivalent to the low pressure (LP) outlet (8.69 MPa, 361 K) and the high pressure (HP) inlet (20.7 MPa, 350 K) of the low temperature recuperator (LTR) in an indirect sCO 2 cycle. The Wilson plot technique was utilized to measure the bulk heat transfer coefficients. For the best performing design (tube A, pin length to tube diameter ratio: 1.33, pin diameter to tube diameter ratio: 0.19), the local heat transfer coefficient increased by 136% relative to the Dittus-Boelter correlation at the LTR low pressure outlet and 194% at the LTR high pressure inlet. These correspond to a 282% and a 271% increase in the product of the heat transfer coefficient and surface area (adjusted for fin efficiency) product, respectively. Large pressure drops across the test articles were observed. For Tube Design A, the average friction factor, across the range of ReD considered, was significantly larger than the McAdams correlation at both the LTR LP outlet and the LTR HP inlet. A thermal performance factor was utilized to express the ratio of material required to build a finned heat exchanger relative to a finless heat exchanger with the same heat duty and pumping power. Tube Design A was estimated to decrease the required heat exchanger material by 13%.

Searle, Matthew↗

Critical heat flux on zircaloy and accident tolerant fuel cladding under prototypical conditions of pressurized and boiling water reactors

First of a kind high spatial and temporal resolution temperature measurements of Pressurized Water Reactor and Boiling Water Reactor simulated fuel pins (high heat flux cosine profile heaters) in annular rod type geometry during nucleate boiling, departure from nucleate boiling and re-wetting has been performed. Distributed fiber optic measurements of the cladding axial temperature with spatial resolutions down to 2.5 mm axially at up to three equally spaced azimuthal locations per rod have demonstrated the similarities and differences of the different accident tolerant fuel cladding materials. These experiments show the performance of accident tolerant cladding with a precise measurement of the location and evolution of dryout followed by the re-wetting phenomena (see supplementary time history of a critical heat flux event). Experimental parameters varied coolant mass fluxes from 1695 to 2712 kg/m 2 s, pressures from 10 to 20 MPa, and inlet subcooling from 10 to 55 °C. It was found that the coatings and different material have little effect on the value of the critical heat flux. Any discrepancies between the local critical heat flux verified among the cladding material are attributed to different departure from nucleate boiling phenomena captured during the experiment with the high-resolution optical fiber temperature data. Substantial differences in the performance of the different cladding material to survive a critical heat flux event were however observed by post-test X-ray and optical inspection. External evident damage and internal cracks to the simulated fuel pin having the bare zircaloy cladding were observed in regions associated with the occurrence of the critical heat flux for longer periods. Such damages were associated with zirconium oxidation, and not noticed on the Cr coated zircaloy or FeCrAl claddings. These results indicate that the use of FeCrAl or a Cr coated zircaloy as cladding in Pressurized Water Reactor and Boiling Water Reactor may increase the survivability of the cladding material for short transients.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗