求助 土壤源热泵运行时周围土壤温度场分布的程序

控制方程：+
控制容积法离散
初始条件: 热源温度75 土壤初始温度16
边界条件：土壤远边界温度t(x, y,)=16
热边界
= 2700;
c = 1720;
= 0.4;
= 1.44;
v=1*10^-6;
时间总长 48h
空间总长-3到3m(x,y方向都是) 热源位于原点位置

类似出这种图 在平面内的温度分布