Scribd
Upload
138 views

Finite Element - Heat Transfer

The document discusses 1-D heat transfer using finite element analysis. It will analyze heat transfer along one dimension, both with and without convection on perimeter surfaces. The analysis will discretize the governing differential equations using Galerkin's method. It provides the governing heat transfer equation and examines boundary conditions for a sample problem involving heat generation and convection in a rod.

Uploaded by

OpenTalk
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
138 views

Finite Element - Heat Transfer

The document discusses 1-D heat transfer using finite element analysis. It will analyze heat transfer along one dimension, both with and without convection on perimeter surfaces. The analysis will discretize the governing differential equations using Galerkin's method. It provides the governing heat transfer equation and examines boundary conditions for a sample problem involving heat generation and convection in a rod.

Uploaded by

OpenTalk
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 13

Lecture 22: 1-D Heat Transfer.

Chapter 13: Heat Transfer and Mass Transport.


We will focus initially on the steady state heat transfer problem. Start by looking at the transfer of thermal energy along one dimension.

With no convection off of the perimeter surface (insulated). With convection off of the perimeter surface.

Our analysis will be based on a discretization of the governing DEs using Galerkins method.

13.4 of Logan applies a variational formulation of of the heat transfer element equations. 13.8 switches to Galerkins method due to the addition of mass transport terms that complicate the variational method. We still use Steps 1 and 2 as illustrated in 13.4. If you read the rest of the 13.4 you will see why we are using a weighted residual technique.
MECH 420: Finite Element Applications

Lecture 22: 1-D Heat Transfer.


13.1. Derivation of the Basic Differential Equation.


First problem addressed is 1-D Heat Conduction with no convection . The finite element is a region in space.

The boundaries of the region are defined by fixed points (or nodes). There are quantities of interest at the boundaries of the region temperatures and heat flux values. Want to calculate the values for this discrete set of entities (the state variables). For 1-D problems this region can be considered as a filament (a line element). The fluctuations of the state variables along the line are assumed to be representative of what occurs over the whole region.

MECH 420: Finite Element Applications

Lecture 22: 1-D Heat Transfer.


Fouriers law of heat conduction:


dT qx = K xx dx

qx heat conducted (heat flux ) kW/m 2 K xx thermal conductivity in x direction kW/(o C m) dT temperature gradient o C/m dx c specific heat capacity kWh/(kg o C) E energy kWh U total internal thermal energy kWh Q rate of internal heat generation kW/m3
MECH 420: Finite Element Applications

Lecture 22: 1-D Heat Transfer.


The 1st law of thermodynamics states that there is a balance of the energy produced, stored, and conducted through the finite element (control volume):
Ein + QTOT = U + Eout

qx Adt + QAdxdt = U + qx + dx Adt dT dx dT d dT qx + dx = K xx + K xx dx dx dx dx U = c ( Adx ) dT qx = K xx

MECH 420: Finite Element Applications

Lecture 22: 1-D Heat Transfer.


Substituting these expressions into the energy balance


d dT dT K Q c + = xx dx dx dt

If the temperature profile can be assumed to be steady:


d dT K xx +Q = 0 dx dx

MECH 420: Finite Element Applications

Lecture 22: 1-D Heat Transfer.


13.3. Typical Units; thermal Conductivities, K; and Heat Transfer Coefficients, h.


Tables 13-1 through 13-3 list common K and h values and define the units of various thermal quantities. We will visit the use of the convective coefficients when we visit 13.2 in Lecture 23. The convective coefficients listed are for free convection. In forced convection there is a velocity of the liquid/gas created by an external source which would increase the effective h value.

MECH 420: Finite Element Applications

Lecture 22: 1-D Heat Transfer.


The first law of thermodynamics has produced a differential equation that defines a temperature distribution. The solution is a function T(x): the temperature distribution along the dimension x. The initial stage of the problem statement reduced us to a 1-D DE. We assumed there was heat transfer in only the x direction. Calculation of the true solution to the differential equation will involve applying boundary conditions to solve for constants of integration. We will look at the FE approximation to the true solution. We expect to have boundary terms exposed in our formulation of the element equations.
MECH 420: Finite Element Applications

Lecture 22: 1-D Heat Transfer.


Consider what the boundary terms might be

Boundary Conditions (case study):


At the left edge of the domain:


T = TB

Insulated boundaries ensuring a 1-D heat conduction.

At the right edge of the domain:


qx = K xx dT = q* x dx
q* x heat flux into element

MECH 420: Finite Element Applications

Lecture 22: 1-D Heat Transfer.


P. 13.4 (modified), pg.# 500.

The rod of 1 in. radius shown below generates heat internally at the rate of Q=10,000 Btu/(hft3) throughout the rod. The left edge and perimeter of the the rod are insulated, and the right edge is exposed to an ambient temperature of 100 deg F. The thermal conductivity of the rod is Kxx=12 Btu/(h-ft-degF) and the length of the rod is L=3 in. Calculate the temperature distribution along the rod axis assuming that there is no convection over the perimeter of the rod and that the convection off of the right end is ideal. Use at least three elements.

TB = T = 100 F

MECH 420: Finite Element Applications

Lecture 22: 1-D Heat Transfer.


Governing differential equation:


d dT K xx +Q = 0 dx dx

Fouriers law defines the heat flux:


qx = K xx dT dx

Provided an approximate temperature distribution we can recover the heat flux. A conforming element requires at least a 1st order approximation to the temperature distribution.

MECH 420: Finite Element Applications

Lecture 22: 1-D Heat Transfer.


Consider the boundary conditions: At the left end of the rod (flux condition):
dT =0 qx x =0 = K xx dx x =0 At the right end of the rod:

T x =3.0 = 100

x
x=0
MECH 420: Finite Element Applications

x=3

Lecture 22: 1-D Heat Transfer.


13.4 One dimensional Finite Element Formulation (using a WR method). We can use the weighted residuals method to get element equations that define a linear segment of an approximate temperature profile.
T2 T ( x) T1

x
TB = T = 100 F

x
x=0 x=3

MECH 420: Finite Element Applications

x = x1

x = x2

Lecture 22: 1-D Heat Transfer.


This problem is very similar to that completed in Lecture #16.


= NU U
2.4

x x N = 1 L L 12 U U = 1 U 2

2.3

2.2

x
L
( x) U

2.1

U(x)

= x x1 x + x1 x=x dx =1 dx

1.9

1.8

1.7

1.2

1.4

1.6

1.8

2.2

2.4

2.6

2.8

x1

x2

MECH 420: Finite Element Applications

You might also like