summaryrefslogtreecommitdiff
path: root/587/CH2/EX2.6/example2_6.sce
blob: fcf80bc41f2a48734ff445db200fb4facd5ec589 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
clear;
clc;

//Example2.6[Heat Conduction in a Short Cylinder]
//Given:-
//Radius R and height h of the small cylinder
T=300;//Temperature of cylinder[degree Celcius]
T_ambient=20;//Temperature of ambient air[degree Celcius]
//Variation is thermal conductivity is negligible
//The cylinder is cooled unifromly from the top and bottom surfaces in the z-direction as well as the lateral surface in the radial r-direction.Also Temperature at any point in the ball changes with time during cooling. Therefore this is a two dimensional transient heat conduction problem i.e. T=T(r,z,t)
disp("The differential equation governing the variation of temperature in the billet is ")
disp("((1/r)(d/dr)(k*r(dT/dr)))+((d/dz)(k(dT/dz)))=rho*c(dT/dt)")
disp("In case of constant thermal conductivity it reduces to")
disp("((1/r)(d/dr)(r(dT/dr)))+(d^2)T/(dz^2)=(1/a)(dT/dt)")