Quench hardening is a mechanical process in which the ferrous metals or alloys are first heated and then quickly cooled down to improve their physical properties and avoid phase transformation. Consider a \(40 \mathrm{~cm} \times 20 \mathrm{~cm}\) block of copper alloy \(\left(k=120 \mathrm{~W} / \mathrm{m} \cdot \mathrm{K}, \alpha=3.91 \times 10^{-6} \mathrm{~m}^{2} / \mathrm{s}\right)\) being heated uniformly until it reaches a temperature of \(800^{\circ} \mathrm{C}\). It is then suddenly immersed into the water bath maintained at \(15^{\circ} \mathrm{C}\) with \(h=100 \mathrm{~W} / \mathrm{m}^{2} \cdot \mathrm{K}\) for quenching process. However, the upper surface of the metal is not submerged in the water and is exposed to air at \(15^{\circ} \mathrm{C}\) with a convective heat transfer coefficient of \(10 \mathrm{~W} / \mathrm{m}^{2} \cdot \mathrm{K}\). Using an explicit finite difference formulation, calculate the temperature distribution in the copper alloy block after \(10 \mathrm{~min}\) have elapsed using \(\Delta t=10 \mathrm{~s}\) and a uniform mesh size of \(\Delta x=\Delta y=10 \mathrm{~cm}\).

Short Answer

Expert verified
In this exercise, we have predicted the temperature distribution of a copper alloy block after quenching it in water for 10 minutes. We used the explicit finite difference method and solved the heat equation with given boundary conditions. The copper block was divided into a mesh, and calculations were made using the given time step of 10 seconds to find the temperature distribution after 10 minutes. Note that the specific temperature values at each location depend on the given material properties and may require further calculation or numerical simulation.

Step by step solution

01

Setting Up the Math Model

We have a 2D transient conduction problem governed by the heat equation: \(\frac{\partial T}{\partial t} = \alpha \cdot \left(\frac{\partial^2 T}{\partial x^2} + \frac{\partial^2 T}{\partial y^2}\right)\) Given additional information: - Block size: \(40 \mathrm{~cm} \times 20 \mathrm{~cm}\) - Initial temperature: \(800^{\circ} \mathrm{C}\) - Convective heat transfer coefficients: \(h_{water} = 100 \mathrm{~W} / \mathrm{m}^{2} \cdot \mathrm{K}\) and \(h_{air} = 10 \mathrm{~W} / \mathrm{m}^{2} \cdot \mathrm{K}\) - Mesh size: \(\Delta x = \Delta y = 10 \mathrm{~cm}\)
02

Apply Explicit Finite Difference Method

The explicit finite difference method approximates the spatial derivatives using centered differences: \(\frac{\partial^2 T}{\partial x^2} \approx \frac{T_{i+1,j}^{n} - 2T_{i,j}^{n} + T_{i-1,j}^{n}}{(\Delta x)^2}\) \(\frac{\partial^2 T}{\partial y^2} \approx \frac{T_{i,j+1}^{n} - 2T_{i,j}^{n} + T_{i,j-1}^{n}}{(\Delta y)^2}\) Substituting these into the heat equation, we get: \(\frac{T_{i,j}^{n+1} - T_{i,j}^{n}}{\Delta t} = \alpha \cdot \left(\frac{T_{i+1,j}^{n} - 2T_{i,j}^{n} + T_{i-1,j}^{n}}{(\Delta x)^2} + \frac{T_{i,j+1}^{n} - 2T_{i,j}^{n} + T_{i,j-1}^{n}}{(\Delta y)^2}\right)\) Rearrange the equation to get the value of temperature at the next time step: \(T_{i,j}^{n+1} = T_{i,j}^{n} + \frac{\alpha \Delta t}{(\Delta x)^2}\left(T_{i+1,j}^{n} - 2T_{i,j}^{n} + T_{i-1,j}^{n}\right) + \frac{\alpha \Delta t}{(\Delta y)^2}\left(T_{i,j+1}^{n} - 2T_{i,j}^{n} + T_{i,j-1}^{n}\right)\)
03

