Heat transient equation not executing

2 visualizzazioni (ultimi 30 giorni)
Reji G
Reji G il 17 Set 2022
Commentato: Torsten il 17 Set 2022
I want to implement this formula in matlab. But i couldn't. So it'd be great if anyone can help me. Where x&t ranges from 0.1 to 10.
  2 Commenti
Torsten
Torsten il 17 Set 2022
Please show us your attempt.
Reji G
Reji G il 17 Set 2022
Spostato: Star Strider il 17 Set 2022
clc;
close all;
clear all;
T1=300;
Ti = 350;
L = 10;
alpha = 100;
syms n
T = symsum((exp(((-n*pi)/(2*L))*(alpha*t)))*sin((n*pi*x)/(2*L)),n,1,9)
T(x,t) = T1+((4/pi)*(Ti-T1)*T);

Accedi per commentare.

Risposte (1)

Torsten
Torsten il 17 Set 2022
The terms in your series do not agree with the formula you posted.
You must decide which one is correct and which one is wrong.
clc;
close all;
clear all;
T1 = 300;
Ti = 350;
L = 10;
alpha = 100;
syms n t x
T = symsum((exp(((-n*pi)/(2*L))*(alpha*t)))*sin((n*pi*x)/(2*L)),n,1,20);
T = T1+4/pi*(Ti-T1)*T;
T = matlabFunction(T);
x = 0:0.1:L;
t = 0:0.1:1;
[X,TIME] = meshgrid(x,t);
surf(TIME,X,T(TIME,X))
  2 Commenti
Reji G
Reji G il 17 Set 2022
I need Temperature value(T(x,t)). an x y plot is required. Not the graph you plotted
Torsten
Torsten il 17 Set 2022
clc;
close all;
clear all;
T1 = 300;
Ti = 350;
L = 10;
alpha = 100;
syms n t x
T = symsum((exp(((-n*pi)/(2*L))*(alpha*t)))*sin((n*pi*x)/(2*L)),n,1,20);
T = T1+4/pi*(Ti-T1)*T;
T = matlabFunction(T);
x = 0:0.01:L;
t = 0.1;
figure(1)
plot(x,T(t,x))
x = L/2;
t = 0:0.001:0.5;
figure(2)
plot(t,T(t,x))

Accedi per commentare.

Prodotti


Release

R2019b

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by