Search NASA⌕ Search

SEARCH · Search NASA

Results for “temperature gradient”

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 343 records · Page 19

Quasilinear Carbon Transport In An Impurity Hole Plasma In LHD

Comprehensive electrostatic gyrokinetic linear stability calculations for ion-scale microinstabilities in an LHD plasma with an ion-ITB and carbon "impurity hole" are used to make quasilinear estimates of particle flux to explore whether microturbulence can explain the observed outward carbon fluxes that flow "up" the impurity density gradient. The ion temperature is not stationary in the ion-ITB phase of the simulated discharge, during which the core carbon density decreases continuously. To fully sample these varying conditions the calculations are carried out at three radial locations and four times. The plasma parameter inputs are based on experimentally measured profiles of electron and ion temperature, as well as electron and carbon density. The spectroscopic line-average ratio of hydrogen and helium densities is used to set the density of these species. Three ion species (H,He,C) and the electrons are treated kinetically, including collisions. Electron instability drive does enhance the growth rate significantly, but the most unstable modes have characteristics of ion temperature gradient (ITG) modes in all cases. As the carbon density gradient is scanned between the measured value and zero, the quasilinear carbon flux is invariably inward when the carbon density profile is hollow, so turbulent transport due to the instabilities considered here does not explain the observed outward flux of impurities in impurity hole plasmas. The stiffness of the quasilinear ion heat flux is found to be 1.7-2.3, which is lower than several estimates in tokamaks.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Thermal transpiration: A molecular dynamics study

Thermal transpiration is a phenomenon where fluid molecules move from the cold end towards the hot end of a channel under the influence of longitudinal temperature gradient alone. Although the phenomenon of thermal transpiration is observed at rarefied gas conditions in macro systems, the phenomenon can occur at atmospheric pressure if the characteristic dimensions of the channel is less than 100 nm. The flow through these nanosized channels is characterized by the free molecular flow regimes and continuum theory is inadequate to describe the flow. Thus a non-continuum method like molecular dynamics (MD) is necessary to study such phenomenon. In the present work, MD simulations were carried out to investigate the occurance of thermal transpiration in copper and platinum nanochannels at atmospheric pressure conditions. The mean pressure of argon gas confined inside the nano channels was maintained around 1 bar. The channel height is maintained at 2nm. The argon atoms interact with each other and with the wall atoms through the Lennard-Jones potential. The wall atoms are modelled using an EAM potential. Further, separate simulations were carried out where a Harmonic potential is used for the atom-atom interaction in the platinum channel. A thermally insulating wall was introduced between the low and high temperature regions and those wall atoms interact with fluid atoms through a repulsive potential. A reduced cut off radius were used to achieve this. Thermal creep is induced by applying a temperature gradient along the channel wall. It was found that flow developed in the direction of the increasing temperature gradient of the wall. An increase in the volumetric flux was observed as the length of the cold and the hot regions of the wall were increased. The effect of temperature gradient and the wall-fluid interaction strength on the flow parameters have been studied to understand the phenomenon better.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Weibel-like instability in magnetohydrodynamics

In magnetohydrodynamics (MHD), a density perturbation perpendicular to an electron temperature gradient generates a magnetic field around itself that acts to increase the perturbation, which can lead to instability. An MHD dispersion relation is obtained for perturbations perpendicular to a fixed electron temperature gradient with an initial in-plane magnetic field, including resistivity, viscosity, and the electrothermal coefficient. Instability occurs for sufficiently small electron temperature-gradient scale lengths determined by the ion collisionless skin depth. Both viscosity and resistivity are required to prevent growth at arbitrarily small spatial scales and to give a physical result for the fastest growing mode. The perpendicular electrothermal coefficient is only significant for a narrow range of low electron Hall parameters, causing a modest reduction in magnetic field growth and modifying the criteria for instability in the presence of viscosity. If the definition of the Weibel instabilityis extended to include all instabilities due to anisotropy in the electron velocity distribution, then this is a Weibel-like instability because an electron temperature gradient implies an anisotropic electron velocity distribution. Here, the implications for the formation of filaments in laser-produced plasmas and for the verification of MHD codes are considered.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Onset of surface-tension-driven Benard convection

