2D热传导练习

原本是想作为Homework 0,但是时间已经过太久了,而且内容又水就不作为作业丢人了。

求解一个长方区域内的非稳态热传导方程:

video

1.初始温度场的太极形状是借用了Lecture 4的演示代码advection.py中画太极图的函数
2.类的构造是参考了zwang同学的优秀Homework 0,lbm_solver中的写法
3.时间步是显式的,没有求解Ax=b的步骤
4.加了一个点击屏幕的event,每点一下就会把周围的温度提高一下

代码放在了这里,欢迎大佬来指点修改~

8 个赞