Azimuthal Temperature Variations Research Articles

Three-dimensional flow structures in turbulent Rayleigh–Bénard convection at low Prandtl number Pr = 0.03

In this paper we report on an experimental study focusing on the manifestation and dynamics of the large-scale circulation (LSC) in turbulent liquid metal convection. The experiments are performed inside a cylinder of aspect ratio $\varGamma = 0.5$ filled with the ternary alloy GaInSn, which has a Prandtl number of $Pr = 0.03$ . The large-scale flow structures are classified and characterized at Rayleigh numbers of ${Ra} = 9.33 \times 10^6, 5.31 \times 10^7$ and $6.02 \times 10^8$ by means of the contactless inductive flow tomography which enables the full reconstruction of the three-dimensional (3-D) flow structures in the entire convection cell. This is complemented with the multi-thermal-probe method for capturing the azimuthal temperature variation induced by the LSC at the sidewall. We use proper orthogonal decomposition (POD) to identify the dominating modes of the turbulent convection. The analysis reveals that a single-roll structure of the LSC alternates in short succession with double-roll structures or a three-roll structure. This is accompanied by dramatic fluctuations of the Reynolds number, whose instantaneous values can deviate by more than 50 % from the time-average value. No coherent oscillations are observed, whereas a correlation analysis indicates a residual contribution of the torsion and sloshing modes. Results of the POD analysis suggest a stabilization of the single-roll LSC with increasing $Ra$ at the expense of flow structures with multiple rolls. Moreover, the relative lifetime of all identified flow states, measured in units of free-fall times, increases with rising $Ra$ .

Azimuthal temperature variations in ISO-Oph 2 from multifrequency ALMA observations

ABSTRACT Environmental effects, such as stellar fly-bys and external irradiation, are thought to affect the evolution of protoplanetary discs in clustered star formation. Previous Atacama Large Millimetre/submillimetre Array (ALMA) images at 225 GHz of the ISO-Oph 2 binary revealed a peculiar morphology in the disc of the primary, perhaps due to a possible fly-by with the secondary. Here, we report on new ALMA continuum observations of this system at 97.5, 145, and 405 GHz, which reveal strong morphological variations. Multifrequency positional alignment allows us to interpret these spectral variations in terms of underlying physical conditions. ISO-Oph 2A is remarkably offset from the centroid of its ring, at all frequencies, and the disc is lopsided, pointing at gravitational interactions. However, the dust temperature also varies in azimuth, with two peaks whose direction connects with HD 147889, the earliest-type star in the Ophiuchus complex, suggesting that it is the dominant heat source. The stellar environment of ISO-Oph 2 appears to drive both its density structure and its thermal balance.

SO and SiS Emission Tracing an Embedded Planet and Compact 12CO and 13CO Counterparts in the HD 169142 Disk

Planets form in dusty, gas-rich disks around young stars, while at the same time, the planet formation process alters the physical and chemical structure of the disk itself. Embedded planets will locally heat the disk and sublimate volatile-rich ices, or in extreme cases, result in shocks that sputter heavy atoms such as Si from dust grains. This should cause chemical asymmetries detectable in molecular gas observations. Using high-angular-resolution ALMA archival data of the HD 169142 disk, we identify compact SO J = 88 − 77 and SiS J = 19 − 18 emission coincident with the position of a ∼ 2 M Jup planet seen as a localized, Keplerian NIR feature within a gas-depleted, annular dust gap at ≈38 au. The SiS emission is located along an azimuthal arc and has a morphology similar to that of a known 12CO kinematic excess. This is the first tentative detection of SiS emission in a protoplanetary disk and suggests that the planet is driving sufficiently strong shocks to produce gas-phase SiS. We also report the discovery of compact 12CO and 13CO J = 3 − 2 emission coincident with the planet location. Taken together, a planet-driven outflow provides the best explanation for the properties of the observed chemical asymmetries. We also resolve a bright, azimuthally asymmetric SO ring at ≈24 au. While most of this SO emission originates from ice sublimation, its asymmetric distribution implies azimuthal temperature variations driven by a misaligned inner disk or planet–disk interactions. Overall, the HD 169142 disk shows several distinct chemical signatures related to giant planet formation and presents a powerful template for future searches of planet-related chemical asymmetries in protoplanetary disks.