An experimental investigation of the onset of convection in shallow fluid layers heated uniformly from below and cooled from above by an air layer has been made. If the depth of the silicone layer is smaller than 2 mm the onset of convection takes place in two stages. There is first a weak pattern, which is characterized by its appearance at ever smaller temperature gradients as the depth of the fluid is decreased. When the temperature difference across the fluid is increased a second strong pattern forms near the predicted critical Marangoni number. The cells in this pattern are hexagonal and seem to be what one has always referred to as Benard cells. The temperature gradient at which this pattern appears increases with decreased depth. The heat transfer through the fluid has been measured. The critical temperature gradient for the formation of the hexagonal pattern has been determined from the break of the heat transfer curve.

Koschmieder, E. L.↗

Developing Soret Transport Coefficients for Lanthanides from Experimental Results

Understanding the nature of lanthanide transport in metallic fuels and its relation to fuel clad chemical interaction (FCCI) has been has been a topic that is long on interest and short on data. Recent experiments have shown a strong temperature-gradient based diffusion (Soret diffusion) and have thus provided vitally important information for modelers. However, details about the actual coefficients and many of the relevant processes that crucially participate in the cerium migration and precipitation process have not yet been disclosed. Using a combination of computing tools (primarily python and BISON) we are exploring methods that may extract useful information from visual data. In this specific example, we are using what is formally called 'blob detection' algorithms to analyzed the precipitate distribution of cerium in a U-Zr metal that has had a fixed temperature gradient enforced upon it for a specified period of time. These algorithms, once properly adapted to the problem at hand, may be useful at extracting information from radial and axial scans of metallic fuels which are much more abundant than direct measurements of lanthanides. We have developed a model that qualitatively works to move mass along a temperature gradient and precipitate heavily in the cold regions. Additionally, we have observed that the solubility of cerium in U-Zr held under this kind of temperature gradient is subject to some amount of U-Zr constituent redistribution. Our results predict movement of the Zr toward the hot end of the rod followed by a U-rich layer with the cooler 3/5 of the rod staying close the the initial U-Zr distribution. A more complete analysis will define groups of correlated parameters that correctly simulate the resulting distribution of precipitates. This restricted parameters space will give smaller bounds to modelers and will expedite the analysis of future works as new and more accurate data are made available. A follow-up experiment that started with 1/10 of the initial concentration of cerium may provide valuable additional information with regard to the effect of the rates of precipitation and dissolution as well as demonstrating the potential rapidity of migration in a lower burnup scenario. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Merging high resolution geophysical and geochemical surveys to reduce exploration risk at glass buttes, Oregon

The engineering and studies phase of the Glass Buttes project was aimed at reducing risk during the early stages of geothermal project development. The project’s inclusion of high resolution geophysical and geochemical surveys allowed Ormat to evaluate the value of these surveys both independently and in combination to quantify the most valuable course of action for exploration in an area where structure, permeability, and temperature are the most pressing questions. The sizes of the thermal anomalies at Glass Buttes are unusually large. Over the course of Phase I Ormat acquired high resolution LIDAR data to accurately map fault manifestations at the surface and collected detailed gravity and aeromagnetic surveys to map subsurface structural features. In addition, Ormat collected airborne hyperspectral data to assist with mapping the rock petrology and mineral alteration assemblages along Glass Buttes faults and magnetotelluric (MT) survey to try to better constrain the structures at depth. Direct and indirect identification of alteration assemblages reveal not only the geochemical character and temperature of the causative hydrothermal fluids but can also constrain areas of upflow along specific fault segments. All five datasets were merged along with subsurface lithologies and temperatures to predict the most likely locations for high permeability and hot fluids. The Glass Buttes temperature anomalies include 2 areas, totaling 60 km2 (23 mi2) of measured temperature gradients over 165° C/km (10° F/100ft). The Midnight Point temperature anomaly includes the Strat-1 well with 90°C (194 °F) at 603 m (1981 ft) with a 164 °C/km (10°F/100ft) temperature gradient at bottom hole and the GB-18 well with 71°C (160 °F) at 396 m (1300 ft) with a 182°C/km (11°F/100ft) gradient. The primary area of alteration and elevated temperature occurs near major fault intersections associated with Brothers Fault Zone and Basin and Range systems. Evidence for faulting is observed in each data set as follows. Field observations include fault plane orientations, complicated fault intersections, and hydrothermal alteration apparently pre-dating basalt flows. Geophysical anomalies include large, linear gradients in gravity and aeromagnetic data with magnetic lows possibly associated with alteration. Resistivity low anomalies also appear to have offsets associated with faulting. Hyperspectral and XRF identified alteration and individual volcanic flow units, respectively. When incorporated into a 3D geologic model, the fault intersections near the highest proven temperature and geophysical anomalies provide the first priority targets at Midnight Point. Ormat geologists selected the Midnight Point 52-33 drilling target based on a combination of pre-existing drilling data, geologic field work, geophysical interpretation, and geochemical analysis. Deep temperatures of well 52-33 was lower than anticipated. Temperature gradients in the well mirrored those found in historical drilling, but they decreased below 1500 ft and were isothermal below 2000 ft.