Implement Method with Given Time and Space Steps

- Divide the copper block into a mesh of \(4 \times 2\) cells, with each cell having dimensions \(10 \mathrm{~cm} \times 10 \mathrm{~cm}\). - Initialize all cells to the temperature \(800^{\circ} \mathrm{C}\) and apply the given boundary conditions. - Use the given time step \(\Delta t = 10 \mathrm{~s}\) and perform the calculation for \(10 \mathrm{~min}\) or \(600 \mathrm{~s}\). This means we need to perform the calculation for \(n = 60\) time steps.
04

Compute Temperature Distribution

At each time step, update the temperature in each cell of the mesh using the explicit finite difference method equation from Step 2. Apply boundary conditions on the submerged and exposed surfaces. After 60 time steps or \(10 \mathrm{~min}\), the temperature distribution in the copper alloy block will have been calculated. Note that the explicit finite difference method may have stability issues if the time step is too large. In this case, it is assumed that the given time step is appropriate for this exercise.

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

Using EES (or other) software, solve these systems of algebraic equations. (a) $$ \begin{array}{r} 3 x_{1}-x_{2}+3 x_{3}=0 \\ -x_{1}+2 x_{2}+x_{3}=3 \\ 2 x_{1}-x_{2}-x_{3}=2 \end{array} $$ (b) $$ \begin{aligned} 4 x_{1}-2 x_{2}^{2}+0.5 x_{3} &=-2 \\ x_{1}^{3}-x_{2}+x_{3} &=11.964 \\ x_{1}+x_{2}+x_{3} &=3 \end{aligned} $$

Starting with an energy balance on a volume element, obtain the steady two- dimensional finite difference equation for a general interior node in rectangular coordinates for \(T(x, y)\) for the case of variable thermal conductivity and uniform heat generation.

Suggest some practical ways of reducing the roundoff error.

Consider a 2-m-long and 0.7-m-wide stainless-steel plate whose thickness is \(0.1 \mathrm{~m}\). The left surface of the plate is exposed to a uniform heat flux of \(2000 \mathrm{~W} / \mathrm{m}^{2}\) while the right surface of the plate is exposed to a convective environment at \(0^{\circ} \mathrm{C}\) with \(h=400 \mathrm{~W} / \mathrm{m}^{2} \cdot \mathrm{K}\). The thermal conductivity of the stainless steel plate can be assumed to vary linearly with temperature range as \(k(T)=k_{o}(1+\beta T)\) where \(k_{o}=48 \mathrm{~W} / \mathrm{m} \cdot \mathrm{K}\) and \(\beta=9.21 \times 10^{-4}{ }^{\circ} \mathrm{C}^{-1}\). The stainless steel plate experiences a uniform volumetric heat generation at a rate of \(8 \times 10^{5} \mathrm{~W} / \mathrm{m}^{3}\). Assuming steady state one-dimensional heat transfer, determine the temperature distribution along the plate thickness.

A circular fin \((k=240 \mathrm{~W} / \mathrm{m} \cdot \mathrm{K})\) of uniform cross section, with diameter of \(10 \mathrm{~mm}\) and length of \(50 \mathrm{~mm}\), is attached to a wall with surface temperature of \(350^{\circ} \mathrm{C}\). The fin tip has a temperature of \(200^{\circ} \mathrm{C}\), and it is exposed to ambient air condition of \(25^{\circ} \mathrm{C}\) and the convection heat transfer coefficient is \(250 \mathrm{~W} / \mathrm{m}^{2} \cdot \mathrm{K}\). Assume steady one-dimensional heat transfer along the fin and the nodal spacing to be uniformly \(10 \mathrm{~mm}\), (a) using the energy balance approach, obtain the finite difference equations to determine the nodal temperatures, and (b) determine the nodal temperatures along the fin by solving those equations and compare the results with the analytical solution.

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