Heat Transfer Analysis

Fares Al Rwaihy

May 16, 2023

12

0

This series of blog posts aims to compare the hand calculation and numerical simulation techniques through technical softwares

4 min read

Steady State Thermal Analysis

Main Image

Introduction:

Heat transfer is an important aspect of many engineering applications, ranging from electronic cooling to thermal management in buildings. Accurately predicting heat transfer rates and temperatures is crucial for optimizing system performance and ensuring reliability.

Objectives:

The main objective of this series of blog posts is to compare hand calculation and numerical simulation techniques in heat transfer. We will discuss the role of technical software, such as finite element analysis (FEA), in accurately simulating complex heat transfer phenomena. By the end of this blog, readers will have a deeper understanding of heat transfer and the techniques used to analyze and optimize heat transfer systems.

Overview

Heat transfer is the exchange of thermal energy between physicals systems. A physical system is a region of space bounded by a closed surface and filled with a physical substance (medium) which could be a solid, liquid, or gaseous.

However, the amount of heat transferred through a cross section are (A) at moment (t) is defined as following:

IQ(t) = lim⁡Δt→0(ΔQΔt) [W]I_Q\left(t\right)\:=\:\lim _{\Delta t\to 0}\left(\frac{\Delta Q}{\Delta t}\right)\:\left[W\right]

In case of Steady-State heat transfer the amount of heat intensity remains constant along the cross section are (A). Furthermore in order for the heat to be transferred, three modes of transforming which are Conduction, Convection, and Radiation.

**Conduction and Convection can be only released, if the transfer space is filled with a medium which is not necessary for radiation which could be transferred in a vacuum

Steady-State thermal conduction

Thermal conduction equation

Heat is transferred from a hotter region of the intervening medium to a colder one through collisions between neighboring objects such as molecules, atoms, or electrons.

Image

It follows this equation

IQ = −k⋅ A⋅ T2−T1xI_Q\:=\:-k\cdot \:A\cdot \:\frac{T2-T1}{x}

where k is the thermal conductivity of the intervening medium, thus the thermal resistance R can be introduced

R = xk⋅A thus===> IQ ⋅ R= −(T2−T1) R\:=\:\frac{x}{k\cdot A}\:thus===>\:I_Q\:\cdot \:R=\:-\left(T2-T1\right)\:

Thermal resistance in series vs parallel

in Series:

ΔT = ΔT1+ ΔT2\Delta T\:=\:\Delta T_1+\:\Delta T_2
IQeq⋅  Req = IQ1⋅  R1 + IQ2⋅  R2I_{Qeq}\cdot \:\:R_{eq}\:=\:I_{Q1}\cdot \:\:R_1\:+\:I_{Q2}\cdot \:\:R_2
Using  IQeq = IQ1= IQ2Using\:\:I_{Qeq}\:=\:I_{Q1}=\:I_{Q2}
Req = R1+R2R_{eq}\:=\:R_1+R_2

in Parallel:

IQeq = IQ1+IQ2I_{Qeq}\:=\:I_{Q1}+I_{Q2}
ΔTeqReq = ΔT1R1+ΔT2R2\frac{\Delta T_{eq}}{R_{eq}}\:=\:\frac{\Delta T_1}{R_1}+\frac{\Delta T_2}{R_2}
Using  ΔTeq = ΔT1=ΔT2 Using\:\:\Delta T_{eq}\:=\:\Delta T_1=\Delta T_2\:
1Req = 1R1+1R2\frac{1}{R_{eq}}\:=\:\frac{1}{R_1}+\frac{1}{R_2}

Steady-State thermal convection

On the other hand, convection is not the transfer of heat within the layers of a material, but rather the transfer of heat between a surface and the surrounding medium through fluid flow.

It follows this equation

IQ = −α⋅A⋅(Tm−T)I_Q\:=\:-\alpha \cdot A\cdot \left(Tm-T\right)

where α is the convection heat transform coefficient of the intervening medium, thus the thermal resistance R can be introduced

Rα = 1α⋅ A thus===> IQ ⋅  Rα= −(Tm−T) R_{\alpha }\:=\:\frac{1}{\alpha \cdot \:A}\:thus===>\:I_Q\:\cdot \:\:R_{\alpha }=\:-\left(Tm-T\right)\:

Steady-State thermal FEA simulation

Problem Setup

The simulation below shows the composite wall of a cold store. The outer layer of the wall is made of structural steel while the inner one is of Aluminum.

Results

To [c]Ti [c]αo [W/k . m^2]αi [W/k . m^2]λ1 [W/m . k]λ2 [W/m . k]A [m^2]x [m]
25 -2 20 1060.5 237.5 100 1

As a result of the analysis, a chart is provided that indicates the variation of temperature with the distance along the wall layers.

Hand Calculation

Total thermal resistance

Req = Rα o+R1+R2+R3+ Rα iR_{eq}\:=\:R_{\alpha \:o}+R_1+R_2+R_3+\:R_{\alpha \:i}
1α0⋅ A+xλ1⋅ A+xλ2⋅A+xλ1⋅A+1α1⋅A\frac{1}{\alpha _0\cdot \:A}+\frac{x}{\lambda _1\cdot \:A}+\frac{x}{\lambda _2\cdot A}+\frac{x}{\lambda _1\cdot A}+\frac{1}{\alpha _1\cdot A}
120⋅100+160.5⋅100+1237.5⋅100+160.5⋅100+110⋅100 = 0.00187 [W]\frac{1}{20\cdot 100}+\frac{1}{60.5\cdot 100}+\frac{1}{237.5\cdot 100}+\frac{1}{60.5\cdot 100}+\frac{1}{10\cdot 100}\:=\:0.00187\:\left[W\right]

Heat intensity through the wall

IQ⋅Req = − (Ti−To)I_Q\cdot R_{eq}\:=\:-\:\left(T_i-To\right)
IQ = − (−2−25)0.00187 = 14417.84 [W]I_Q\:=\:\frac{-\:\left(-2-25\right)}{0.00187}\:=\:14417.84\:\left[W\right]

Calculate T1, T2, T3 surface temperature

we will first calculate the temperature drop between the ambient air and the surface of the composite wall, from T0 to T1 then T1 to T2....etc

IQ⋅ Rαo = − (T1−T0) = T0−T1I_Q\cdot \:R_{\alpha o}\:=\:-\:\left(T_1-T_0\right)\:=\:T_0-T_1
T1 = T0−IQ⋅ Rα 0 =17.79 [C∘]T_1\:=\:T_0-I_Q\cdot \:R_{\alpha \:0}\:=17.79\:\left[C^{\circ }\right]

Using the same principle

T2 = T1−IQ⋅ R1 =15.41[C∘]T_2\:=\:T_1-I_Q\cdot \:R_1\:=15.41\left[C^{\circ }\right]
T3 = T2−IQ⋅R2 =14.8 [C∘]T_3\:=\:T_2-I_Q\cdot R_2\:=14.8\:\left[C^{\circ }\right]
T4 = T3−IQ⋅R3 =12.42 [C∘]T_4\:=\:T_3-I_Q\cdot R_3\:=12.42\:\left[C^{\circ }\right]

Discussion

The agreement between the software simulation and hand calculation in predicting temperatures is a significant achievement. It validates the accuracy and reliability of both methods, highlighting the effectiveness of numerical techniques in modeling heat transfer. This convergence of results enhances our understanding of heat transfer mechanisms and supports the practical use of numerical simulations in engineering applications.