15 GEOTHERMAL ENERGY↗

CONDUCTION IN LOW MACH NUMBER FLOWS. I. LINEAR AND WEAKLY NONLINEAR REGIMES

Thermal conduction is an important energy transfer and damping mechanism in astrophysical flows. Fourier's law, in which the heat flux is proportional to the negative temperature gradient, leading to temperature diffusion, is a well-known empirical model of thermal conduction. However, entropy diffusion has emerged as an alternative thermal conduction model, despite not ensuring the monotonicity of entropy. This paper investigates the differences between temperature and entropy diffusion for both linear internal gravity waves and weakly nonlinear convection. In addition to simulating the two thermal conduction models with the fully compressible Navier-Stokes equations, we also study their effects in the reduced ''soundproof'' anelastic and pseudoincompressible (PI) equations. We find that in the linear and weakly nonlinear regime, temperature and entropy diffusion give quantitatively similar results, although there are some larger errors in the PI equations with temperature diffusion due to inaccuracies in the equation of state. Extrapolating our weakly nonlinear results, we speculate that differences between temperature and entropy diffusion might become more important for strongly turbulent convection.

79 ASTRONOMY AND ASTROPHYSICS↗

What lies beneath: Vertical temperature heterogeneity in a Mediterranean woodland savanna

As the availability of satellite and airborne thermal infrared remote sensing (TIR-RS) data increases and their spatial, temporal, and spectral resolutions improve, researchers are finding diverse applications for TIR-RS measurements. TIR-RS is now commonly applied in regional- and continental-scale analyses, such as those focused on fire and surface energy balance. However, its application lags in plant physiology and ecology, for which a finer-scale understanding of plant canopy temperatures would be useful to elucidate plant water dynamics, for example. In particular, while methods to disaggregate TIR-RS pixels in horizontal space have advanced, possible vertical stratification of plant canopy temperature and its implications for understanding the correspondence between TIR-RS and finer-scale, field-based thermal measurements (e.g. made with a thermal camera) remain unexplored. Here, we use data from a thermal camera deployed concurrently with the recent ECOSTRESS mission to quantify vertical temperature gradients within tree canopies and temperatures of over- vs. under-story plants in a Mediterranean woodland savanna. We then leverage diverse ancillary data to maximize the geometric comparability of ECOSTRESS and thermal camera measurements, in order to assess the extent to which the two forms of thermal measurements correspond. Specifically, we ask: (1) What are the patterns of intra-canopy and over- vs. under-story vertical temperature in a Mediterranean woodland savanna?, and (2) How can vertically-resolved, but spatially-limited field-based temperature measurements be reconciled with spatially-extensive, but surface-only, temperature measurements of a space-borne remote sensor? Further, we found consistent patterns of vertical thermal heterogeneity both within tree canopies and between ecosystem over- and under-stories. The daytime difference between the top and bottom thirds of blue oak canopies was, on average, 0.48 ° C – and sometimes several times larger. Notably, canopy tops are cooler, likely associated with the under-story grass reaching daytime temperatures often exceeding over-story temperatures by 10° C. Given the consistency of the intra-canopy temperature gradients, we expected the ECOSTRESS sensor would be in better agreement with camera measurements of canopy tops than bulk canopies or canopy bottoms. However, within-canopy gradients were overwhelmed by other sources of disagreement between the measurements, in part associated with upscaling camera measurements across space. Overall, thermal camera and ECOSTRESS measurements were largely in agreement at night (pixel RMSE = 1.1°C), but they were more divergent during times of low (but >0 W/m2) and high incoming solar radiation (daytime pixel RMSE = 3.5°C).

