Example of Heat Equation – Problem with Solution

Example of Heat Equation – Problem with Solution

Heat Conduction in a Large Plane Wall

Example of Heat Equation – Problem with Solution

Consider the plane wall of thickness 2L, in which there is uniform and constant heat generation per unit volume, qV [W/m3]. The centre plane is taken as the origin for x and the slab extends to + L on the right and – L on the left. For constant thermal conductivity k, the appropriate form of the heat equation, is:

heat equation - heat generation - equation

The general solution of this equation is:

heat conduction equation - general solution

where C1 and C2 are the constants of integration.

1)

Heat Conduction in a Large Plane WallCalculate the temperature distribution, T(x), through this thick plane wall, if:

  • the temperatures at both surfaces are 15.0°C
  • the thickness this wall is 2L = 10 mm.
  • the materials conductivity is k = 2.8 W/m.K (corresponds to uranium dioxide at 1000°C)
  • the volumetric heat rate is qV = 106 W/m3

In this case, the surfaces are maintained at given temperatures Ts,1 and Ts,2. This corresponds to the Dirichlet boundary condition. Moreover, this problem is thermally symmetric and therefore we may use also thermal symmetry boundary condition. The constants may be evaluated using substitution into the general solution and are of the form:

heat conduction equation - boundary conditions

The resulting temperature distribution and the centerline (x = 0) temperature (maximum) in this plane wall at these specific boundary conditions will be:

heat conduction equation - solution

The heat flux at any point, qx [W.m-2], in the wall may, of course, be determined by using the temperature distribution and with the Fourier’s law. Note that, with heat generation the heat flux is no longer independent of x, the therefore:

heat condution through wall

Example of Heat Equation – Problem with Solution

Heat Conduction in a Fuel Rod

Most of PWRs use the uranium fuel, which is in the form of uranium dioxide. Uranium dioxide is a black semiconducting solid with very low thermal conductivity. On the other hand the uranium dioxide has very high melting point and has well known behavior. The UO2 is pressed into cylindrical pellets, these pellets are then sintered into the solid.

These cylindrical pellets are then loaded and encapsulated within a fuel rod (or fuel pin), which is made of zirconium alloys due to its very low absorption cross-section (unlike the stainless steel). The surface of the tube, which covers the pellets, is called fuel cladding.

See also: Thermal Conduction of Uranium Dioxide

Thermal and mechanical behavior of fuel pellets and fuel rods constitute one of three key core design disciplines. Nuclear fuel is operated under very inhospitable conditions (thermal, radiation, mechanical) and must withstand more than normal conditions operation. For example temperatures in the centre of fuel pellets reach more than 1000°C (1832°F) accompanied by fission-gas releases. Therefore detailed knowledge of temperature distribution within a single fuel rod is essential for safe operation of nuclear fuel. In this section we will study heat conduction equation in cylindrical coordinates using Dirichlet boundary condition with given surface temperature (i.e. using Dirichlet boundary condition). Comprehensive analysis of fuel rod temperature profile will be studied in separate section.

Temperature in the centerline of a fuel pellet

Consider the fuel pellet of radius rU = 0.40 cm, in which there is uniform and constant heat generation per unit volume, qV [W/m3]. Instead of volumetric heat rate qV [W/m3], engineers often use the linear heat rate, qL [W/m], which represents the heat rate of one meter of fuel rod. The linear heat rate can be calculated from the volumetric heat rate by:

linear heat rate vs volumetric heat rate

The centreline is taken as the origin for r-coordinate. Due to symmetry in z-direction and in azimuthal direction, we can separate of variables and simplify this problem to one-dimensional problem. Thus, we will solve for the temperature as function of radius, T(r), only. For constant thermal conductivity, k, the appropriate form of the cylindrical heat equation, is:

heat equation - cylindrical - 2

The general solution of this equation is:

heat equation - cylindrical - general solution

where C1 and C2 are the constants of integration.

Thermal conduction - fuel pelletCalculate the temperature distribution, T(r), in this fuel pellet, if:

  • the temperatures at the surface of the fuel pellet is TU = 420°C
  • the fuel pellet radius rU = 4 mm.
  • the averaged material’s conductivity is k = 2.8 W/m.K (corresponds to uranium dioxide at 1000°C)
  • the linear heat rate is qL = 300 W/cm and thus the volumetric heat rate is qV = 597 x 106 W/m3

In this case, the surface is maintained at given temperatures TU. This corresponds to the Dirichlet boundary condition. Moreover, this problem is thermally symmetric and therefore we may use also thermal symmetry boundary condition. The constants may be evaluated using substitution into the general solution and are of the form:

heat equation - cylindrical - boundary conditions

The resulting temperature distribution and the centerline (r = 0) temperature (maximum) in this cylindrical fuel pellet at these specific boundary conditions will be:

heat equation - cylindrical - solution

The radial heat flux at any radius, qr [W.m-1], in the cylinder may, of course, be determined by using the temperature distribution and with the Fourier’s law. Note that, with heat generation the heat flux is no longer independent of r.

The following figure shows the temperature distribution in the fuel pellet at various power levels.

Temperature distribution - nuclear fuel

______

The temperature in an operating reactor varies from point to point within the system. As a consequence, there is always one fuel rod and one local volume, that are hotter than all the rest. In order to limit these hot places the peak power limits must be introduced. The peak power limits are associated with a boiling crisis and with the conditions which could cause fuel pellet melt. However, metallurgical considerations place an upper limits on the temperature of the fuel cladding and the fuel pellet. Above these temperatures there is a danger that the fuel may be damaged. One of the major objectives in the design of a nuclear reactors is to provide for the removal of the heat produced at the desired power level, while assuring that the maximum fuel temperature and the maximum cladding temperature are always below these predetermined values.

