import numpy as np
import matplotlib.pyplot as plt
#datos
ep=0.8
D=0.001
Re=0.4
Ta=300
h=100
Too=300
I=5.2
c=390
rho=8960
sig=5.67*(10**(-8))
pi=3.14159265359
A=pi*(D**2)/4
#Condiciones iniciales
t=0
Tinc=300
u=Tinc
#Campo de direcciones
def f(u,t):
return ((Re*I**2)/(rho*c*A))-((pi*D*h*(u-Too))/(rho*c*A))- ((pi*D*ep*sig*((u**4)-(Ta**4)))/(rho*c*A))
#Solution
tsol=[t]
Tsol=[u]
dt=0.001
tfin=5
while t<tfin:
k1=f(u,t)
k2=f(u+dt+k1/2,t+dt/2)
k3=f(u+dt*k2/2,t*dt/2)
k4=f(u+dt*k3,t+dt)
u=u+(k1+2*k2+2*k3+k4)*dt/6
t=t+dt
Tsol.append(u)
Tsol.append(t)
plt.plot(tsol,Tsol)
plt.show()
Cuando le doy correr me sale el siguiente error:
ValueError: x and y must have same first dimension, but have shapes (1,) and (10001,)
La ecuación es de transferencia de calor para un alambre delgado, y está basada en el primer ejemplo del siguiente video
¡Gracias!