54 ENVIRONMENTAL SCIENCES↗

How does negative triangularity mitigate ITG turbulence and transport?

Improved confinement in negative triangularity (NT) experiments is attributed to reduced fluxes driven by micro-turbulence. The physical mechanism of why thermal confinement improves in NT relative to PT is unknown. This study employs gyrokinetic flux tube simulations using the GENE code with local Miller equilibrium to elucidate the physical mechanisms behind the beneficial effects of NT flux surface shapes. The focus is on collisionless ion temperature gradient (ITG) driven turbulence with adiabatic electrons. The kinetic profiles are held fixed across a scan of triangularity values, thus enabling comparisons on a level playing field. The reduced linear growth rates for NT is shown to be due to a reduced eigenmode averaged magnetic drift frequency and a wider, stronger negative local magnetic shear region about the outboard mid-plane. The nonlinear heat flux is lower for NT than that for PT, due to reduced radial correlation length and increased correlation time (τ c ) of fluctuations. These, in turn, are due to a comparatively higher level of self-generated zero-frequency E × B zonal shearing rate ω E in NT as compared to PT. Though the linear zonal potential residual is lower for NT, the nonlinearly generated E × B zonal shearing rate is higher for NT than for PT. This outcome is linked to the distinctive features of the radial wavenumber spectra of the zonal potential and the zonal shearing rate. The dimensionless parameter ω E τ c is suggested as a figure of merit. This is higher for NT than for PT. Thus, the reduced heat diffusivity for NT is linked to increased ω E τ c . Self-generated temperature corrugations (i.e. zonal temperature gradients) are much weaker than the background mean temperature gradient. Nevertheless, temperature corrugations are more pronounced in NT than in PT.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Developmental testing of a programmable multizone furnace

A multizone furnace was evaluated for its potential utilization for process experimentation on board the Space Shuttle. A temperature gradient can be created through the use of a series of connected temperature zones and can be translated by the coordinated sequencing of zone temperatures. The Bridgman-Stockbarger thermal configuration for directional solidification was implemented so that neither the sample nor furnace was translated. The thermal behavior of the furnace was measured and characterized. Limitations due to both thermal and electronic (computer) factors are identified. The results indicate that the multizone design is limited to low temperature gradients because of the indirect furnace-to-sample thermal coupling needed to blend the discrete thermal zones. The multizone furnace design inherently consumes more power than a similar (two temperature) conventional Bridgman type directional solidification furnace because every zone must be capable of the high cooling rates needed to produce the maximum desired temperature drop. Typical achievable static temperature gradients for the furnace tested were between 6 and 75 C/in. The maximum gradient velocity was approximately 10 in./hr. Several aspects of the tested system could be improved, but the dependence of the multizone design on high heat loss will limit Space Shuttle applications in the form tested unless additional power is available. The multizone furnace offers great flexibility but requires a high level of operator understanding for full advantage to be obtained.

Ting, E. Y.↗

The Climate Response to Stratospheric Aerosol Geoengineering Can Be Tailored Using Multiple Injection Locations

By injecting different amounts of SO 2 at multiple different latitudes, the spatial pattern of aerosol optical depth (AOD) can be partially controlled. This leads to the ability to influence the climate response to geoengineering with stratospheric aerosols, providing the potential for design. In this work, we use simulations from the fully coupled whole-atmosphere chemistry climate model CESM1(WACCM) to demonstrate that by appropriately combining injection at just four different locations, 30°S, 15°S, 15°N, and 30°N, then three spatial degrees of freedom of AOD can be achieved: an approximately spatially uniform AOD distribution, the relative difference in AOD between Northern and Southern Hemispheres, and the relative AOD in high versus low latitudes. For forcing levels that yield 1–2°C cooling, the AOD and surface temperature response are sufficiently linear in this model so that the response to different combinations of injection at different latitudes can be estimated from single-latitude injection simulations; nonlinearities associated with both aerosol growth and changes to stratospheric circulation will be increasingly important at higher forcing levels. Optimized injection at multiple locations is predicted to improve compensation of CO 2 -forced climate change relative to a case using only equatorial aerosol injection (which overcools the tropics relative to high latitudes). The additional degrees of freedom can be used, for example, to balance the interhemispheric temperature gradient and the equator to pole temperature gradient in addition to the global mean temperature. Further research is needed to better quantify the impacts of these strategies on changes to long-term temperature, precipitation, and other climate parameters.

