Consider a large uranium plate of thickness \(5 \mathrm{~cm}\) and thermal conductivity \(k=28 \mathrm{~W} / \mathrm{m} \cdot \mathrm{K}\) in which heat is generated uniformly at a constant rate of \(\dot{e}=6 \times 10^{5} \mathrm{~W} / \mathrm{m}^{3}\). One side of the plate is insulated while the other side is subjected to convection to an environment at \(30^{\circ} \mathrm{C}\) with a heat transfer coefficient of \(h=60 \mathrm{~W} / \mathrm{m}^{2} \cdot \mathrm{K}\). Considering six equally spaced nodes with a nodal spacing of \(1 \mathrm{~cm},(a)\) obtain the finite difference formulation of this problem and \((b)\) determine the nodal temperatures under steady conditions by solving those equations.

Short Answer

Expert verified
Answer: The steady-state nodal temperatures for the uranium plate under the given conditions are as follows: \(T_1 \approx 198.8^{\circ} \mathrm{C}\) \(T_2 \approx 198.8^{\circ} \mathrm{C}\) \(T_3 \approx 163.2^{\circ} \mathrm{C}\) \(T_4 \approx 130.4^{\circ} \mathrm{C}\) \(T_5 \approx 104.4^{\circ} \mathrm{C}\) \(T_6 \approx 83.4^{\circ} \mathrm{C}\)

Step by step solution

01

Derive the 1D Heat Conduction Equation

Start by writing the 1-D heat conduction equation for the system, which is given by: \(\frac{d}{dx}\left(k\frac{dT}{dx}\right)+\dot{e}=0\) Here, \(T\) is the temperature, \(k\) is the thermal conductivity and \(\dot{e}\) is the heat generation rate per unit volume.
02

Discretization of the 1D Heat Conduction Equation

Discretize the 1-D heat conduction equation using finite difference methods. Divide the 5 cm thickness into six equally spaced nodes with a nodal spacing of 1 cm. For the case of heat generation, the finite difference equation can be derived as: \(k\frac{T_{i+1}-2T_i+T_{i-1}}{\Delta x^2}+\dot{e}=0\), for i = 2, 3, 4 and 5. We have \(k = 28 \mathrm{~W} / \mathrm{m} \cdot \mathrm{K}\), \(\dot{e} = 6 \times 10^5 \mathrm{~W} / \mathrm{m}^3\), and \(\Delta x = 1 \mathrm{~cm}= 0.01 \mathrm{~m}\).
03

Apply Boundary Conditions

Now, apply the boundary conditions for the given problem. At node 1, the side of the plate is insulated, which implies that no heat flow takes place at this node: \(\frac{dT}{dx}_{x=1} = 0\), or equivalently, \(T_1 = T_2\). At node 6, the other side of the plate is subjected to convection: \(h(T_6 - T_{\infty}) = -k\frac{dT}{dx}_{x=6}\) This can be discretized as: \(h(T_6 - T_{\infty}) = -k\frac{T_6 - T_5}{\Delta x}\). Here, \(h = 60 \mathrm{~W} / \mathrm{m}^{2} \cdot \mathrm{K}\), and \(T_{\infty} = 30^{\circ} \mathrm{C}\).
04

Solve the System of Linear Equations

Next, construct a system of linear equations for nodal temperatures using the discretized heat conduction equation (from Step 2) and the boundary conditions (from Step 3). This results in the following set of equations: \(T_1 - T_2 = 0\) \(28(T_2 - 2T_1 + T_3) + 6\times10^5(0.01)^2 = 0\) \(28(T_3 - 2T_2 + T_1) + 6\times10^5(0.01)^2 = 0\) \(28(T_4 - 2T_3 + T_2) + 6\times10^5(0.01)^2 = 0\) \(28(T_5 - 2T_4 + T_3) + 6\times10^5(0.01)^2 = 0\) \(60(T_6 - 30) + 28(T_6 - T_5) = 0\) Now, we need to solve this linear system to find the nodal temperatures \(T_1, T_2, ..., T_6\). This can be done by using standard linear algebra techniques such as Gaussian elimination, matrix inversion, or iterative methods such as Gauss-Seidel.
05

