1D Transient Heat Conduction Problem in Cylindrical Coordinates Using FTCS Finite Difference Method

[@babeksultanov7562]

% Diffusion equation - ID - Explicit Method - Cylindrical I Polar % Coordinates % Dirichlet BC conditions - Constant temperatures at boundaries clear clc clf % Inputs alpha = 1e-4; % Thermal diffiisivity of copper, m2/s (Heat Conduction - Ozisik & Hahn) r_in = 0.1; % Inside Radius of hollow cylinder, r_in, say m" r_out = 1.1; % Outside Radius of hollow cylinder, r_out, say m" m = 10; % no. of sections divided between r_in and r_out eg 5, I0 delta_r = (r_out - r_in)/m; % section length, m" nr = m+1; % total no. ofradial points % phi = linspace(O, 2*pi, 500);" % total angle = 2*pi nphi = 36; % no. of angle steps delta_phi = 2*pi/nphi; % angle step, rad t = 9000; % total time, s eg 200, 9000 nt = 300; % no. of time steps eg 2, 300 delta_t = t/nt; % timestep, s" Trin = 200; % BCl, Temperature at r_in, deg.C Trout = 300; % BC2, Temperature at r_out, deg.C" Tin = 100; % IC, Initial Temperature at r_in

[@gabrielbetancourt5657]

such a master explanation, thanks a lot

[@shengyueshan3593]

Hello Sam, I am really interested in how to make the plots that you showed in this video, would you mind giving me some instructions on making the plots? Thank you very much!

[@muhammadreshaeel3061]

How to model such a problem if you want it for phase change material?

[@vivekananddhumal7853]

Hi Sam.. thank you very informative videos...can you make a video for 2d cylindrical coordinates?

[@samseleven1]

I think I already made one. Check into that. That is a heart sign profile. I think there is one more in my website

[@vivekananddhumal7853]

@@samseleven1 may I have the link

[@samseleven1]

@@vivekananddhumal7853 engineering-stream.com under Numerical methods

[@samseleven1]

I posted the code in engineering-stream.com under Numerical Methods Section. Hope that helps.

[@Saliqo911]

I could'nt find numerical section on your site sir. Can you please send a direct link to that section?

[@jingyuan2990]

only PPT slides, no code is found in the link

[@samseleven1]

@@jingyuan2990 check engineering-stream.com under Numerical Methods codes

[@samseleven1]

Codes may be obtained from engineering-stream.com under Numerical Methods section

[@هانيايوب-ض8د]

Can you help me with contour

[@ariespacudan9243]

Would it be possible to get a copy of the matlab code? please and thank you!

[@samseleven1]

See below

[@wholemonkey]

Is it possible to see the matlab code? It will be very helpful to understand the graphical plot, but I cannot find it

[@samseleven1]

Look for numerical methods at the website providef

[@samseleven1]

Look for Numerical Methods in the website provided.

[@wholemonkey]

@@samseleven1 can you please provide me the matlab code?

[@wholemonkey]

@@samseleven1 the matlab code to put the matrix in the cylindrical coordinate?

[@luzzyrogue]

What if r=0 at the centre? What do we get?

[@samseleven1]

Temp at the center is finite. So dT/dr =0 at r=0

[@luzzyrogue]

@@samseleven1 but in your equation we have 1/r

[@samseleven1]

@@luzzyrogue that is correct. The equation will blow up if r=0. One way to tackle this problem is to keep r closer to 0 but not zero.

[@luzzyrogue]

In your code you have not included 1/r I reckon

[@ayoublamsiah347]

Good evening my teacher, I hope that you are well, good sir I need a pdf file of this court if you like, and I also need resolution in 3D, good night professor.

[@sourabhdas5062]

Can you please share the code, it is not available on the given website link

[@samseleven1]

engineering-stream.com Under Numerical methods

[@vivekbagul2217]

What about surface and centre node

[@samseleven1]

This is FDM not FVM

[@vivekbagul2217]

@@samseleven1 thank you for replying. Can you pls help me with surface and centre node as well

[@samseleven1]

Are you doing finite volume nethod?

[@vivekbagul2217]

@@samseleven1actually I am working on quenching

[@samseleven1]

@@vivekbagul2217 there are boundary nodes and interior nodes.

[@Rubyzzzzzz]

I think the plot should be made by Matlab and the used function name might be ’surf’