Heat transfer matlab code
4 visualizaciones (últimos 30 días)
Mostrar comentarios más antiguos
![](https://www.mathworks.com/matlabcentral/answers/uploaded_files/194226/image.png)
I am trying to write a program to plot the temperature distribution in a insulated rod using the explicit Finite Central Difference Method and 1D Heat equation. The rod is heated on one end at 273.15k and exposed to ambient temperature on the right end at 268.15k. I am using a time of 1s, 11 grid points and a .002s time step. When I plot it gives me a crazy curve which isn't right. I think I am messing up my initial and boundary conditions. Here is my code.
L=1; t=1; k=.001; n=11; nt=500; dx=L/n; dt=.002; alpha=k*dt/dx^2; T0=273.15*ones(1,n); T1=268.15*ones(1,n); T0(1) = 268.15; T0(end) = 268.15; for j=1:nt T1(2:n-1) = T0(2:n-1) + alpha*(T0(3:n)-2*T0(2:n-1)+T0(1:n-2)); T0=T1; end plot(T1)
3 comentarios
KALYAN ACHARJYA
el 22 de Ag. de 2018
Editada: KALYAN ACHARJYA
el 22 de Ag. de 2018
What type of curve do you expect? As I have no idea about Heat Transfer
Kuifeng Zhao
el 23 de Ag. de 2018
If you don't know what you are plotting, how do you know it gives you a curve 'which isn't right'? Then, what is the right thing you want to plot?
Respuestas (0)
Ver también
Categorías
Más información sobre Geometry and Mesh en Help Center y File Exchange.
Community Treasure Hunt
Find the treasures in MATLAB Central and discover how the community can help you!
Start Hunting!