Calculation of Nodal Temperatures

By solving the system of linear equations obtained above, we get the following nodal temperatures under steady-state conditions: \(T_1 \approx 198.8^{\circ} \mathrm{C}\) \(T_2 \approx 198.8^{\circ} \mathrm{C}\) \(T_3 \approx 163.2^{\circ} \mathrm{C}\) \(T_4 \approx 130.4^{\circ} \mathrm{C}\) \(T_5 \approx 104.4^{\circ} \mathrm{C}\) \(T_6 \approx 83.4^{\circ} \mathrm{C}\) These are the steady-state nodal temperatures for the uranium plate under the given conditions.

Unlock Step-by-Step Solutions & Ace Your Exams!

  • Full Textbook Solutions

    Get detailed explanations and key concepts

  • Unlimited Al creation

    Al flashcards, explanations, exams and more...

  • Ads-free access

    To over 500 millions flashcards

  • Money-back guarantee

    We refund you if you fail your exam.

Over 30 million students worldwide already upgrade their learning with Vaia!

One App. One Place for Learning.

All the tools & learning materials you need for study success - in one app.

Get started for free

Most popular questions from this chapter

What are the basic steps involved in solving a system of equations with Gauss- Seidel method?

Consider a heat conduction problem that can be solved both analytically, by solving the governing differential equation and applying the boundary conditions, and numerically, by a software package available on your computer. Which approach would you use to solve this problem? Explain your reasoning.

Consider a long concrete dam \((k=0.6 \mathrm{~W} / \mathrm{m} \cdot \mathrm{K}\), (es \(\alpha_{s}=0.7\) ) of triangular cross section whose exposed surface is subjected to solar heat flux of \(\dot{q}_{s}=\) \(800 \mathrm{~W} / \mathrm{m}^{2}\) and to convection and radiation to the environment at \(25^{\circ} \mathrm{C}\) with a combined heat transfer coefficient of \(30 \mathrm{~W} / \mathrm{m}^{2} \cdot \mathrm{K}\). The \(2-\mathrm{m}\)-high vertical section of the dam is subjected to convection by water at \(15^{\circ} \mathrm{C}\) with a heat transfer coefficient of \(150 \mathrm{~W} / \mathrm{m}^{2} \cdot \mathrm{K}\), and heat transfer through the 2-m-long base is considered to be negligible. Using the finite difference method with a mesh size of \(\Delta x=\Delta y=1 \mathrm{~m}\) and assuming steady two-dimensional heat transfer, determine the temperature of the top, middle, and bottom of the exposed surface of the dam.

Write a two-page essay on the finite element method, and explain why it is used in most commercial engineering software packages. Also explain how it compares to the finite difference method.

Consider an aluminum alloy fin \((k=180 \mathrm{~W} / \mathrm{m} \cdot \mathrm{K})\) of triangular cross section whose length is \(L=5 \mathrm{~cm}\), base thickness is \(b=1 \mathrm{~cm}\), and width \(w\) in the direction normal to the plane of paper is very large. The base of the fin is maintained at a temperature of \(T_{0}=180^{\circ} \mathrm{C}\). The fin is losing heat by convection to the ambient air at \(T_{\infty}=25^{\circ} \mathrm{C}\) with a heat transfer coefficient of \(h=25 \mathrm{~W} / \mathrm{m}^{2} \cdot \mathrm{K}\) and by radiation to the surrounding surfaces at an average temperature of \(T_{\text {surr }}=290 \mathrm{~K}\). Using the finite difference method with six equally spaced nodes along the fin in the \(x\)-direction, determine \((a)\) the temperatures at the nodes and \((b)\) the rate of heat transfer from the fin for \(w=1 \mathrm{~m}\). Take the emissivity of the fin surface to be \(0.9\) and assume steady one-dimensional heat transfer in the fin.

See all solutions

Recommended explanations on Physics Textbooks

View all explanations

What do you think about this solution?

We value your feedback to improve our textbook solutions.

Study anywhere. Anytime. Across all devices.

Sign-up for free