A physically motivated “charge-exchange method” for measuring electron temperatures within H ii regions

Aims. Temperature uncertainties plague our understanding of abundance variations within the interstellar medium. Using the PHANGS-MUSE large program, we develop and apply a new technique to model the strong emission lines arising from H ii regions in 19 nearby spiral galaxies at ~50 pc resolution and infer electron temperatures for the nebulae. Methods. Due to the charge-exchange coupling of the ionization fraction of the atomic oxygen to that of hydrogen, the emissivity of the observed [O i]λ6300/Hα line ratio can be modeled as a function of the gas phase oxygen abundance (O/H), ionization fraction (fion), and electron temperature (Te). We measure O/H using a strong-line metallicity calibration and identify a correlation between fion and [S iii]λ9069/[S ii]λ6716,6730, tracing ionization parameter variations. Results. We solve for Te and test the method by reproducing direct measurements of Te([N ii]λ5755) based on auroral line detections to within ~600 K. We apply this“charge-exchange method” of calculating Te to 4129 H ii regions across 19 PHANGS-MUSE galaxies. We uncover radial temperature gradients, increased homogeneity on small scales, and azimuthal temperature variations in the disks that correspond to established abundance patterns. This new technique for measuring electron temperatures leverages the growing availability of optical integral field unit spectroscopic maps across galaxy samples, increasing the statistics available compared to direct auroral line detections.

Analysis of fuel rod behavior during loss-of-coolant accidents using the BISON code: Cladding modeling developments and simulation of separate-effects experiments

BISON is a finite-element based, multidimensional fuel performance code developed at Idaho National Laboratory. In this paper, we present BISON modeling developments for Zircaloy nuclear fuel cladding behavior under loss-of-coolant accident (LOCA) conditions and code validation to separate-effects experiments. Code developments include models for cladding high-temperature creep, crystallographic phase transition, high-temperature steam oxidation and failure due to burst. Code validation is performed against a substantial number of experimental tests for ballooning and burst behavior of Zircaloy claddings under LOCA conditions. Validation calculations are completed using 2D axisymmetric models. Additionally, we present a 3D calculation that demonstrates simulation of the cladding response in presence of azimuthal temperature variations. Calculated results are systematically compared to experimental data of cladding temperature, inner pressure and time at burst failure. Several of the presented BISON validation cases were developed in the framework of the IAEA Coordinated Research Project FUMAC. BISON modeling developments for UO2 fuel behavior under LOCA conditions and code validation to integral fuel rod experiments are presented in a complementary paper.

Deep neural network based prediction of burst parameters for Zircaloy-4 fuel cladding during loss-of-coolant accident

BackgroundUnderstanding the behaviour of nuclear fuel claddings by conducting burst test on single cladding tube under simulated loss-of-coolant accident conditions and developing theoretical cum empirical predictive computer codes have been the focus of several investigations. The developed burst criterion (a) assumes symmetrical deformation of cladding tube in contrast to experimental observation (b) interpolates the properties of Zircaloy-4 cladding in mixed α+ β phase (c) does not account for azimuthal temperature variations. In order to overcome all these drawbacks of burst criterion, it is reasoned that artificial intelligence technique may be a better option to predict the burst parameters. MethodsArtificial neural network models based on feedforward backpropagation algorithm with logsig transfer function are developed. ResultsNeural network architecture of 2-4-4-3, that is model with two hidden layers having four nodes in each layer is found to be the most suitable. The mean, maximum, and minimum prediction errors for this optimised model are 0.82%, 19.62%, and 0.004%, respectively. ConclusionThe burst stress, burst temperature, and burst strain obtained from burst criterion have average deviation of 19%, 12%, and 53% respectively whereas the developed neural network model predicted these parameters with average deviation of 6%, 2%, and 8%, respectively.

