Pre-Analysis: Example: Steady One-Dimensional Heat Conduction in A Bar
Pre-Analysis: Example: Steady One-Dimensional Heat Conduction in A Bar
Pre-Analysis: Example: Steady One-Dimensional Heat Conduction in A Bar
1. Mathematical model
2. Numerical solution procedure
3. Hand-calculations of expected results/trends
x
z
= 0,
Boundary conditions
0 =
q L = q =
Numerical Solution:
Discretization
Reduce the problem to determining temperature
values at selected locations (nodes)
T
1
x
We have assumed a shape for ( ) consisting of piecewise
polynomials
Mathematical Model
(Boundary Value
Problem)
Piecewise
polynomial
approximation for T
={ }
Each algebraic
equation will
relate a nodal
temperature to
its neighbors
Nodal
temperatures
Post
processing
( )
=0
System of algebraic
eqs. in nodal
temperatures
Piecewise polynomial
approximation for T
Weighted Integral
Form
=0
( ) is an
arbitrary function
=0
(x) is an
arbitrary piecewise
polynomial function
dx = 0
(x) is an
arbitrary piecewise
polynomial function
Integration by Parts
w k
=0
w k
dT
1
w
=0
=0
dT
+ 0.5
+ 0.5 +
=0
w k
dT
+w
w k
dT
=0
=0
+ 0.5 +
dT
k
3
dT
w k
1
w
dT
dT
=0
=0
3
+ 0.5
+ 0.5 +
=0
={ }
w k
1
dT
k
3
dT
={ }
dT
w k
1
w
+ 0.5
+
+ 0.5 +
dT
=0
=0
w k
dT
=0
k
2
dT
4
={ }
= 0.5
= 0.5 +
=T
= 0.5 +
= 0.5Q
=
=
5.5 W/m
Energy is
conserved for
the bar
10
Original Mesh
2
Refined Mesh
1
Second-Order Element
6 elements
11
12