This preview has intentionally blurred sections. Sign up to view the full version.
View Full DocumentThis preview has intentionally blurred sections. Sign up to view the full version.
View Full Document
Unformatted text preview: City University of Hong Kong Electronic Engineering Semester B 20072008 EE3108 Engineering Analysis Tutorial 9 The cooling a copper rod is investigated using the heat equation. The rod is uniformed heated to 100°C before t =0. It is then put into contact with two cold objects at 0 °C. The problem is very similar to that in the lecture notes, but the boundary conditions are changed. The implicit method for solving the partial differential equation (PDE) is considered as follows. where Q(x,t) is again the temperature in Celsius and D = 1.16 cm2/s is the heat diffusion constant of copper. ∂Q ∂ 2Q =D 2 , ∂t ∂x t=0 x Cold contact L/2 t>0 t= 0 Hot Cooling −L/2 Cold contact The boundary conditions are modified as: Q x=± L / 2 = 0 .
The initial condition is Q( x,0) = 100 for −L/2< x <L/2. The length L is 10 cm. We are again interested in the time window from t = 0 to T, where T=10 s. The time axis is discretized into N−1 = 5000 steps and the length axis is discretized into M−1 = 100 steps so that: tn = (n − 1)∆t ∆t = T /( N − 1) xm = (m − 1)∆x − L / 2 ∆x = L /( M − 1)
m Qn = Q( xm , tn ) for n = 1, 2, …, N and m = 1, 2, …, M. EE 3108 Semster B 2007/2008 8
Equilibrium T91 S C Chan I. Propagation From the lecture notes, the difference equation is:
m m m m Qn = −bQn ++1 + (1 + 2b)Qn − bQn +−1 1 1 where b = D∆t /( ∆x ) . However, for m = 1 and M, the equations have to be modified to fit the boundary conditions. How?
2 In matrix notation, write down the propagation from n to n+1: T92 II. Matlab Modify the file below to implement the simulation. %File: t0902 clear; L=10; T=10; M=101; N=5001; dx=L/(M1); dt=T/(N1); x=((1:M)1)*dxL/2; t=((1:N)1)*dt; D=1.16; %Initialize the input at t=0 Q=zeros(M,N); %Find the propagation matrix S b=D*dt/(dx^2) S=zeros(M,M); for r=1:M if r==1 elseif r==M else end end %Propagate in time for n=1:N1 end %Plot mesh(t,x,Q); T93 T94 ...
View
Full
Document
This note was uploaded on 01/11/2011 for the course EE 3108 taught by Professor Nelsonszechunchan during the Spring '07 term at City University of Hong Kong.
 Spring '07
 NelsonSzeChunChan
 Gate

Click to edit the document details