Large-scale flow and Reynolds numbers in the presence of boiling in locally heated turbulent convection

We report on an experimental study of the large-scale flow (LSF) and Reynolds numbers in turbulent convection in a cylindrical sample with height equal to its diameter and heated locally around the center of its bottom plate (locally heated convection). The sample size and shape are the same as those of Narezo Guzman et al. [D. Narezo Guzman, J. Fluid Mech. 787, 331 (2015)JFLSA70022-112010.1017/jfm.2015.701; D. Narezo Guzman, J. Fluid Mech. 795, 60 (2016)JFLSA70022-112010.1017/jfm.2016.178]. Measurements are made at a nearly constant Rayleigh number as a function of the mean temperature, both in the presence of controlled boiling (two-phase flow) and for the superheated fluid (one-phase flow). Superheat values Tb-Ton up to about 11 K (Tb is the bottom-plate temperature and Ton is the lowest Tb at which boiling is observed) are used. The LSF is less organized than it is in (uniformly heated) Rayleigh-Benard convection (RBC), where it takes the form of a single convection roll. Large-scale-flow-induced sinusoidal azimuthal temperature variations (like those found for RBC) could be detected only in the lower portion of the sample, indicating a less organized flow in the upper portions. Reynolds numbers are determined using the elliptic model (EM) of He and Zhang [G.-W. He and J.-B. Zhang, Phys. Rev. E 73, 055303(R) (2006)PLEEE81539-375510.1103/PhysRevE.73.055303]. We found that for our system the EM is applicable over a wide range of space and time displacements, as long as these displacements are within the inertial range of the temporal and spatial spectrum. At three locations in the sample the results show that the vertical mean-flow velocity component is reduced while the fluctuation velocity is enhanced by the bubbles of the two-phase flow. Enhancements of velocity fluctuations up to about 60% are found at the largest superheat values. Local temperature measurements within the sample reveal temperature oscillations that also used to determine a Reynolds number. These results are generally consistent with the mean-flow EM results and show a two-phase-flow enhancement of up to about 30%.

Proof-of-principle of high-fidelity coupled CRUD deposition and cycle depletion simulation

A multiphysics framework for the high-fidelity simulation of CRUD deposition is developed to better understand the coupled physics and their respective feedback mechanisms. This framework includes the primary physics of lattice depletion, computational fluid dynamics, and CRUD chemistry. The three physics are coupled together via the operator-splitting technique, where predictor–corrector and fixed-point iteration schemes are utilized to converge the nonlinear solution. High-fidelity simulations may provide a means to predict and assess potential operating issues, including CRUD induced power shift and CRUD induced localized corrosion, known as CIPS and CILC, respectively. As a proof-of-principle, a coupled 500-day cycle depletion simulation of a pressurized water reactor fuel pin cell was performed using the coupled code suite; a burnup of 31MWd/kgHM was reached. The simulation recreated the classic striped CRUD pattern often seen on pulled fuel rods containing CRUD. It is concluded that the striping is caused by the flow swirl induced by spacer grid mixing vanes. Two anti-correlated effects contribute to the striping: (1) the flow swirl yields significant azimuthal temperature variations, which impact the locations where CRUD deposits, and (2) the flow swirl is correlated to increased shear stress along the cladding surface and subsequent erosion of the CRUD layer. The CIPS condition of the core is concluded to be primarily controlled by lithium tetraborate precipitation, referred to as boron hideout, which occurs in regions experiencing subcooled nucleate boiling as soluble boron and lithium species reach their solubility limit within the CRUD layer. Subsequently, a localized reduction in power occurs due to the high neutron absorption cross section of boron-10.

Turbulent Rayleigh–Bénard convection for a Prandtl number of 0.67

