Solving the Heat Diffusion Equation (1D PDE) in Python

最新推荐文章于 2024-07-12 21:55:13 发布

哗啦cl

最新推荐文章于 2024-07-12 21:55:13 发布

阅读量383

点赞数

分类专栏： cfd

本文链接：https://blog.csdn.net/qq_43439240/article/details/108908324

版权

cfd 专栏收录该内容

12 篇文章 0 订阅

订阅专栏

Solving the Heat Diffusion Equation (1D PDE) in Python

https://www.youtube.com/watch?v=6-2Wzs0sXd8

import numpy as np
import matplotlib.pyplot as plt


L = 0.1  #Grosor de la pared en metros
n = 10   #Numero de nodos utilizados   
T0 = 20  #Temperatura inicial
T1s = 40 #Temperatura superficial en cara 1
T2s = 20 #Temperatura superficial en cara 2
dx = L/n
alpha = 0.000000054713 #Difusividad termica K/(Rho*C_p) 
t_final = 1800  #Tiempo final en segundos
dt = 60


x = np.linspace(dx/2, L-dx/2, n)

T = np.ones(n)*T0
dTdt = np.empty(n)

t = np.arange(0, t_final, dt)


for j in range(1,len(t)):
    plt.clf()
    for i in range(1, n-1):
        dTdt[i] = alpha*(-(T[i]-T[i-1])/dx**2+(T[i+1]-T[i])/dx**2)
    dTdt[0] = alpha*(-(T[0]-T1s)/dx**2+(T[1]-T[0])/dx**2)
    dTdt[n-1] = alpha*(-(T[n-1]-T[n-2])/dx**2+(T2s-T[n-1])/dx**2)
    T = T + dTdt*dt
    plt.figure(1)
    plt.plot(x,T)
    plt.axis([0, L, 0, 50])
    plt.xlabel('Distance (m)')
    plt.ylabel('Temperature (C)')
    plt.show()
    plt.pause(0.01)