Temperature distribution in Fuel Cladding

Cladding is the outer layer of the fuel rods, standing between the reactor coolant and the nuclear fuel (i.e. fuel pellets). It is made of a corrosion-resistant material with low absorption cross section for thermal neutrons, usually zirconium alloyCladding prevents radioactive fission products from escaping the fuel matrix into the reactor coolant and contaminating it. Cladding constitute one of barriers in ‘defence-in-depth‘ approach.

Consider the fuel cladding of inner radius rZr,2 = 0.408 cm and outer radius rZr,1 = 0.465 cm. In comparison to fuel pellet, there is almost no heat generation in the fuel cladding (cladding is slightly heated by radiation). All heat generated in the fuel must be transferred via conduction through the cladding and therefore the inner surface is hotter than the outer surface.

To find the temperature distribution through the cladding we must solve the heat conduction equation. Due to symmetry in z-direction and in azimuthal direction, we can separate of variables and simplify this problem to one-dimensional problem. Thus, we will solve for the temperature as function of radius, T(r), only. In this example, we will assume that there is strictly no heat generation within the cladding. For constant thermal conductivity, k, the appropriate form of the cylindrical heat equation, is:

heat equation - cladding

The general solution of this equation is:

heat equation - cladding - general solution

where C1 and C2 are the constants of integration.

1)

Thermal conduction - fuel pelletCalculate the temperature distribution, T(r), in this fuel cladding, if:

  • the temperature at the inner surface of the cladding is TZr,2 = 360°C
  • the temperature of reactor coolant at this axial coordinate is Tbulk = 300°C
  • the heat transfer coefficient (convection; turbulent flow) is h = 41 kW/m2.K.
  • the averaged material’s conductivity is k = 18 W/m.K
  • the linear heat rate of the fuel is qL = 300 W/cm and thus the volumetric heat rate is qV = 597 x 106 W/m3

From the basic relationship for heat transfer by convection we can calculate the outer surface of the cladding as:

newtons law - cladding

As can be seen, also in this case we have given surface temperatures TZr,1 and TZr,2. This corresponds to the Dirichlet boundary condition. The constants may be evaluated using substitution into the general solution and are of the form:

heat equation - general solution - cladding

Solving for C1 and C2 and substituting into the general solution, we then obtain:

heat equation - cladding - solution

∆T – cladding surface – coolant

Detailed knowledge of geometry, outer radius of cladding, linear heat rate, convective heat transfer coefficient and the coolant temperature determines ∆T between the coolant (Tbulk) and the cladding surface (TZr,1). Therefore we can calculate the cladding surface temperature (TZr,1) simply using the Newton’s Law:

∆T - cladding surface - coolant - newtons law

∆T in fuel cladding

Detailed knowledge of geometry, outer and inner radius of cladding, linear heat rate, and the cladding surface temperature (TZr,1) determines ∆T between outer and inner surfaces of cladding. Therefore we can calculate the inner cladding surface temperature (TZr,2) simply using the Fourier’s Law:

∆T in fuel cladding - fouriers law

 
References:
Heat Transfer:
  1. Fundamentals of Heat and Mass Transfer, 7th Edition. Theodore L. Bergman, Adrienne S. Lavine, Frank P. Incropera. John Wiley & Sons, Incorporated, 2011. ISBN: 9781118137253.
  2. Heat and Mass Transfer. Yunus A. Cengel. McGraw-Hill Education, 2011. ISBN: 9780071077866.
  3. Fundamentals of Heat and Mass Transfer. C. P. Kothandaraman. New Age International, 2006, ISBN: 9788122417722.
  4. U.S. Department of Energy, Thermodynamics, Heat Transfer and Fluid Flow. DOE Fundamentals Handbook, Volume 2 of 3. May 2016.

Nuclear and Reactor Physics:

  1. J. R. Lamarsh, Introduction to Nuclear Reactor Theory, 2nd ed., Addison-Wesley, Reading, MA (1983).
  2. J. R. Lamarsh, A. J. Baratta, Introduction to Nuclear Engineering, 3d ed., Prentice-Hall, 2001, ISBN: 0-201-82498-1.
  3. W. M. Stacey, Nuclear Reactor Physics, John Wiley & Sons, 2001, ISBN: 0- 471-39127-1.
  4. Glasstone, Sesonske. Nuclear Reactor Engineering: Reactor Systems Engineering, Springer; 4th edition, 1994, ISBN: 978-0412985317
  5. W.S.C. Williams. Nuclear and Particle Physics. Clarendon Press; 1 edition, 1991, ISBN: 978-0198520467
  6. G.R.Keepin. Physics of Nuclear Kinetics. Addison-Wesley Pub. Co; 1st edition, 1965
  7. Robert Reed Burn, Introduction to Nuclear Reactor Operation, 1988.
  8. U.S. Department of Energy, Nuclear Physics and Reactor Theory. DOE Fundamentals Handbook, Volume 1 and 2. January 1993.
  9. Paul Reuss, Neutron Physics. EDP Sciences, 2008. ISBN: 978-2759800414.

Advanced Reactor Physics:

  1. K. O. Ott, W. A. Bezella, Introductory Nuclear Reactor Statics, American Nuclear Society, Revised edition (1989), 1989, ISBN: 0-894-48033-2.
  2. K. O. Ott, R. J. Neuhold, Introductory Nuclear Reactor Dynamics, American Nuclear Society, 1985, ISBN: 0-894-48029-4.
  3. D. L. Hetrick, Dynamics of Nuclear Reactors, American Nuclear Society, 1993, ISBN: 0-894-48453-2. 
  4. E. E. Lewis, W. F. Miller, Computational Methods of Neutron Transport, American Nuclear Society, 1993, ISBN: 0-894-48452-4.

See above:

Heat Equation