For the Rayleigh-number range 107 ≲ Ra ≲ 1011 we report measurements of the Nusselt number Nu and of properties of the large-scale circulation (LSC) for cylindrical samples of helium gas (Prandtl number Pr = 0.674) that have aspect ratio Γ ≡ D/L = 0.50 (D and L are the diameter and the height respectively) and are heated from below. The results for Nu are consistent with recent direct numerical simulations. We measured the amplitude δ of the azimuthal temperature variation induced by the LSC at the sidewall, and the LSC circulation-plane orientation θ0, at three vertical positions. For the entire Ra range the LSC involves a convection roll that is coherent over the height of the system. However, this structure frequently collapses completely at irregular time intervals and then reorganizes from the incoherent flow. At small δ the probability distribution p(δ) increases linearly from zero; for Γ = 1 and Pr = 4.38 this increase is exponential. No evidence of a two-roll structure, with one above the other, was observed. This differs from recent direct numerical simulations for Γ = 0.5 and Pr = 0.7, where a one-roll LSC was found to exist only for Ra ≲ 109 to 1010, and from measurements for Γ = 0.5 and Pr ≃ 5, where one- and two-roll structures were observed with transitions between them at random time intervals.

XMM‐NewtonandChandraObservations of the Galaxy Group NGC 5044. I. Evidence for Limited Multiphase Hot Gas

Using new XMM and Chandra observations, we present an analysis of the temperature structure of the hot gas within a radius of 100 kpc of the bright nearby galaxy group NGC 5044. A spectral deprojection analysis of data extracted from circular annuli reveals that a two-temperature model (2T) of the hot gas is favored over single-phase or cooling flow ( = 4.5 ? 0.2 M? yr-1) models within the central ~30 kpc. Alternatively, the data can be fitted equally well if the temperature within each spherical shell varies continuously from ~Th to Tc ~ Th/2, but no lower. The high spatial resolution of the Chandra data allows us to determine that the temperature excursion Th ? Tc required in each shell exceeds the temperature range between the boundaries of the same shell in the best-fitting single-phase model. This is strong evidence for a multiphase gas having a limited temperature range. We do not find any evidence that azimuthal temperature variations within each annulus on the sky can account for the range in temperatures within each shell. We provide a detailed investigation of the systematic errors on the derived spectral models considering the effects of calibration, plasma codes, bandwidth, variable NH, and background rate. We find that the RGS gratings and the EPIC and ACIS CCDs give fully consistent results when the same models are fitted over the same energy ranges for each instrument. The cooler component of the 2T model has a temperature (Tc ~ 0.7 keV) similar to the kinetic temperature of the stars. The hot phase has a temperature (Th ~ 1.4 keV) characteristic of the virial temperature of the ~1013 M? halo expected in the NGC 5044 group. However, in view of the morphological disturbances and X-ray holes visible in the Chandra image within R ? 10 kpc, bubbles of gas heated to ~Th in this region may be formed by intermittent AGN feedback. Some additional heating at larger radii may be associated with the evolution of the cold front near R ~ 50 kpc, as suggested by the sharp edge in the EPIC images.

Impact of Wind Velocity on the Performance of the RVACS Decay Heat Removal System

The impact of wind velocity on the performance of the reactor vessel auxiliary cooling system (RVACS) of an advanced liquid-metal reactor design is analyzed, and design modifications that mitigate adverse wind effects are investigated. In the reference design, the reactor is served by four communicating RVACS stacks, and each stack has two air inlets. In this two-inlet stack design, winds blowing in a direction 90 deg from the axis formed by the two stack inlets result in pressure distributions around the stacks that drastically change the desired airflow pattern in the RVACS. This leads to significantly elevated RVACS air temperatures and significant azimuthal guard vessel temperature variations. For example, a 27 m/s (60 mph) wind leads to an air temperature at the exit of the RVACS heated section that is ∼115°C higher than that under no-wind conditions. The addition of two more inlets per stack, one inlet per stack side, significantly improves RVACS performance. The air temperature at the exit of the heated RVACS section is significantly reduced below that of the two-inlet design, and this temperature decreases as the wind speed increases. An increase in wind speed from 3 to 27 m/s leads to an air temperature change from 186 to 165°C. The azimuthal temperature variation is also improved. At the top of the guard vessel, this variation is reduced from 62.5 to 8.5°Cat the low wind speed of 3 m/s and from 85.0 to 30.5°C at the high wind speed of 27 m/s.

