‭MM 218 Programming Assignment‬

‭GROUP 40:‬
‭ ANSY JENA:‬‭22B2513‬
T
‭ASHUTOSH GANDHE:‬‭22B2409‬
‭PRATHAMESH YEOLE:‬‭22B2541‬

‭PROBLEM 4:‬
‭ onsider‬‭two‬‭plates,‬‭A‬‭and‬‭B,‬‭that‬‭are‬‭each‬‭initially‬‭isothermal‬‭and‬‭each‬‭of‬‭thickness‬‭L‬‭=‬‭5mm.‬‭The‬
C
‭faces‬ ‭of‬ ‭the‬ ‭plates‬ ‭are‬ ‭suddenly‬ ‭brought‬ ‭into‬ ‭contact‬ ‭in‬ ‭a‬ ‭joining‬ ‭process.‬ ‭Material‬ ‭A‬ ‭is‬ ‭acrylic,‬
‭initially‬‭at‬‭T‬‭i,A‬‭=‬‭20‬‭o‬‭C‬‭with‬‭ρ‭A‬ ‬ ‭=‬‭1990‬‭kg/m‬‭3‬‭,‬‭c‬‭A‬ ‭=‬‭1470‬‭J/kg.K‬‭and‬‭k‬‭A‬ ‭=‬‭0.21‬‭W/m.K.‬‭Material‬‭B‬‭is‬‭steel‬‭initially‬
‭at‬‭Ti,‬‭B‬ ‭=‬‭300‬‭o‬‭C‬‭with‬‭ρ‬‭B‬ ‭=‬‭7800‬‭kg/m‬‭3‭,‬ ‬‭c‬‭B‬ ‭=‬‭500‬‭J/kg.K‬‭and‬‭k‭B‬ ‬ ‭=‬‭45‬‭W/m.K.‬‭The‬‭external‬‭(back)‬‭surfaces‬‭of‬
‭the‬ ‭acrylic‬ ‭and‬ ‭steel‬ ‭are‬ ‭insulated.‬ ‭Neglecting‬ ‭the‬ ‭thermal‬ ‭contact‬ ‭resistance‬‭between‬‭the‬‭plates,‬
‭determine‬ ‭how‬ ‭long‬ ‭it‬ ‭will‬ ‭take‬ ‭for‬ ‭the‬ ‭external‬ ‭surface‬ ‭of‬ ‭the‬ ‭acrylic‬ ‭to‬ ‭reach‬ ‭its‬ ‭softening‬
‭temperature,‬ ‭T‬‭soft‬ ‭=‬ ‭90‬‭o‬‭C.‬ ‭Plot‬ ‭the‬ ‭acrylic’s‬ ‭external‬ ‭surface‬ ‭temperature‬ ‭as‬ ‭well‬ ‭as‬ ‭the‬ ‭average‬
‭temperatures of both materials over the time span . Use 20 equally spaced nodal points.‬
MATLAB CODE‬
‭

%we will make an array of 300/delt x 20 shape‬

‭

%the 20 columns will denote the temperature of each node‬

‭

%at the time instance denoted by the row number‬

‭

%rows denote time instances‬

‭

%we can set delt = 0.0001 seconds (Using the stablility condition)‬
‭

%Total time instances = 3000001‬

‭

%3000000 starting from t = 0.0 and the 1 extra for when we reach t=300‬
‭

T= zeros(3000001, 20);‬
‭

%array for calculating average temp‬

‭

T_a_avg = zeros(3000001, 1);‬

‭

T_b_avg = zeros(3000001, 1);‬

‭

%Constants‬
‭

delt = 0.0001;‬
‭

delx = 5 * 10^(-3)/9.5;‬
‭

%For A‬
‭

ka = 0.21;‬
‭

ca = 1470;‬
‭

pa = 1990;‬
‭