54 ENVIRONMENTAL SCIENCES↗

Proposed Drill Sites

Proposed drill sites for intermediate depth temperature gradient holes and/or deep resource confirmation wells. Temperature gradient contours based on shallow TG program and faults interpreted from seismic reflection survey are shown, as are two faults interpreted by seismic contractor Optim but not by Oski Energy, LLC.

15 GEOTHERMAL ENERGY↗

Phenotypic plasticity of fungal traits in response to moisture and temperature

Phenotypic plasticity of traits is commonly measured in plants to improve understanding of organismal and ecosystem responses to climate change but is far less studied for microbes. Specifically, decomposer fungi are thought to display high levels of phenotypic plasticity and their functions have important implications for ecosystem dynamics. Assessing the phenotypic plasticity of fungal traits may therefore be important for predicting fungal community response to climate change. Here, we assess the phenotypic plasticity of 15 fungal isolates (12 species) from a Southern California grassland. Fungi were incubated on litter at five moisture levels (ranging from 4–50% water holding capacity) and at five temperatures (ranging from 4–36 °C). After incubation, fungal biomass and activities of four extracellular enzymes (cellobiohydrolase (CBH), β -glucosidase (BG), β -xylosidase (BX), and N-acetyl- β - D -glucosaminidase (NAG)) were measured. We used response surface methodology to determine how fungal phenotypic plasticity differs across the moisture-temperature gradient. We hypothesized that fungal biomass and extracellular enzyme activities would vary with moisture and temperature and that the shape of the response surface would vary between fungal isolates. We further hypothesized that more closely related fungi would show more similar response surfaces across the moisture-temperature gradient. In support of our hypotheses, we found that plasticity differed between fungi along the temperature gradient for fungal biomass and for all the extracellular enzyme activities. Plasticity also differed between fungi along the moisture gradient for BG activity. These differences appear to be caused by variation mainly at the moisture and temperature extremes. We also found that more closely related fungi had more similar extracellular enzymes activities at the highest temperature. Altogether, this evidence suggests that with global warming, fungal biodiversity may become increasingly important as functional traits tend to diverge along phylogenetic lines at higher temperatures.

59 BASIC BIOLOGICAL SCIENCES↗

Thermal transport dynamics in the quasi-single helicity state

A dynamical model describing oscillations between multiple and single helicity configurations in the quasi-single helicity (QSH) state of the reversed field pinch [P. W. Terry and G. G. Whelan, Plasma Phys. Controlled Fusion 56, 094003 (2014)] is extended to include electron temperature profile dynamics. It is shown that QSH dynamics is linked to the electron temperature profile because the suppression of mode coupling between tearing modes proposed to underlie QSH also suppresses magnetic-fluctuation-induced thermal transport. Above the threshold of dominant-mode shear that marks the transition to QSH, the model produces temperature-gradient steepening in the strong shear region. Oscillations of the dominant and secondary mode amplitudes give rise to oscillations of the temperature gradient. The phasing and amplitude of temperature gradient oscillations relative to those of the dominant mode are in agreement with experiment. This provides further evidence that the model, while heuristic, captures key physical aspects of the QSH state.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Numerical Investigation of the characteristics of pressurized biomass-oxy-combustion

As a second generation of oxy-fuel combustion technologies, pressurized oxy-combustion has a potential to improve the process efficiency and economy. With the modern decarbonization needs, biomass combustion is considered to be nearly CO2 neutral, and thus co-firing biomass with coal in an oxy-combustion process can achieve a negative CO2 balance. The size of the biomass fuel particles is a key parameter in an entrained flow combustor because large biomass particles, exceeding some threshold size, may only partially be burnt due to the restrictions on the heat transport, thereby diminishing the combustion efficiency. While such a threshold particle size has been scrutinized at an atmospheric pressure, the present study extends the analysis to elevated pressures. Specifically, two scenarios were set up, for experimental-scale and full-scale pressurized oxy-combustors, with the threshold sizes of biomass particles appearing 0.5 mm and 5 mm, respectively, without temperature gradient in biomass particles. With temperature gradient, the parametric study included the particles in the range from 0.15 mm to 5 mm. It is shown that the thermal-thin model is applicable to spherical biomass particles not exceeding 0.15 mm. The impact of the heating rate on biomass particle devolatilization was also investigated and the biomass devolatilization could happen in less than 1 second; the same amount of heat energy could induce different volatile release rates.