Azimuthal temperature variation in the cladding of PWR fuel elements

This paper is concerned with the prediction of the azimuthal temperature distribution in the cladding of the fuel elements of a pressurized water reactor (PWR), in the event of a loss of coolant accident (LOCA). It is assumed that ballooning of the fuel elements results in local restrictions of flow area as the cladding of adjacent rods contact each other to form singly connected passages of the 4-cusp shape. A “fin-type” analysis of the circumferential heat transfer in the walls (the cladding) of the channel is made, account being taken of the variation of convective heat transfer coefficient around the channel perimeter. Predicted results for the maximum azimuthal temperature difference in the cladding are presented for a range of variables and percentage blockage, and a comparison is made with the predictions of other workers in this field of study.

An Analytical Study of Azimuthal Temperature Distribution in Pressurized Water Reactor Fuel Rods Under Steady-State and Blowdown Conditions

An analytical study was conducted to quantify the azimuthal temperature variations that might occur around a light water reactor nuclear fuel rod under steady-state and loss-of-coolant accident (LOCA) blowdown conditions. Significant azimuthal temperature variations on the cladding surface will result in a reduction of the total circumferential elongation of the cladding deforming in the alpha phase of the Zircaloy and thus a reduction in coolant channel blockage. Power skewing across the fuel rod and off-center location of the fuel pellet within the cladding are considered as possible causes for an asymmetrical temperature distribution on the cladding. The conclusions reached were as follows. Off-center location of a fuel pellet produces significant azimuthal temperature variations in the fuel during steady-state operation and small circumferential temperature differences in the cladding during the LOCA blowdown. Highly localized small strains-to-failure in Zircaloy out-of-pile ballooning experiments may not be representative of in-pile results where cladding temperatures are calculated to be relatively uniform.

Effect of fuel pellet eccentricity with cladding on fuel rod thermal behavior under reactivity initiated accident condition.

Fuel rod behavior under Reactivity Initiated Accident (RIA) conditions has been studied in the Nuclear Safety Research Reactor (NSRR), JAERI. In the experiments, cladding thermal behavior was observed to be influenced by the fuel pellet eccentricity to produce large azimuthal temperature variation in the cladding. The maximum azimuthal cladding temperature difference was measured to be as large as 150°C by thermocouples attached to opposite sides of the cladding around the circumference, though the thermocouples did not always detect the maximum temperature difference around the circumference. The actual temperature differences in the fuel rods subjected to less than 290 cal/g•UO2 were estimated to be 350°C at maximum based on metallographies. A simple calculation considering gap conductance variations also showed that the maximum temperature difference became 350°C under fully eccentrical condition in the fuel rod subjected to 260 cal/g•UO2. Moreover, as the rod damage such as cladding deformation, melting and failure occurs unevenly around the circumference due to the fuel pellet eccentricity in general, the fuel pellet eccentricity should influence the fuel rod failure under RIA conditions.

Steady-state, local temperature fields with turbulent liquid sodium flow in nominal and disturbed bundle geometries with spacer grids

Abstract The operating reliability of nuclear reactors calls for a reliable strength analysis of the highly loaded core elements, one of its prerequisites being the reliable determination of the three-dimensional velocity and temperature fields. To verify thermohydraulics computer programs, extensive local temperature measurements in the rod claddings of the critical bundle zone were performed on a heated 19-rod bundle model with sodium flow and provided with spacer grids (P/D = 1.30; W/D = 1.19). The essential results are: - Outside the spacer grids, the azimuthal temperature variations of the side and corner rods are approximately 10-fold those of rods in the central bundle zone. - The spacer grids investigated give rise to great local temperature peaks and correspondingly great temperature gradients in the axial and azimuthal directions immediately around the support points. - Continuous reduction of a subchannel by rod bowing results in substantial rises of temperature which, however, are limited to adjacent cladding tubes.