ฉันกำลังพยายามสร้างแบบจำลองของการให้ความร้อนแบบไม่คงที่บนจานซึ่งฉันต้องใช้วิธีการเชิงตัวเลขอย่างชัดเจนในการแก้แบบจำลอง ฉันตัดสินใจว่าแผ่นจะทำจากเหล็ก (0,5% C) และฉันพบว่าค่าสัมประสิทธิ์การแพร่กระจายความร้อน $ k / \ rho c_p $ คือ:
$$ \ alpha_ {steel} = 14.74 \ times 10 ^ {- 6} \: \ mathrm {m ^ 2 / s} $$ $$ \ alpha_ {air} = 1.9 \ times 10 ^ {- 5} \: \ mathrm {m ^ 2 / s} $$
ฉันต้องใช้แหล่งความร้อนด้วย แต่ฉันไม่รู้ว่าควรใช้สัมประสิทธิ์อะไร
เนื่องจากเป็นปัญหา 2 มิติรหัสการแก้ไขของฉันจึงเป็นดังนี้:
def heat_source(x, y, matrica):
if (x >= 0.8 and x <= 1.2) and (y >= 0.8 and y <= 1.2):
if matrica[10,10]<300.:
return 280 + randn(1)*20
else:
return 0.
else:
return 0.0
dmx = (M0[iY,iX+1] + M0[iY,iX-1] - 2.0 * M0[iY,iX])/dx**2.0 # conduction in x-direction
dmy = (M0[iY+1,iX] + M0[iY-1,iX] - 2.0 * M0[iY,iX])/dy**2.0 # conduction in y-direction
M_new[iY,iX] = M0[iY, iX] + k_diffusion*dt*(dmx+dmy) + dt*k_heating*heat_source(X[iX], Y[iY], M0) - dt*k_cooling*(M0[iY,iX]-T_ok)
คำถามของฉันคือ:
- ฉันใช้สัมประสิทธิ์ที่ถูกต้องหรือไม่? ถ้าใช่ฉันควรใช้สัมประสิทธิ์อะไรสำหรับแหล่งความร้อน?
- ความแตกต่างระหว่างแหล่งที่ไม่หยุดนิ่งและที่หยุดนิ่งคืออะไร