Li, Lei↗

Thermal-Mechanical Stability of Single Crystal Oxide Refractive Concentrators for High-Temperature Solar Thermal Propulsion

Single crystal oxides such as yttria-stabilized zirconia (Y2O3-ZrO2), yttrium aluminum garnet (Y3Al5O12, or YAG), magnesium oxide (MgO) and sapphire (Al2O3) are candidate refractive secondary concentrator materials for high temperature solar propulsion applications. However, thermo-mechanical reliability of these components in severe thermal environments during the space mission sun/shade transition is of great concern. Simulated mission tests are important for evaluating these candidate oxide materials under a variety of transient and steady-state heat flux conditions, and thus provide vital information for the component design. In this paper, a controlled heat flux thermal shock test approach is established for the single crystal oxide materials using a 3.0 kW continuous wave CO2 laser, with a wavelength 10.6 micron. Thermal fracture behavior and failure mechanisms of these oxide materials are investigated and critical temperature gradients are determined under various temperature and heating conditions. The test results show that single crystal sapphire is able to sustain the highest temperature gradient and heating-cooling rate, and thus exhibit the best thermal shock resistance, as compared to the yttria-stabilized zirconia, yttrium aluminum garnet, and magnesium oxide.

Zhu, Dongming↗

Thermal-Mechanical Stability of Single Crystal Oxide Refractive Concentrators for High-Temperature Solar Thermal Propulsion

Single crystal oxides such as yttria-stabilized zirconia (Y2O3-ZrO2), yttrium aluminum garnet (Y3Al5O12, or YAG), magnesium oxide (MgO) and sapphire (Al2O3) are candidate refractive secondary concentrator materials for high temperature solar propulsion applications. However, thermo-mechanical reliability of these components in severe thermal environments during the space mission sun/shade transition is of great concern. Simulated mission tests are important for evaluating these candidate oxide materials under a variety of transient and steady-state heat flux conditions, and thus provide vital information for the component design. In this paper, a controlled heat flux thermal shock test approach is established for the single crystal oxide materials using a 3.0 kW continuous wave CO2 laser, with a wavelength 10.6 micron. Thermal fracture behavior and failure mechanisms of these oxide materials are investigated and critical temperature gradients are determined under various temperature and heating conditions. The test results show that single crystal sapphire is able to sustain the highest temperature gradient and heating-cooling rate, and thus exhibit the best thermal shock resistance, as compared to the yttria-stabilized zirconia, yttrium aluminum garnet and magnesium oxide.

Zhu, Dongming↗

Modeling fuel assembly distortion in light water reactors using BISON

Individual fuel rods in a light water reactor fuel assembly (FA) often have different power profiles that result in a spatial temperature gradient in the FA. This temperature gradient causes some rods to elongate more than others which could in turn distort the entire FA. Such FA distortion has been known to cause operational difficulties with refueling and with the insertion of control rods [1]. Numerical modeling of FA distortion due to differential elongation of fuel rods requires: (i) modeling the thermal elongation as well as the creep and irradiation induced elongation of the fuel rods, and (ii) modeling the bending response of the FA due to the differential elongation of the fuel rods. Past studies on this topic have been focused on the mechanical deformation of the FA with simplistic modeling of fuel rod elongation. In this paper, a simulation approach is presented that directly couples 2D axisymmetric fuel performance models of the individual fuel rods in a FA with structural beam models of the fuel rods in a FA using the Bison code [2]. This allows all of the mechanisms internal to the fuel rod leading to cladding elongation to be accounted for in the FA distortion model. To demonstrate this capability, the bowing of a simplified and idealized 6x6 FA (Fig. 1a), modeled using beam elements, is analyzed under different spatially varying power profiles. For this study, the guide tubes are rigidly connected to the spacer grid. Also, the fuel rods are rigidly connected to the spacer grid in all directions except the vertical direction (Y) along which slipping can occur between the fuel rods and spacer grids, depending on the friction coefficient between them. Coupling the thermo-mechanical response of individual fuel rods to the global mechanical behavior of the FA can shed light on the mechanisms that cause distortion of the FA and potential solutions to address this issue.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