ME-460 Solid State Joining
Transient problems with MATLAB
Transient 2D in plate
Lets consider a square plate as shown in fig.
The entire plate is maintained at a temp. of
250 C except the edges and temp. at the
edges are shown in the fig. Determine the
temp at the nodes considering a 13*13 grid
mesh.
250 C
We know 2D transient heat conduction
equation will be the governing equation for
the case.
Transient 2D in plate
So as problem has to be solved for transient state , then the generalized heat
conduction equation in 2 D for Cartesian coordinate system reduces to
𝜕2 𝑇 𝜕2 𝑇 1 𝜕𝑇
+ 2 = ……………..(1)
𝜕𝑥 2 𝜕𝑦 α 𝜕𝑡
where α is thermal diffusivity of material.
Transient 2D in plate
The discretized equation can be written as
𝜕2 𝑇 𝑇(𝑖−1,𝑗) −2𝑇(𝑖,𝑗) +𝑇(𝑖+1,𝑗)
By central difference =
𝜕𝑥 2 ∆𝑥 2
𝜕2 𝑇 𝑇(𝑖,𝑗−1) −2𝑇(𝑖,𝑗) +𝑇(𝑖,𝑗−1)
By central difference = at time t
𝜕𝑦 2 ∆𝑦 2
𝑡+1 𝑡
𝜕𝑇 𝑇𝑖,𝑗 − 𝑇𝑖,𝑗
By forward difference in time =
𝜕𝑡 ∆𝑡
Transient 2D in plate
250 C
Transient 2D in plate
Where 𝑇𝑖,𝑗
𝑡+1
= 𝑇𝑛 (𝑖, 𝑗)
Effect of iteration on Temp Values
Tolerance =10^(-1) , No. of iteration (K)=499
Effect of iteration on Temp Values
Tolerance =10^(-3) , No. of iteration (K)=1201
Temp. profile at time-step t=1000
Temp Values at time-step t=1000
Temp. profile at time-step t=2000
Temp Values at time-step t=2000
Temp. profile at time-step t=5000
Temp Values at time-step t=5000
Thank you!!