A high-temperature, gas-cooled nuclear reactor consists of a composite, cylindrical wall for which a thorium fuel element (kth = 57 W/m-K) is encased in graphite (kgr = 3 W/m-K) and a very high flowrate of gaseous helium over the outer wall of the graphite. Consider conditions for which the helium flow where T[infinity] = 600 K and the heat transfer convection coefficient h = 2,000 W/m2K. The radius values are as follows: r1 = 8 mm, r2 = 11 mm, and r3 = 14 mm. Assume the insulation material in the center of the vessel is perfect. Also assume steady-state conditions.

a) T_1 = 938 K , T_2 = 931 K

b) To prevent softening of the materials, which would occur below their  melting points, the reactor should not be operated much above:

                                      q = 3*10^8 W/m^3

Explanation:

Given:

- See the attachment for the figure for this question.

- Melting point of Thorium T_th = 2000 K

- Melting point of Thorium T_g = 2300 K

Find:

a) If the thermal energy is uniformly generated in the fuel element at a rate q = 10^8 W/m^3 then what are the temperatures T_1 and T_2 at the inner and outer surfaces, respectively, of the fuel element?

b) Compute and plot the temperature distribution in the composite wall for selected values of q.  What is the maximum allowable value of q.

Solution:

part a)

- The outer surface temperature of the fuel, T_2, may be determined from the rate equation:

                                 q*A_th = T_2 - T_inf / R'_total

Where,

           A_th: Area of the thorium section

           T_inf: The temperature of coolant = 600 K

           R'_total: The resistance per unit length.

- Calculate the resistance per unit length R' from thorium surface to coolant:

           R'_total = Ln(r_3/r_2) / 2*pi*k_g + 1 / 2*pi*r_3*h

Plug in values:

           R'_total = Ln(14/11) / 2*pi*3 + 1 / 2*pi*0.014*2000

           R'_total = 0.0185 mK / W

- And the heat rate per unit length may be determined by applying an energy balance to a control surface  about the fuel element. Since the interior surface of the element is essentially adiabatic, it follows that:

           q' = q*A_th = q*pi*(r_2^2 - r_1^2)

           q' = 10^8*pi*(0.011^2 - 0.008^2) = 17,907 W / m

Hence,

           T_2 = q' * R'_total + T_inf

           T_2 = 17,907*0.0185 + 600

          T_2 = 931 K

- With zero heat flux at the inner surface of the fuel element, We will apply the derived results for boundary conditions as follows:

 T_1 = T_2 + (q*r_2^2/4*k_th)*( 1 - (r_1/r_2)^2) - (q*r_1^2/2*k_th)*Ln(r_2/r_1)

Plug values in:

 T_1 = 931+(10^8*0.011^2/4*57)*( 1 - (.8/1.1)^2) - (10^8*0.008^2/2*57)*Ln(1.1/.8)

 T_1 = 931 + 25 - 18 = 938 K

part b)

The temperature distributions may be obtained by using the IHT model for one-dimensional, steady state conduction in a hollow tube. For the fuel element (q > 0),  an adiabatic surface condition is  prescribed at r_1 while heat transfer from the outer surface at r_2 to the coolant is governed by the thermal  resistance:

                              R"_total = 2*pi*r_2*R'_total

                              R"_total = 2*pi*0.011*0.0185 = 0.00128 m^2K/W

- For the graphite ( q = 0 ), the value of T_2 obtained from the foregoing solution is prescribed as an inner boundary condition at r_2, while a convection condition is prescribed at the outer surface (r_3).

- For 5*10^8 < q and q > 5*10^8, the distributions are given in attachment.

The graphs obtained:

- The comparatively large value of k_t yields small temperature variations across the fuel element,  while the small value of k_g results in large temperature variations across the graphite.

Operation  at q = 5*10^8 W/^3  is clearly unacceptable, since the melting points of thorium and graphite are exceeded  and approached, respectively. To prevent softening of the materials, which would occur below their  melting points, the reactor should not be operated much above:

                                      q = 3*10^8 W/m^3

frequency or wavelenth of em radiation in a vacuum