%Note: lambda = k/(c*p) * (delt/(delx)^2‬

‭

lam_a = ka/(ca*pa) * delt/(delx)^2;‬

‭

%For B‬
‭
kb = 45;‬
‭

cb = 500;‬
‭

pb = 7800;‬
‭

lam_b = kb/(cb*pb) * delt/(delx)^2;‬

‭

%resistance at the interface‬

‭

R = delx/(2*ka) + delx/(2*kb);‬
‭

%now for the initial temperatures of the two materials‬

‭

%nodes 1 to 10 belonged to material A and hence their‬

‭

%initial temp = 20 deg celcius‬

‭

T(:, 1:10) = 20+273.15;‬

‭

%similar for material B, from nodes 11 to 20‬

‭

%initial temp = 300 deg celcius‬

‭

T(:, 11:20) = 300+273.15;‬

‭

%Implement 1 round of gauss seidel below‬

‭

%going by time instances denoted by L‬

‭

%counter to print 90 degree celcius‬

‭

z=0;‬
‭

for‬‭
‭ L = 1:300000‬
‭0‬

%now to
‭ column by column to perform gauss seidel‬

%at each node‬

‭

for‬‭
‭ i=1:20‬

%for node 1‬
‭
 if‬‭
‭ i ==1‬

T(L+1,i) = 2*lam_a * (T(L, i+1) - T(L,i)) + T(L,i);‬

‭

end‬
‭

%for nodes 2 to 9‬
‭

if‬‭
‭ i >1 && i<10‬

T(L+1,i) = (lam_a * (T(L, i+1) + - 2*T(L,i) + T(L, i-1))+ T(L,i));‬

‭

end‬
‭

%for node 10‬

‭

if‬‭
‭ i == 10‬

T(L+1,i) = (lam_a*delx/(ka*R)) * (T(L, i+1) - T(L,i)) + lam_a* (T(L,‬

‭
i-1) - T(L,i))+ T(L,i);‬
‭

end‬
‭

%for node 11‬

‭

if‬‭
‭ i == 11‬

T(L+1,i) = (lam_b*delx/(kb*R)) * (T(L, i-1) - T(L,i)) + lam_b*(T(L,‬

‭
i+1) - T(L,i))+ T(L,i);‬
‭

end‬
‭

%for nodes 12 to 19‬

‭

if‬‭
‭ i >11 && i<20‬

T(L+1,i) = (lam_b * (T(L, i+1) + - 2*T(L,i) + T(L, i-1))+ T(L,i));‬

‭

end‬
‭

%for node 20‬

‭
 if‬‭
‭ i == 20‬

T(L+1,i) = 2*lam_b * (T(L, i-1) - T(L,i)) + T(L,i);‬

‭

end‬
‭

end‬
‭

%to print the instance at which‬

‭

%temperature of external acrylic surface becomes 90 degrees celcius‬

‭

%we set a tolerable value around which we wish to see‬

‭

%the temperature to reach nearly 90 degrees‬

‭

if‬‭
‭ (abs(T(L, 1) - (90+273)) < 0.001 && z == 0)‬

fprintf(‬
‭ "The time at which temperature is 90 degrees = %9.9f seconds \n"‬
‭ ,‬
‭
L * 0.0001);‬
‭

z=1;‬
‭

end‬
‭

end‬
‭

%calculating average temp‬

‭

for‬‭
‭ j = 1:3000001‬

%to calculate the sum‬

‭

sum_a = 0;‬
‭

sum_b = 0;‬
‭

for‬‭
‭ k = 1: 20‬
%for material A‬
‭

if‬‭
‭ (k<11)‬

%delx/2 for node 1‬

‭

if‬‭
‭ (k == 1)‬

sum_a = sum_a + 0.5*T(j,k);‬

‭

%using delx for node 2 to 10‬

‭

else‬
‭

sum_a = sum_a + T(j,k);‬

‭

end‬
‭

end‬
‭

if‬‭
‭ k>10‬

%using delx/2 for node 20‬

‭

if‬‭
‭ k == 20‬

sum_b = sum_b + 0.5*T(j,k);‬

‭

%using delx for node 11-19‬

‭

else‬
‭

sum_b = sum_b + T(j,k);‬

‭
 end‬
‭

end‬
‭

end‬
‭

T_a_avg(j) = sum_a/9.5;‬
‭

T_b_avg(j) = sum_b/9.5;‬
‭

end‬
‭

%now for plotting purposes‬

‭

x = linspace(1, 3000001, 3000001);‬

‭

time_range = x*0.0001;‬
‭

%we plot the surface temperature‬

‭

plot(time_range(1:3000000), T(1:3000000, 1));‬

‭

xlabel(‬
‭ 'time in seconds'‬
‭ );‬
‭

ylabel(‬
‭ 'Temperature (K)'‬
‭ );‬
‭

title(‬
‭ 'External Temperature of Acrylic'‬
‭ );‬
‭

grid‬‭
‭ on‬
;‬
‭

print(‬
‭ 'External Temp of Acrylic'‬
‭ ,‬‭
‭ '-dpng'‬
);‬
‭

figure;‬
‭
%now to plot the average surface temperatures‬
‭

%for material A‬
‭

plot(time_range(1:3000000), T_a_avg(1:3000000));‬
‭

hold‬‭
‭ on‬
;‬
‭

%for material B‬
‭

plot(time_range(1:3000000), T_b_avg(1:3000000));‬
‭

xlabel(‬
‭ 'Time in seconds'‬
‭ );‬
‭

ylabel(‬
‭ 'Temperature of Surfaces (K)'‬
‭ );‬
‭

legend(‬
‭ 'Material A'‬
‭ ,‬‭
‭ 'Material B'‬
);‬
‭

title(‬
‭ 'Average Temperatures of Materials'‬
‭ )‬
‭

grid‬‭
‭ on‬
;‬
‭

print(‬
‭ 'Average Surface Temp'‬
‭ ,‬‭
‭ '-dpng'‬
);‬
‭